Car−Parrinello Molecular Dynamics Simulation of Liquid Formic Acid

Solubility of CO2 in Aqueous Formic Acid Solutions and the Effect of NaCl Addition: A Molecular Simulation Study

There is a growing interest in the development of routes to produce formic acid from CO₂, such as the electrochemical reduction of CO₂ to formic acid. The solubility of CO₂ in the electrolyte influences the production rate of formic acid. Here, the dependence of the CO₂ solubility in aqueous HCOOH solutions with electrolytes on the composition and the NaCl concentration was studied by Continuous Fractional Component Monte Carlo simulations at 298.15 K and 1 bar. The chemical potentials of CO₂, H₂O, and HCOOH were obtained directly from single simulations, enabling the calculation of Henry coefficients and subsequently considering salting in or salting out effects. As the force fields for HCOOH and H₂O may not be compatible due to the presence of strong hydrogen bonds, the Gibbs-Duhem integration test was used to test this compatibility. The combination of the OPLS/AA force field with a new set of parameters, in combination with the SPC/E force field for water, was selected. It was found that the solubility of CO₂ decreases with increasing NaCl concentration in the solution and increases with the increase of HCOOH concentration. This continues up to a certain concentration of HCOOH in the solution, after which the CO₂ solubility is high and the NaCl concentration has no significant effect.

Bulk Self-Assembly of Giant, Unilamellar Vesicles

The desire to create cell-like models for fundamental science and applications has spurred extensive effort toward creating giant unilamellar vesicles (GUVs). However, a route to selectively self-assemble GUVs in bulk has remained elusive. In bulk solution, membrane-forming molecules such as phospholipids, single-tailed surfactants, and block copolymers typically self-assemble into multilamellar, onion-like structures. So although self-assembly processes can form nanoscale unilamellar vesicles, scaffolding by droplets or surfaces is required to create GUVs. Here we show that it is possible to bulk self-assemble cell-sized GUVs with almost complete selectivity over other vesicle topologies. The seemingly paradoxical pair of features that enables this appears to be having very dynamic molecules at the nanoscale that create unusually rigid membranes. The resultant self-assembly pathway enables encapsulation of molecules and colloids and can also generate model primitive cells that can grow and divide.

Formic acid dimers in a nitrogen matrix

Formic acid (HCOOH) dimers are studied by infrared spectroscopy in a nitrogen matrix and by ab initio calculations. We benefit from the use of a nitrogen matrix where the lifetime of the higher-energy (cis) conformer is very long (∼11 h vs. 7 min in an argon matrix). As a result, in a nitrogen matrix, a large proportion of the cis conformer can be produced by vibrational excitation of the lower-energy (trans) conformer. Three trans-trans, four trans-cis, and three cis-cis dimers are found in the experiments. The spectroscopic information on most of these dimers is enriched compared to the previous studies in an argon matrix. The cis-cis dimers of ordinary formic acid (without deuteration) are reported here for the first time. Several conformational processes are obtained using selective excitation by infrared light, some of them also for the first time. In particular, we report on the formation of cis-cis dimers upon vibrational excitation of trans-cis dimers. Tunneling decays of several dimers have been detected in the dark. The tunneling decay of cis-cis dimers of formic acid as well as the stabilization of cis units in cis-cis dimers is also observed for the first time.

Mapping the Free Energy of Lithium Solvation in the Protic Ionic Liquid Ethylammonuim Nitrate: A Metadynamics Study

Understanding lithium solvation and transport in ionic liquids is important due to their possible application in electrochemical devices. Using first-principles simulations aided by a metadynamics approach we study the free-energy landscape for lithium ions at infinite dilution in ethylammonium nitrate, a protic ionic liquid. We analyze the local structure of the liquid around the lithium cation and obtain a quantitative picture in agreement with experimental findings. Our simulations show that the lowest two free energy minima correspond to conformations with the lithium ion being solvated either by three or four nitrate ions with a transition barrier between them of 0.2 eV. Other less probable conformations having different solvation pattern are also investigated.

Floating orbital molecular dynamics simulations

We introduce an alternative ab initio molecular dynamics simulation as a unification of Hartree-Fock molecular dynamics and the floating orbital approach. The general scheme of the floating orbital molecular dynamics method is presented. Moreover, a simple but sophisticated guess for the orbital centers is provided to reduce the number of electronic structure optimization steps at each molecular dynamics step. The conservation of total energy and angular momentum is investigated in order to validate the floating orbital molecular dynamics approach with and without application of the initial guess. Finally, a water monomer and a water dimer are simulated, and the influence of the orbital floating on certain properties like the dipole moment is investigated.

A theoretical study on the red- and blue-shift hydrogen bonds of cis-trans formic acid dimer in excited states

The excited states of cis-trans formic acid dimer and its monomers have been investigated by time-dependent density functional theory (TDDFT) method. The formation of intermolecular hydrogen bonds O1-H1...O2=C2 and C2-H2...O4=C1 induces bond length lengthening of the groups related to the hydrogen bond, while that of the C2-H2 group is shortened. It is demonstrated that the red-shift hydrogen bond O1-H1...O2=C2 and blue-shift hydrogen bond C2-H2...O4=C1 are both weakened when excited to the S1 state. Moreover, it is found that the groups related to the formation of red-shift hydrogen bond O1-H1...O2=C2 are both strengthened in the S1 state, while the groups related to the blue-shift hydrogen bond C2-H2...O4=C1 are both weakened. This will provide information for the photochemistry and photophysical study of red- and blue-shift hydrogen bond.

Electronic Structures of Formic Acid (HCOOH) and Formate (HCOO(-)) in Aqueous Solutions

The electronic structures of formic acid (HCOOH) and formate (HCOO(-)) have been determined in aqueous solutions over a pH range of 1.88-8.87 using a combination of X-ray photoelectron spectroscopy (XPS), partial electron-yield X-ray absorption spectroscopy (PEY XAS), and density functional theory (DFT). The carbon 1s XPS measurements reveal a binding energy shift of -1.3 eV for deprotonated HCOO(-) compared with neutral HCOOH. Such distinction between neutral HCOOH and deprotonated HCOO(-) cannot be made based solely on the respective carbon K-edge PEY XA spectra. Independent of pH, the C1s → π* state excitations occur at 288.0 eV and may lead to the incorrect conclusion that the energy levels of the π* state are the same for both species. The DFT calculations are consistent with the experimental observations and show a shift to higher energy for both the occupied C1s (lower binding energy) and unoccupied π* orbitals of deprotonated HCOO(-) compared to neutral HCOOH in aqueous solutions.

Car−Parrinello Molecular Dynamics Study of Anharmonic Systems: A Mannich Base in Solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schrödinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.