Water Dimer under Electric Fields: An Ab Initio Investigation up to Quantum Accuracy.

It is well-established that strong electric fields (EFs) can align water dipoles, partially order the H-bond network of liquid water, and induce water splitting and proton transfers. To illuminate the fundamental behavior of water under external EFs, we present the first benchmark, to the best of our knowledge, of DFT calculations of the water dimer exposed to intense EFs against coupled cluster calculations. The analyses of the vibrational Stark effect and electron density provide a consistent picture of the intermolecular charge transfer effects driven along the H-bond by the increasing applied field at all theory levels. However, our findings prove that at extreme field regimes (∼1-2 V/Å) DFT calculations significantly exaggerate by ∼10-30% the field-induced strengthening of the H-bond, both within the GGA, hybrid GGA, and hybrid meta-GGA approximations. Notably, a linear correlation emerges between the vibrational Stark effect on OH stretching and H-bond strengthening: a 1 kcal mol-1 increase corresponds to an 80 cm-1 red-shift in OH stretching frequency.

Electric-field induced entropic effects in liquid water.

Externally applied electric fields in liquid water can induce a plethora of effects with wide implications in electrochemistry and hydrogen-based technologies. Although some effort has been made to elucidate the thermodynamics associated with the application of electric fields in aqueous systems, to the best of our knowledge, field-induced effects on the total and local entropy of bulk water have never been presented so far. Here, we report on classical TIP4P/2005 and ab initio molecular dynamics simulations measuring entropic contributions carried by diverse field intensities in liquid water at room temperature. We find that strong fields are capable of aligning large fractions of molecular dipoles. Nevertheless, the order-maker action of the field leads to quite modest entropy reductions in classical simulations. Albeit more significant variations are recorded during first-principles simulations, the associated entropy modifications are small compared to the entropy change involved in the freezing phenomenon, even at intense fields slightly beneath the molecular dissociation threshold. This finding further corroborates the idea that electrofreezing (i.e., the electric-field-induced crystallization) cannot take place in bulk water at room temperature. In addition, here, we propose a molecular-dynamics-based analysis (3D-2PT) that spatially resolves the local entropy and the number density of bulk water under an electric field, which enables us to map their field-induced changes in the environment of reference H2O molecules. By returning detailed spatial maps of the local order, the proposed approach is capable of establishing a link between entropic and structural modifications with atomistic resolution.

A Computational Quantum-Based Perspective on the Molecular Origins of Life's Building Blocks.

The search for the chemical origins of life represents a long-standing and continuously debated enigma. Despite its exceptional complexity, in the last decades the field has experienced a revival, also owing to the exponential growth of the computing power allowing for efficiently simulating the behavior of matter-including its quantum nature-under disparate conditions found, e.g., on the primordial Earth and on Earth-like planetary systems (i.e., exoplanets). In this minireview, we focus on some advanced computational methods capable of efficiently solving the Schro¨dinger equation at different levels of approximation (i.e., density functional theory)-such as ab initio molecular dynamics-and which are capable to realistically simulate the behavior of matter under the action of energy sources available in prebiotic contexts. In addition, recently developed metadynamics methods coupled with first-principles simulations are here reviewed and exploited to answer to old enigmas and to propose novel scenarios in the exponentially growing research field embedding the study of the chemical origins of life.

Molecular dissociation and proton transfer in aqueous methane solution under an electric field.

Methane-water mixtures are ubiquitous in our solar system and they have been the subject of a wide variety of experimental, theoretical, and computational studies aimed at understanding their behaviour under disparate thermodynamic scenarios, up to extreme planetary ice conditions of pressures and temperatures [Lee and Scandolo, Nat. Commun., 2011, 2, 185]. Although it is well known that electric fields, by interacting with condensed matter, can produce a range of catalytic effects which can be similar to those observed when material systems are pressurised, to the best of our knowledge, no quantum-based computational investigations of methane-water mixtures under an electric field have been reported so far. Here we present a study relying upon state-of-the-art ab initio molecular dynamics simulations where a liquid aqueous methane solution is exposed to strong oriented static and homogeneous electric fields. It turns out that a series of field-induced effects on the dipoles, polarisation, and the electronic structure of both methane and water molecules are recorded. Moreover, upon increasing the field strength, increasing fractions of water molecules are not only re-oriented towards the field direction, but are also dissociated by the field, leading to the release of oxonium and hydroxyde ions in the mixture. However, in contrast to what is observed upon pressurisation (∼50 GPa), where the presence of the water counterions triggers methane ionisation and other reactions, methane molecules preserve their integrity up to the strongest field explored (i.e., 0.50 V Å-1). Interestingly, neither the field-induced molecular dissociation of neat water (i.e., 0.30 V Å-1) nor the proton conductivity typical of pure aqueous samples at these field regimes (i.e., 1.3 S cm-1) are affected by the presence of hydrophobic interactions, at least in a methane-water mixture containing a molar fraction of 40% methane.

Ab initio molecular dynamics simulations and experimental speciation study of levofloxacin under different pH conditions.

Levofloxacin is an extensively employed broad-spectrum antibiotic belonging to the fluoroquinolone class. Despite the extremely wide usage of levofloxacin for a plethora of diseases, the molecular characterization of this antibiotic appears quite poor in the literature. Moreover, the acid-base properties of levofloxacin - crucial for the design of efficient removal techniques from wastewaters - have never extensively been investigated so far. Here we report on a study on the behavior of levofloxacin under standard and diverse pH conditions in liquid water by synergistically employing static quantum-mechanical calculations along with experimental speciation studies. Furthermore, with the aim of characterizing the dynamics of the water solvation shells as well as the protonation and deprotonation mechanisms, here we present the unprecedented quantum-based simulation of levofloxacin in aqueous environments by means of state-of-the-art density-functional-theory-based molecular dynamics. This way, we prove the cooperative role played by the aqueous hydration shells in assisting the proton transfer events and, more importantly, the key place held by the nitrogen atom binding the methyl group of levofloxacin in accepting excess protons eventually present in water. Finally, we also quantify the energetic contribution associated with the presence of a H-bond internal to levofloxacin which, on the one hand, stabilizes the ground-state molecular structure of this antibiotic and, on the other, hinders the first deprotonation step of this fluoroquinolone. Among other things, the synergistic employment of quantum-based calculations and speciation experiments reported here paves the way toward the development of targeted removal approaches of drugs from wastewaters.

Removal of As(III) from Biological Fluids: Mono- versus Dithiolic Ligands.

Arsenic is one of the inorganic pollutants typically found in natural waters, and its toxic effects on the human body are currently of great concern. For this reason, the search for detoxifying agents that can be used in a so-called "chelation therapy" is of primary importance. However, to the aim of finding the thermodynamic behavior of efficient chelating agents, extensive speciation studies, capable of reproducing physiological conditions in terms of pH, temperature, and ionic strength, are in order. Here, we report on the acid-base properties of meso-2,3-dimercaptosuccinic acid (DMSA) at different temperatures (i.e., T = 288.15, 298.15, 310.15, and 318.15 K). In particular, its capability to interact with As(III) has been investigated by experimentally evaluating some crucial thermodynamic parameters (ΔH and TΔS), stability constants, and its speciation model. Additionally, in order to gather information on the microscopic coordination modalities of As(III) with the functional groups of DMSA and, at the same time, to better interpret the experimental results, a series of state-of-the-art ab initio molecular dynamics simulations have been performed. For the sake of completeness, the sequestering capabilities of DMSA-a simple dithiol ligand-toward As(III) are directly compared with those recently emerged from similar analyses reported on monothiol ligands.

Enhanced conductivity of water at the electrified air-water interface: a DFT-MD characterization.

DFT-based molecular dynamics simulations of the electrified air-liquid water interface are presented, where a homogeneous field is applied parallel to the surface plane. We unveil the field intensity for the onset of proton transfer and molecular dissociation; the protonic current/proton conductivity is measured as a function of the field intensity/voltage. The air-water interface is shown to exhibit a proton conductivity twice the one in the liquid water for field intensities below 0.40 V Å-1. We show that this difference arises from the very specific organization of water in the binding interfacial layer (BIL, i.e. the air-water interface region) into a 2D-HBond-network that is maintained and enforced at the electrified interface. Beyond fields of 0.40 V Å-1, water in the BIL and in the bulk liquid are aligned in the same way by the rather intense fields, hence leading to the same proton conductivity in both BIL and bulk water.

Electric-Field-Induced Effects on the Dipole Moment and Vibrational Modes of the Centrosymmetric Indigo Molecule.

Intense static electric fields can strongly perturb chemical bonds and induce frequency shifts of the molecular vibrations in the so-called vibrational Stark effect. Based on a density functional theory (DFT) approach, here, we report a detailed investigation of the influence of oriented external electric fields (OEEFs) on the dipole moment and infrared (IR) spectrum of the nonpolar centrosymmetric indigo molecule. When an OEEF as intense as ∼0.1 V Å-1 is applied, several modifications in the IR spectrum are observed. Besides the notable frequency shift of some modes, we observe the onset of new bands-forbidden by the selection rules in the zero-field case. Such a neat field-induced modification of the vibrational selection rules, and the subsequent variations of the peaks' intensities in the IR spectrum, paves the way toward the design of smart tools employing centrosymmetric molecules as proxies for mapping local electric fields. In fact, here, we show that the ratio between the IR and the Raman intensities of selected modes is proportional to the square of the local field. This indicator can be used to quantitatively measure local fields, not only in condensed matter systems under standard conditions but also in field-emitting-tip apparatus.

Ab Initio Molecular Dynamics Study of Methanol-Water Mixtures under External Electric Fields.

Intense electric fields applied on H-bonded systems are able to induce molecular dissociations, proton transfers, and complex chemical reactions. Nevertheless, the effects induced in heterogeneous molecular systems such as methanol-water mixtures are still elusive. Here we report on a series of state-of-the-art ab initio molecular dynamics simulations of liquid methanol-water mixtures at different molar ratios exposed to static electric fields. If, on the one hand, the presence of water increases the proton conductivity of methanol-water mixtures, on the other, it hinders the typical enhancement of the chemical reactivity induced by electric fields. In particular, a sudden increase of the protonic conductivity is recorded when the amount of water exceeds that of methanol in the mixtures, suggesting that important structural changes of the H-bond network occur. By contrast, the field-induced multifaceted chemistry leading to the synthesis of e.g., hydrogen, dimethyl ether, formaldehyde, and methane observed in neat methanol, in 75:25, and equimolar methanol-water mixtures, completely disappears in samples containing an excess of water and in pure water. The presence of water strongly inhibits the chemical reactivity of methanol.

Ab initio spectroscopy of water under electric fields.

Whereas a broad range of literature exists on the spectroscopy of water in disparate conditions, infrared (IR) and Raman spectra of water subjected to electric fields have never extensively been investigated so far. Based on ab initio molecular dynamics simulations, here we present IR and Raman spectra of bulk liquid water under the effect of static electric fields. A contraction of the entire frequency range is recorded upon increasing the field intensity both in the IR and in the Raman spectra. Whilst the OH stretching band is progressively shifted toward lower frequencies - indicating a field-induced strengthening of the H-bond network - all the other bands are up-shifted by the field. Furthermore, an evident modification of the librational mode band appears in all the spectra. Finally, the order-maker action of the field emerges also from the increase of the water orientational tetrahedral order. Upon field exposure, the water structure becomes more "ice like".

Monte Carlo simulation and integral equation study of Hertzian spheres in the low-temperature regime.

We investigate the behavior of Hertzian spheres in the fluid phase and in proximity of the freezing threshold by using Monte Carlo (MC) simulations and integral equation theories, based on the Ornstein-Zernike (OZ) approach. The study is motivated by the importance of the Hertzian model in representing a large class of systems interacting via soft interactions, such as star polymers or microgels. Radial distribution functions, structure factors, and excess entropy clearly show the reentrant behavior typical of the Hertzian fluid, well caught by both MC simulations and OZ theory. Then, we make use of some phenomenological one-phase criteria for testing their reliability in detecting the freezing threshold. All criteria provide evidence of the fluid-solid transition with different degrees of accuracy. This suggests the possibility to adopt these empirical rules to provide a quick and reasonable estimate of the freezing transition in model potentials of wide interest for soft matter systems.

Mobilities of iodide anions in aqueous solutions for applications in natural dye-sensitized solar cells.

Dye-sensitized solar cells (DSSCs) composed of aqueous electrolytes represent an environmentally friendly, low-cost, and concrete alternative to standard DSSCs and typical solar cells. Although flammable and toxic organic-solvent-based electrolytes have so far been employed more than simpler (iodide) aqueous solutions, recently recorded efficiencies of water-based DSSCs suggest a trend inversion in the near future. Here, we present a study, based on both experiments and ab initio molecular dynamics simulations, in which assessments on the efficiencies of three water electrolytes commonly employed in DSSCs (i.e., LiI, NaI, and KI) are reported. In particular, by atomistically tracing the ability of the iodides as charge carriers and by experimentally measuring the generated currents, we demonstrate that NaI aqueous solutions are more efficient electrolytes than LiI and KI - in descending order - in transporting electrons in DSSCs under bias. Monitoring the role played by the hydration shells of the ionic species under an electric field, we interpret, by first-principles, the various iodide mobilities. This finding, when combined with general considerations on the cation-induced effects on the TiO2 electronic structure, is able to account for the distinct efficiencies of the investigated electrolytes.

Stability of hydrolytic arsenic species in aqueous solutions: As3+vs. As5.

Notwithstanding the fact that arsenic compounds are ubiquitous in the As3+ and As5+ forms in aqueous solutions, most of the microscopic features underlying the conditions of the hydrolysis steps are completely unknown. This way, a first-principles description of the fundamental behaviour of common arsenic species in natural waters and biological fluids is still lacking. Here we report on a synergistic computational and experimental investigation on As3+ and As5+ speciation in aqueous solution under both standard and sizably different alkaline circumstances. If, on the one hand, ab initio molecular dynamics simulations have been used to microscopically trace the different hydrolysis steps of As3+ and As5+ by explicitly taking into account the solvent contribution, on the other hand, they have been able to identify - and predict - the most stable hydrolytic species. In addition, by means of potentiometric and calorimetric measurements, the thermodynamic parameters (log K, ΔH, and TΔS) have been determined at different ionic strength values (0 < I ≤ 1 mol L-1). By comparing the computational and the experimental findings of the species distribution under conditions of some biological fluids, a qualitative agreement on the compounds formed by As3+ and As5+ is thoroughly recorded and, therefore, the stable hydrolytic arsenic species present in natural waters and other biosystems are fully characterised.

Residual Multiparticle Entropy for a Fractal Fluid of Hard Spheres.

The residual multiparticle entropy (RMPE) of a fluid is defined as the difference, Δs, between the excess entropy per particle (relative to an ideal gas with the same temperature and density), sex, and the pair-correlation contribution, s2. Thus, the RMPE represents the net contribution to sex due to spatial correlations involving three, four, or more particles. A heuristic "ordering" criterion identifies the vanishing of the RMPE as an underlying signature of an impending structural or thermodynamic transition of the system from a less ordered to a more spatially organized condition (freezing is a typical example). Regardless of this, the knowledge of the RMPE is important to assess the impact of non-pair multiparticle correlations on the entropy of the fluid. Recently, an accurate and simple proposal for the thermodynamic and structural properties of a hard-sphere fluid in fractional dimension 1

Ionic diffusion and proton transfer in aqueous solutions of alkali metal salts.

We report on a series of ab initio molecular dynamics investigations on LiCl, NaCl, and KCl aqueous solutions under the effect of static electric fields. We have found that although in low-to-moderate field intensity regimes the well-known sequence of cationic mobilities µ(K+) > µ(Na+) > µ(Li+) (i.e., the bigger the cation the higher the mobility) is recovered, from intense field strengths this intuitive rule is no longer verified. In fact, field-induced water molecular dissociations lead to more complex phenomena regulating the standard migration properties of the simplest monovalent cations. The water dissociation threshold is lowered from 0.35 V Å-1 to 0.25 V Å-1 by the presence of charged species in all samples. However, notwithstanding a one-stage process of water ionization and proton conduction takes place at 0.25 V Å-1 in the electrolyte solutions where "structure maker" cations are present (i.e., LiCl and NaCl), the KCl aqueous solution shows some hindrance in establishing a proton conductive regime, which is characterized by the same proton conduction threshold of neat water (i.e., 0.35 V Å-1). In addition, it turns out that protons flow easily in the LiCl (σp = 3.0 S cm-1) solution and then - in descending order - in the NaCl (σp = 2.5 S cm-1) and KCl (σp = 2.3 S cm-1) electrolyte solutions. The protonic conduction efficiency is thus inversely proportional to the ionic radii of the cations present in the samples. Moreover, Cl- anions act as a sort of "protonic well" for high field intensities, further lowering the overall proton transfer efficiency of the aqueous solutions. As a consequence, all the recorded protonic conductivities are lower than that for neat water (σp = 7.8 S cm-1), which strongly indicates that devices exploiting the proton transfer ability should be designed so as to minimize the presence of ionic impurities.

Stability of 2',3' and 3',5' cyclic nucleotides in formamide and in water: a theoretical insight into the factors controlling the accumulation of nucleic acid building blocks in a prebiotic pool.

Synthesis of the first RNAs represents one of the cornerstones of the emergence of life. Recent studies demonstrated powerful scenarios of prebiotic synthesis of cyclic nucleotides in aqueous and formamide environments. This raised a question about their thermodynamic stability, a decisive factor determining their accumulation in a prebiotic pool. Here we performed ab initio molecular dynamics simulations at various temperatures in formamide and water to study the relative stabilities of the 2',3' and 3',5' isomers of cyclic nucleotides. The computations show that in an aqueous environment 2',3' cyclic nucleotides are more stable than their 3',5' counterparts at all temperatures up to the boiling point. In contrast, in formamide higher temperatures favor the accumulation of the 3',5' cyclic form, whereas below about 400 K the 2',3' cyclic form becomes more stable. The latter observation is consistent with a formamide-based origin scenario, suggesting that 3',5' cyclic nucleotides accumulated at higher temperatures subsequently allowed oligomerization reactions after fast cooling to lower temperatures. A statistical analysis of the geometrical parameters of the solutes indicates that thermodynamics of cyclic nucleotides in aqueous and formamide environments are dictated by the floppiness of the molecules rather than by the ring strain of the cyclic phosphodiester linkages.

Virial coefficients, equation of state, and demixing of binary asymmetric nonadditive hard-disk mixtures.

Values of the fifth virial coefficient, compressibility factors, and fluid-fluid coexistence curves of binary asymmetric nonadditive mixtures of hard disks are reported. The former correspond to a wide range of size ratios and positive nonadditivities and have been obtained through a standard Monte Carlo method for the computation of the corresponding cluster integrals. The compressibility factors as functions of density, derived from canonical Monte Carlo simulations, have been obtained for two values of the size ratio (q = 0.4 and q = 0.5), a value of the nonadditivity parameter (Δ = 0.3), and five values of the mole fraction of the species with the biggest diameter (x1 = 0.1, 0.3, 0.5, 0.7, and 0.9). Some points of the coexistence line relative to the fluid-fluid phase transition for the same values of the size ratios and nonadditivity parameter have been obtained from Gibbs ensemble Monte Carlo simulations. A comparison is made between the numerical results and those that follow from some theoretical equations of state.

Miller experiments in atomistic computer simulations.

The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.

Density and structural anomalies in soft-repulsive dimeric fluids.

We report Monte Carlo results for the fluid structure of a system of dimeric particles interacting via a core-softened potential. More specifically, dimers interact through a repulsive pair potential of an inverse-power form, modified in such a way that the repulsion strength is softened for a given range of distances. The aim of such a study is to investigate how both the elongation of the dimers and the softness of the potential affect some features of the model. Our results show that the dimeric fluid exhibits both density and structural anomalies, even if the interaction is not characterized by two length scales. Upon increasing the aspect ratio of the dimers, such anomalies are progressively hindered, with the structural anomaly surviving even after the disappearance of the density anomaly. These results shed light on the peculiar behaviour of molecular systems of non-spherical shape, showing how geometrical and interaction parameters play a fundamental role in determining the presence of anomalies.

Ab initio molecular dynamics study of an aqueous NaCl solution under an electric field.

We report on an ab initio molecular dynamics study of an aqueous NaCl solution under the effect of static electric fields. We found that at low-to-moderate field intensity regimes chlorine ions have a greater mobility than sodium ions which, being a sort of "structure makers", are able to drag their own coordination shells. However, for field strengths exceeding 0.15 V Å(-1) the mobility of sodium ions overcomes that of chlorine ions as both types of ions do actually escape from their respective hydration cages. The presence of charged particles lowers the water dissociation threshold (i.e., the minimum field strength which induces a transfer of protons) from 0.35 V Å(-1) to 0.25 V Å(-1); moreover, a protonic current was also recorded at the estimated dissociation threshold of the solution. The behaviour of the current-voltage diagram of the protonic response to the external electric field is Ohmic as in pure water, with a resulting protonic conductivity of about 2.5 S cm(-1). This value is approximately one third of that estimated in pure water (7.8 S cm(-1)), which shows that the partial breaking of hydrogen bonds induced by the solvated ions hinders the migration of protonic defects. Finally, the conductivity of Na(+) and Cl(-) ions (0.2 S cm(-1)) is in fair agreement with the available experimental data for a solution molarity of 1.7 M.