Conformational effects on circular dichroism in the photoelectron angular distribution.

The B-spline density-functional method has been applied to the conformers of the (1R, 2R)-1,2-dibromo-1,2-dichloro-1,2-difluoroethane molecule. The cross section, asymmetry, and dichroic parameters relative to core and valence orbitals, which do not change their nature along the conformational curve, have been systematically studied. While the cross section and the asymmetry parameter are weakly affected, the dichroic parameter appears to be rather sensitive to the particular conformer of the molecule, suggesting that this dynamical property could be a useful tool for conformational analysis. The computational method has also been applied to methyl rotation in methyloxirane. Unexpected and dramatic sensitivity of the dichroic-parameter profile to the methyl rotation, both in the core and valence states, has been found. Boltzmann averaging over the conformers reproduces quite closely the profiles previously obtained for the minimum-energy conformation, which is in good agreement with the experimental results.     

Emergence of spatio-temporal patterns in forest-fire sequences

Temporal variation of spatial clustering in fire data recorded from 1997 to 2003 in Tuscany region, central Italy, has been investigated using the Voronoï polygon area, the Morishita index and the fractal dimension. Our findings reveal that the spatial clustering of fire events changes with time, showing an enhancement of the clustering degree before the largest events.     

Angle-resolved photoelectron spectroscopy of randomly oriented 3-hydroxytetrahydrofuran enantiomers.

Circular dichroism in the angular distribution of valence photoelectrons emitted from randomly oriented 3-hydroxytetrahydrofuran enantiomers (ThS and ThR) has been observed in gas-phase experiments using circularly polarized vacuum ultraviolet (VUV) light. The measured dichroism for both ThS and ThR, acquired at the single magic angle theta=234.73 degrees and at photon energies of 22, 19, 16, and 14 eV, points to an asymmetric forward-backward scattering of the photoelectrons from their highest occupied molecular orbitals (HOMO) HOMO-1 and HOMO-2, of up to 5%, depending on the photon energy. The asymmetry reverses on exchange of either the helicity of the radiation or the configuration of Th. The photoionization dichroic D parameters of ThS and ThR have been measured and their values discussed in the light of LCAO B-spline density functional theory (DFT) predictions. While an acceptable agreement is found between the dichroic parameter measured and calculated at the highest photon energy for the HOMO and HOMO-2 orbitals of Th, a significant discrepancy is observed for the HOMO-1 state which is attributed to the floppiness of Th, in particular to the comparatively large sensitivity of the size and shape of its HOMO-1 on nuclear vibrational motion.     

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl₂) solutions at different concentrations (0.6–2.8 mol kg⁻¹) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl₂ has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.     

The onset of calcium carbonate nucleation: a density functional theory molecular dynamics and hybrid microsolvation/continuum study

Density functional theory (Perdew-Burke-Ernzerhof) based methods have been used to study the structure and hydration environment of the building blocks of CaCO 3 in aqueous solutions of calcium bicarbonate and calcium carbonate. Car-Parrinello molecular dynamics simulations of Ca(2+)/CO3(2-) and Ca (2+)/HCO3(-) in explicit water were performed to investigate the formation of CaCO3 and the hydration shell of the solvated hetero-ion pair. Our simulations show that the formation of the monomer of CaCO3 occurs with an associative mechanism and that the dominant building block of calcium (bi)carbonate in aqueous solution is Ca[eta(1)-(H)CO3](H2O)5, i.e., the preferred hydration number is five, while the (bi)carbonate is coordinated to the calcium in a monodentate mode. This result agrees with static calculations, where a hybrid approach using a combination of explicit solvent molecules and a polarizable continuum model has been applied to compute the solvation free energies of calcium bicarbonate species. Furthermore, the discrete-continuum calculations predict that the Ca(HCO3)2 and Ca(HCO3)3(-) species are stable in an aqueous environment preferentially as Ca(HCO3)2(H2O)4 and Ca(HCO3)3(H2O)2(-), respectively.     

trans-Fe(II)(H)2(diphosphine)(diamine) complexes as alternative catalysts for the asymmetric hydrogenation of ketones? A DFT study

New insights into the structural, electronic and catalytic properties of Fe complexes are provided by a density functional theory study of model as well as real [Fe(II)(H)(2)(diphosphine)(diamine)] systems. Calculations conducted using several different functionals on the trans- and cis-isomers of [Fe(II)(H)(2)(S-xylbinap)(S,S-dpen)] complexes show that, as with the [Ru(II)(H)(2)(diphosphine)(diamine)] complexes, the trans-[Fe(II)(H)(2)(diphosphine)(diamine)] complex is the more stable isomer. Analysis of the spin states of the trans-[Fe(II)(H)(2)(diphosphine)(diamine)] complexes also shows that the singlet state is significantly more stable than the triplet and the quintet, as with the [Ru(II)(H)(2)(diphosphine)(diamine)] complexes. Calculations of the catalytic cycle for the hydrogenation of ketones using two model trans-[M(II)(H)(2)(PH(3))(2)(en)] catalysts, where M = Ru and Fe, show that the mechanism of reaction as well as the activation energies are very similar, in particular: (i) the ketone/alcohol hydrogen transfer reaction occurs through the metal-ligand bifunctional mechanism, with energy barriers of 3.4 and 3.2 kcal mol(-1) for the Ru- and Fe-catalysed reactions, respectively; (ii) the heterolytic splitting of H(2) across the M[partial double bond, bottom dashed]N bond for the regeneration of the Ru and Fe catalysts has an activation barrier of 13.8 and 12.8 kcal mol(-1), respectively, and is expected to be the rate determining step for both catalytic systems. The reduction of acetophenone by trans-[M(II)(H)(2)(S-xylbinap)(S,S-dpen)] complexes along two competitive reaction pathways, shows that the intermediates for the Fe catalytic system are similar to those responsible for the high enantioselectivity of (R)-alcohol in those proposed trans-[Ru(II)(H)(2)(S-xylbinap)(S,S-dpen)] catalysed acetophenone hydrogenation reaction. Thus the high enantiomeric excess in the hydrogenation of acetophenone could, in principle, be achieved using Fe catalysts.     

A multilateral mechanistic study into asymmetric transfer hydrogenation in water

The mechanism of aqueous-phase asymmetric transfer hydrogenation (ATH) of acetophenone (acp) with HCOONa catalyzed by Ru-TsDPEN has been investigated by stoichiometric reactions, NMR probing, kinetic and isotope effect measurements, DFT modeling, and X-ray structure analysis. The chloride [RuCl(TsDPEN)(p-cymene)] (1), hydride [RuH(TsDPEN)(p-cymene)] (3), and the 16-electorn species [Ru(TsDPEN-H)(p-cymene)] (4) were shown to be involved in the aqueous ATH, with 1 being the precatalyst, and 3 as the active catalyst detectable by NMR in both stoichiometric and catalytic reactions. The formato complex [Ru(OCOH)(TsDPEN)(p-cymene)] (2) was not observed; its existence, however, was demonstrated by its reversible decarboxylation to form 3. Both 1 and 3 were protonated under acidic conditions, leading to ring opening of the TsDPEN ligand. 4 reacted with water, affording a hydroxyl species. In a homogeneous DMF/H(2)O solvent, the ATH was found to be first order in the concentration of catalyst and acp, and inhibited by CO(2). In conjunction with the NMR results, this suggests that hydrogen transfer to ketone is the rate-determining step. The addition of water stabilized the ruthenium catalyst and accelerated the ATH reaction; it does so by participating in the catalytic cycle. DFT calculations revealed that water hydrogen bonds to the ketone oxygen at the transition state of hydrogen transfer, lowering the energy barrier by about 4 kcal mol(-1). The calculations also suggested that the hydrogen transfer is more step-wise in nature rather than concerted. This is supported to some degree by the kinetic isotope effects, which were obscured by extensive H/D scrambling.     

Gel‐Like Structure of a Hexadecyl Derivative of Hyaluronic Acid for the Treatment of Osteoarthritis

Hyaluronic acid is a polysaccharide with viscoelastic and mechanical properties that are crucial for the normal functioning of osteoarticular junctions. It is demonstrated that introduction of a hexadecyl side chain into HA yields an injectable polysaccharide capable of forming physical hydrogels, which are stable at very low polymer concentrations, whereas native hyaluronic acid forms viscous solutions at concentrations that are ten times higher. Characterization of this system showed that the driving force for its gel‐like behavior is the occurrence of hydrophobic interactions involving aliphatic side chains, despite the low degree of substitution, as confirmed by molecular dynamics simulations of HYADD4 and HA hydrogels.

     

Full configuration interaction approach to the few-electron problem in artificial atoms

We present a new high performance configuration interaction code optimally designed for the calculation of the lowest-energy eigenstates of a few electrons in semiconductor quantum dots (also called artificial atoms) in the strong interaction regime. The implementation relies on a single-particle representation, but it is independent of the choice of the single-particle basis and, therefore, of the details of the device and configuration of external fields. Assuming no truncation of the Fock space of Slater determinants generated from the chosen single-particle basis, the code may tackle regimes where Coulomb interaction very effectively mixes many determinants. Typical strongly correlated systems lead to very large diagonalization problems; in our implementation, the secular equation is reduced to its minimal rank by exploiting the symmetry of the effective-mass interacting Hamiltonian, including square total spin. The resulting Hamiltonian is diagonalized via parallel implementation of the Lanczos algorithm. The code gives access to both wave functions and energies of first excited states. Excellent code scalability in a parallel environment is demonstrated; accuracy is tested for the case of up to eight electrons confined in a two-dimensional harmonic trap as the density is progressively diluted up to the Wigner regime, where correlations become dominant. Comparison with previous quantum Monte Carlo simulations in the Wigner regime demonstrates power and flexibility of the method.