Simulating fluid-phase equilibria of water from first principles

Efficient Monte Carlo algorithms and a mixed-basis set electronic structure program were used to compute from first principles the vapor-liquid coexistence curve of water. A water representation based on the Becke-Lee-Yang-Parr exchange and correlation functionals yields a saturated liquid density of 900 kg/m3 at 323 K and normal boiling and critical temperatures of 350 and 550 K, respectively. An analysis of the structural and electronic properties of the saturated liquid phase shows an increase of the asymmetry of the local hydrogen-bonded structure despite the persistence of a 4-fold coordination and decreases of the molecular dipole moment and of the spread of the lowest unoccupied molecular orbital with increasing temperature.     

Harmonic Fourier beads method for studying rare events on rugged energy surfaces

We present a robust, distributable method for computing minimum free energy paths of large molecular systems with rugged energy landscapes. The method, which we call harmonic Fourier beads (HFB), exploits the Fourier representation of a path in an appropriate coordinate space and proceeds iteratively by evolving a discrete set of harmonically restrained path points-beads-to generate positions for the next path. The HFB method does not require explicit knowledge of the free energy to locate the path. To compute the free energy profile along the final path we employ an umbrella sampling method in two generalized dimensions. The proposed HFB method is anticipated to aid the study of rare events in biomolecular systems. Its utility is demonstrated with an application to conformational isomerization of the alanine dipeptide in gas phase.     

Conformation and solvation structure for an isolated n-octadecane chain in water, methanol, and their mixtures

Configurational-bias Monte Carlo simulations in the isobaric-isothermal ensemble (T = 323 K and p = 10 atm) were carried out to probe structural properties of an isolated n-octadecane chain solvated in water, methanol, water-rich, or methanol-rich mixtures and, for comparison, of an isolated chain in the gas phase and for neat liquid n-octadecane. The united-atom version of the TraPPE (transferable potentials for phase equilibria) force field was used to represent n-octadecane and methanol and the TIP-4P model was used for water. In all six environments, broad conformational distributions are observed and the n-octadecane chains are found to predominantly adopt extended, but not all-trans conformations. In addition, a small fraction of more collapsed conformations in which the chain ends approach each other is observed for aqueous hydration, the water-rich solvent mixture and the gas phase, but the simulation data do not support a simple two-state picture with folded and unfolded basins of attraction. For chains in these three "poor" solvent environments, the dihedral angles near the center of the chain show an enhancement of the gauche population. The ensemble of water-solvated chains with end-to-end contacts is preferentially found in a U-shaped conformation rather than a more globular state. An analysis of the local solvation structures in the water-methanol mixtures shows, as expected, an enrichment of the methyl group of methanol near the methylene and methyl segments of the n-octadecane chain. Interestingly, these local bead fractions are enhanced by factors of 2.5 and 1.5 for methyl and methylene segments reflecting the more hydrophobic nature of the former segments.     

Phase behavior of elemental aluminum using monte carlo simulations

Monte Carlo simulations are presented for two models of aluminum: an embedded-atom model and an explicit many-body model. Vapor/liquid coexistence curves are determined using Gibbs ensemble Monte Carlo simulations. The normal boiling points predicted by both models are somewhat higher (by about 10%) than the experimental value. Isothermal constant-stress simulations are used to simulate solid Al from 300 K to the triple point. The solid structures are at least metastable in the face-centered cubic configuration, and the specific heat is determined to be lower than the experimental value. The melting point for the embedded-atom model determined via thermodynamic integration along a pseudo-supercritical path is approximately 20% higher than the experimental value.     

Assessing group-based cutoffs and the Ewald method for electrostatic interactions in clusters and in saturated,superheated, and supersaturated vapor phases of dipolar molecules

Monte Carlo simulations in the canonical, isobaric-isothermal, grand canonical, and Gibbs ensembles were used to assess whether the computationally expensive Ewald summation method for the computation of the first-order electrostatic energy can be replaced with a simpler truncation approach for accurate simulations of the saturated, superheated, and supersaturated vapor phases of dipolar and hydrogen-bonding molecules. Rotationally averaged hydrogen fluoride dimer and trimer energies, thermophysical properties and aggregation in the superheated vapor phase of hydrogen fluoride, nucleation free energy barriers for water, and the vapor–liquid coexistence properties of hydrogen fluoride and water were investigated over a wide range of state points. We find that for densities not too close to the critical density, results obtained from simulations using a spherical potential truncation based on neutral groups (molecules or fragments) for the Coulomb interactions are statistically identical to those obtained using the Ewald summation method. Use of the simpler spherical truncation results in a significant reduction of the computational effort for simulations employing molecular mechanics force fields and also allows for straightforward implementation of many-body expansion methods to compute the potential energy from electronic structure calculations of subsystems of the entire vapor-phase system.     

Liquid structures of water, methanol, and hydrogen fluoride at ambient conditions from first principles molecular dynamics simulations with a dispersion corrected density functional

Using first principles molecular dynamics simulations in the isobaric-isothermal ensemble (T = 300 K, p = 1 atm) with the Becke-Lee-Yang-Parr exchange/correlation functional and a dispersion correction due to Grimme, the hydrogen bonding networks of pure liquid water, methanol, and hydrogen fluoride are probed. Although an accurate density is found for water with this level of electronic structure theory, the average liquid densities for both hydrogen fluoride and methanol are overpredicted by 50 and 25%, respectively. The radial distribution functions indicate somewhat overstructured liquid phases for all three compounds. The number of hydrogen bonds per molecule in water is about twice as high as for methanol and hydrogen fluoride, though the ratio of cohesive energy over number of hydrogen bonds is lower for water. An analysis of the hydrogen-bonded aggregates revealed the presence of mostly linear chains in both hydrogen fluoride and methanol, with a few stable rings and chains spanning the simulation box in the case of hydrogen fluoride. Only an extremely small fraction of smaller clusters was found for water, indicating that its hydrogen bond network is significantly more extensive. A special form of water with on average about two hydrogen bonds per molecule yields a hydrogen-bonding environment significantly different from the other two compounds.     

Palladium Complexes Based on Ylide-Functionalized Phosphines (YPhos): Broadly Applicable High-Performance Precatalysts for the Amination of Aryl Halides at Room Temperature

Palladium allyl, cinnamyl, and indenyl complexes with the ylide-substituted phosphines Cy₃P⁺−C⁻(R)PCy₂ (with R=Me ( L1 ) or Ph ( L2 )) and Cy₃P⁺−C⁻(Me)PtBu₂ ( L3 ) were prepared and applied as defined precatalysts in C−N coupling reactions. The complexes are highly active in the amination of 4-chlorotoluene with a series of different amines. Higher yields were observed with the precatalysts in comparison to the in situ generated catalysts. Changes in the ligand structures allowed for improved selectivities by shutting down β-hydride elimination or diarylation reactions. Particularly, the complexes based on L2 (joYPhos) revealed to be universal precatalysts for various amines and aryl halides. Full conversions to the desired products are reached mostly within 1 h reaction time at room temperature, thus making L2 to one of the most efficient ligands in C−N coupling reactions. The applicability of the catalysts was demonstrated for aryl chlorides, bromides and iodides together with primary and secondary aryl and alkyl amines, including gram-scale applications also with low catalyst loadings of down to 0.05 mol %. Kinetic studies further demonstrated the outstanding activity of the precatalysts with TOF over 10.000 h⁻¹.     

Efficient Pd-Catalyzed Direct Coupling of Aryl Chlorides with Alkyllithium Reagents

Organolithium compounds are amongst the most important organometallic reagents and frequently used in difficult metallation reactions. However, their direct use in the formation of C−C bonds is less established. Although remarkable advances in the coupling of aryllithium compounds have been achieved, Csp²−Csp³ coupling reactions are very limited. Herein, we report the first general protocol for the coupling or aryl chlorides with alkyllithium reagents. Palladium catalysts based on ylide-substituted phosphines (YPhos) were found to be excellently suited for this transformation giving high selectivities at room temperature with a variety of aryl chlorides without the need for an additional transmetallation reagent. This is demonstrated in gram-scale synthesis including building blocks for materials chemistry and pharmaceutical industry. Furthermore, the direct coupling of aryllithiums as well as Grignard reagents with aryl chlorides was also easily accomplished at room temperature.     

Nanoscale reversed-phase liquid chromatography–mass spectrometry of permethylated N-glycans

Reversed-phase liquid chromatography on the nanoscale coupled to electrospray tandem mass spectrometry was used to analyse a mixture of four commercial glycan standards, and the method was further adapted to N-glycans enzymatically released from alpha-1-acid glycoprotein and immunoglobulin gamma. Glycans were permethylated to enable their separation by reversed-phase chromatography and to facilitate interpretation of fragmentation data. Prior to derivatization of glycans by permethylation, they were reduced to cancel anomerism because, although feasible, it was not desired to separate α- and β-anomers. The effect of supplementing chromatographic solvent with sodium hydroxide to guide adduct formation was investigated. Raising the temperature in which the separation was performed improved chromatographic resolution and affected retention times as expected. It was shown by using the tetrasaccharides sialyl Lewis X and sialyl Lewis A that reversed-phase chromatography could achieve the separation of methylated isobaric glycan analytes. Isobaric glycans were detected among the N-glycans of immunoglobulin gamma and further analysed by tandem mass spectrometry.