Car-Parrinello Research Articles

Electrochemical Oxidation of UV Filters: A First-Principles Molecular Dynamics Study.

A theoretical model is proposed to study the oxidation mechanisms of the organic UV filters BP3 and BP4 during electrochemical water treatment utilizing Car-Parrinello molecular dynamics. Factors such as the amount of solvent to be included and how to design the system with the least possible intervention are discussed. The proposed model consist of the optimization of the geometries by density functional theory methods, the equilibration of the structure immersed in a water box, the inclusion of the reactive species, and the analysis of the reaction energies of each reaction pathway. The ab-initio molecular dynamics simulations lead to several products, and some trends can be identified, in accordance with the well-known reactivity rules of organic chemistry. The products proposed in this work are intermediates in longer oxidative pathways.

Does Pb2+ Form Holodirected or Hemidirected Solvation Geometries in Water?

The hydration properties of the Pb2+ ion in aqueous solution have been investigated by using a synergic approach based on Classical and Car-Parrinello molecular dynamics (CPMD) simulations and extended X-ray absorption fine structure (EXAFS) spectroscopy. A definite answer has been given to the main question on the Pb2+ hydration structure, which concerns the formation of either holodirected or hemidirected solvation geometries, such terms referring to the arrangements of the ligands that can be directed either throughout the surface of an encompassing globe around Pb2+ or throughout only part of the globe, respectively. Our CPMD results show that the Pb2+ ion in water forms a hemidirected 4-fold cluster with a well-defined distorted pyramidal geometry. This cluster is directed throughout one side of the first shell globe, while the other side contains either two or three very mobile water molecules that do not form a well-defined geometry around the Pb2+ ion. The Pb2+ first shell structural arrangement determined from the CPMD simulation was confirmed by the EXAFS experimental results. A homodirected Pb2+ first shell complex like the 8-fold SAP structure obtained from the classical MD simulations cannot be reconciled with the EXAFS experimental data. These findings represent a significant step forward in the understanding of the solvation chemistry of the Pb2+ ion, which is fundamental to improving the efficiency of lead removal procedures that are crucial to the safety of water resources.

On the Simulation of Photoreactions Using Restricted Open-Shell Kohn-Sham Theory.

It is a well-established standard to describe ground-state chemical reactions at an ab initio level of multi-electron theory. Fast reactions can be directly simulated. The most widely used approach is density functional theory for the electronic structure in combination with molecular dynamics for the nuclear motion. This approach is known as ab initio molecular dynamics. In contrast, the simulation of excited-state reactions at this level of theory is significantly more difficult. It turns out that the self-consistent solution of the Kohn-Sham equations is not easily reached in excited-state simulations. The first program that solved this problem was the Car-Parrinello molecular dynamics code, using restricted open-shell Kohn-Sham theory. Meanwhile, there are alternatives, most prominently the Q-Chem code, which widens the range of applications. The present study investigates the suitability of both codes for the molecular dynamics simulation of excited-state motion and presents applications to photoreactions.

Synthesis and study of nonlinear optical properties of an enaminone derived from dibenzoylmethane and N,N-diethylaminoaniline

The compound, (Z)-3-((4-(diethylamino)phenyl)amino)-1,3-diphenylprop-2-en-1-one, is synthesized by the reaction of dibenzoylmethane and 4-N,N-diethylaniline. The relative stabilities of the possible tautomers of the molecule are studied via the DFT B3LYP-D3BJ, CAM-B3LYP, M062X, and ωB97XD functionals in conjunction with the 6-311++G(d,p) basis set. The results showed that the enaminone tautomer with the intramolecularly hydrogen bonded chelated ring is the most stable. This is further confirmed by the Car–Parrinello MD calculations in the gas phase as well as the NCI analysis. The electronic spectrum is calculated by the TD DFT B3LYP/c-pVDZ level in ethanol, and the hole–electron analysis is carried out for the interpretation of the bands, which revealed that the longest one at 430 nm is of charge transfer origin while the others are of local transition origin. Atoms-in-molecules calculations in several media and levels of theory predicted that the ρBCP at the hydrogen bond in the gas phase to be 0.03791–0.04255 e/a0 3 which is a characteristic of a medium strong hydrogen bond. Researchers investigated the enaminone’s nonlinear optical (NLO) characteristics when it was exposed to a low power (<1 Watt), single fundamental transverse mode laser beam at 473 nm. By using diffraction patterns (DPs) and Z-scan methods, we calculated the nonlinear refractive index (NLRI) of the enaminone up to 4.597 × 10−11 m2 W−1 using DPs. The resulting DPs are numerically investigated using the Fraunhofer (F.) approximation and the Fresnel-Kirchhoff (F.K.) diffraction integral, showing excellent agreement with experimental findings. We successfully explored all-optical switching (AOS) in enaminone using two laser beams.

Very Strong Hydrogen Bond in Nitrophthalic Cocrystals.

This work presents the studies of a very strong hydrogen bond (VSHB) in biologically active phthalic acids. Research on VSHB comes topical due to its participation in many biological processes. The studies cover the modelling of intermolecular interactions and phthalic acids with 2,4,6-collidine and N,N-dimethyl-4-pyridinamine complexes with aim to obtain a VSHB. The four synthesized complexes were studied by experimental X-ray, IR, and Raman methods, as well as theoretical Car-Parrinello Molecular Dynamics (CP-MD) and Density Functional Theory (DFT) simulations. By variation of the steric repulsion and basicity of the complex' components, a very short intramolecular hydrogen bond was achieved. The potential energy curves calculated by the DFT method were characterized by a low barrier (0.7 and 0.9 kcal/mol) on proton transfer in the OHN intermolecular hydrogen bond for 3-nitrophthalic acid with either 2,4,6-collidine or N,N-dimethyl-4-pyridinamine cocrystals. Moreover, the CP-MD simulations exposed very strong bridging proton dynamics in the intermolecular hydrogen bonds. The accomplished crystallographic and spectroscopic studies indicate that the OHO intramolecular hydrogen bond in 4-nitrophthalic cocrystals is VSHB. The influence of a strong steric effect on the geometry of the studied cocrystals and the stretching vibration bands of the carboxyl and carboxylate groups was elaborated.

Ab Initio Molecular Dynamics Investigation on the Permeation of Sodium and Chloride Ions Through Nanopores in Graphene and Hexagonal Boron Nitride Membranes.

Nanoporous membranes promise energy-efficient water desalination. Hexagonal boron nitride (h-BN), like graphene, exhibits outstanding physical and chemical properties, making it a promising candidate for water treatment. We employed Car-Parrinello molecular dynamics simulations to establish an accurate modeling of Na+ and Cl- permeation through hydrogen passivated nanopores in graphene and h-BN membranes. We demonstrate that ion separation works well for the h-BN system by imposing a barrier of 0.13 eV and 0.24 eV for Na+ and Cl- permeation, respectively. In contrast, for permeation of the graphene nanopore, the Cl- ion faces a minimum of energy of 0.68 eV in the nanopore plane and is prone toward blockade of the nanopore, while the Na+ ion experiences a slight minimum of 0.03 eV. Overall, the desalination performance of h-BN nanopores surpasses that of their graphene counterparts.

Renewed spectroscopic and theoretical research of hydrogen bonding in ascorbic acid

The studies of two isomers of ascorbic acid and their deuteroanalogues, presented in the paper, have been accomplished by vibrational spectroscopy methods and quantum-chemical simulations. The spectroscopic research of L-ascorbic and D-isoascorbic acids have been carried out by the infrared (IR) and Raman (R) techniques. On the basis of the obtained results the spectral interpretation of the hydrogen bonded groups of ascorbic acids has been performed. Car-Parrinello Molecular Dynamics (CPMD) and Density Functional Theory (DFT) have been employed to support spectroscopic experimental findings and shed light onto the bridged proton dynamics in the L- and D- isomers of ascorbic acids. The accurate assignments of the hydrogen bond modes have been accomplished with the application of deuterosubstitution, CPMD-solid state simulations and Potential Energy Distribution (PED) analysis. The spectral and structural results have shown that dependency ν(OH) = f(γ(OH)) is the most common for the OHO hydrogen bond, whereas dependency d(OO) = f(γ(OH)) differs as for the ionic and resonance assisted hydrogen bonds.

A holistic understanding of optical properties in amorphous H-terminated Si-nanostructures: Combining TD-DFT with AIMD

Silicon, traditionally known as an indirect band gap semiconductor, unveils intriguing properties at the nanoscale, stemming from deviations from k-conservation rules within nanostructures. In our study, we scrutinized four hydrogenated Si 0D-nanostructures—Si10H16, Si14H20, Si18H24, and Si22H28—to unravel their dynamic stability under thermal fluctuations and optical characteristics. We initiated our exploration by employing the TD-DFT framework to generate and analyze the optical properties of these geometrically optimized nanostructures. Simultaneously, we conducted ab initio molecular dynamics simulations to examine the structural robustness and thermal stability of the four structures. Leveraging the Car-Parrinello molecular dynamics approach within the Quantum ESPRESSO open software suite, we observed temperature evolution and stability differences among the nanostructures at targeted temperatures 40 and 300 K. Our subsequent investigation delved into the Turbo-Lanczos time-dependent DFT method, unraveling the optical properties and excited-state dynamics of hydrogenated Si nanostructures. The results unveiled shifts towards higher energy absorption edges E0, accompanied by alterations in the permittivity tensor, complex refractive index, oscillator strength, and reflectivity. Notably, the analysis revealed an enlarged HOMO-LUMO gap, distinctive from bulk Si. Furthermore, our models predicted the elimination of phase-dependent E1/E2 optical transition peaks in the imaginary part of the dielectric function, and a gradual decrease in the low-frequency dielectric response with increased hydrogenation of the amorphous structures. These findings underscore the promising applications of hydrogenating Si nanostructures in diverse technological domains such as optoelectronics, memristors, sensors, and quantum computing. Their tunable optical properties, size-dependent behaviors, and compatibility with existing silicon-based devices make them particularly appealing for next-generation technologies.

Dynamic and relativistic effects on Pt-Pt indirect spin-spin coupling in aqueous solution studied by abinitio molecular dynamics and two- vs four-component density functional NMR calculations.

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent-solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin-orbit Kohn-Sham (2c-ZKS) and four-component Dirac-Kohn-Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1-Pt2(NH3)4(Am)2-L2]n+ (Am = α-pyrrolidonate and pivalamidate; L = H2O, Cl-, and Br-) are studied. The computed Pt-Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N-C-O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X-) and aquahalo (L1 = H2O; L2 = X-) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Binding of Carbon Monoxide to Hemoglobin in an Oxygen Environment: Force Field Development for Molecular Dynamics.

Carbon monoxide (CO) is a byproduct of the incomplete combustion of carbon-based fuels, such as wood, coal, gasoline, or natural gas. As incomplete combustion in a fire accident or in an engine, massively produced CO leads to a serious life threat because CO competes with oxygen (O2) binding to hemoglobin and makes people suffer from hypoxia. Although there is hyperbaric O2 therapy for patients with CO poisoning, the nanoscale mechanism of CO dissociation in the O2-rich environment is not completely understood. In this study, we construct the classical force field parameters compatible with the CHARMM for simulating the coordination interactions between hemoglobin, CO, and O2, and use the force field to reveal the impact of O2 on the binding strength between hemoglobin and CO. Density functional theory and Car-Parrinello molecular dynamics simulations are used to obtain the bond energy and equilibrium geometry, and we used machine learning enabled via a feedforward neural network model to obtain the classical force field parameters. We used steered molecular dynamics simulations with a force field to characterize the mechanical strength of the hemoglobin-CO bond before rupture under different simulated O2-rich environments. The results show that as O2 approaches the Fe2+ of heme at a distance smaller than ∼2.8 Å, the coordination bond between CO and Fe2+ is reduced to 50% bond strength in terms of the peak force observed in the rupture process. This weakening effect is also shown by the free energy landscape measured by our metadynamics simulation. Our work suggests that the O2-rich environment around the hemoglobin-CO bond effectively weakens the bonding, so that designing of O2 delivery vector to the site is helpful for alleviating CO binding, which may shed light on de novo drug design for CO poisoning.

Metal-Determined Explosive Characteristics of M(NO3)2(1-AT)x of Thermal and Laser Ignition (M2+ = Cr2+, Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+; x = 2 or 3).

Optical and thermal ignition are two common pathways to initiate the explosion of primary explosives, where laser ignition is a more reliable and safer initiation method. Caused by the current-applied laser igniter with the wavelength of 1064 or 915 nm, the energetic complexes with strong absorption in the near-infrared (NIR) region are possibly applied as laser-ignited explosives. Recently, [Cu(NO3)2(1-AT)3] complex has been synthesized with excellent NIR absorption properties, where 1-amino-5H-tetrazole (1-AT) has been proved to be a promising laser-ignited energetic ligand. To confirm the structure-thermal/optical explosive characteristics, based on the structure of synthesized [Cu(NO3)2(1-AT)3], the commonly used transition-metal cations (M2+ = Cr2+, Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+) have been selected to construct the series of complexes of M(NO3)2(1-AT)x (x = 2 or 3) theoretically. Car-Parrinello molecular dynamics (CPMD) method has been applied to unveil the role of center metals in the initiation and growth pathways. Time-dependent density functional theory (TD-DFT) method is used to explore their charge-transfer (CT) characteristics. The optical characteristic of the metal complex is mainly determined by the behaviors of the 3d electrons of center metals in excitation, where the activity of β-d electrons is an important factor to affect the NIR characteristic of complexes.

13 C NMR as an analytical tool for the detection of carbonic acid and pKa determination.

NMR spectroscopy has become a standard technique in studies both on carbon capture and storage. 13 C NMR allows the detection of two peaks for carbonated aqueous samples: one for CO2(aq) and another one for the species H2 CO3 , HCO3 - , and CO3 2- -herein collectively named Hx CO3 x-2 . The chemical shift of this second peak depends on the molar fraction of the three species in equilibrium and has been used to assess the equilibrium between HCO3 - and CO3 2- . The detection of H2 CO3 at low pH solutions is hindered, because of the concurrent liberation of CO2 when the medium is acidified. Herein, a valved NMR tube facilitates the detection of the Hx CO3 x-2 peak across a wide pH range, even at pH1.8 where the dominant species is H2 CO3 . The method employed the formation of frozen layers of NaH13 CO3 and acid solutions within the tube, which are mixed as the tube reaches room temperature. At this point, the tube is already securely sealed, preventing any loss of CO2 to the atmosphere. A spectrophotometry approach allowed the measurement of the actual pH inside the pressurized NMR tube. The chemical shift for H2 CO3 was determined as 160.33 ± 0.03 ppm, which is in good agreement with value obtained by DFT calculations combined with Car-Parrinello molecular dynamics. The H2 CO3 pKa value determined by the present method was 3.41 ± 0.03, for 15% D2 O aqueous medium and 0.8mol/L ionic strength. The proposed method can be extended to studies about analogs such as alkyl carbonic and carbamic acids.

An investigation of the 2D Rb2O monolayer intrinsic properties from a theoretical perspective using DFT

In this research, Density Functional Theory (DFT) was employed to investigate the 1T phase of the Rb2O monolayer. The structural, electronic, optical, and thermal properties were studied utilizing the WIEN2k software. The band gaps found using GGA, LDA, mBJ, and hybrid correlation were 0.52 eV, 0.49 eV, 1.52 eV, and 1.17 eV. The phonon dispersion curve is implemented to assess the Rb2O monolayer's dynamical stability. The thermal stability is assessed using the Car-Parrinello Molecular Dynamics (CPMD) code of Quantum espresso. The absorption coefficient attains its peak value of 35.74 cm−1 at 9.72 eV. Based on the optical property results, the Rb2O monolayer shows promise in designing photovoltaic devices. At 400K, the Figure of merit found to be 0.85 which indicates that the Rb2O monolayer can be used for waste heat recovery systems. The excitonic effects are ascertained by the effective masses of the holes and electrons.

Double-proton transfer triggering mechanism of nitromethane caused by intermolecular hydrogen bonding at low temperature

The initial triggering mechanism of nitromethane is investigated by first-principles molecular dynamics simulations and temperature-dependent Raman spectroscopy experiment. The Raman spectra at different temperature are obtained by Car-Parrinello molecular dynamics simulation to explore the activated vibration modes caused by low-temperature heat. The temperature-dependent Raman spectra experiment is constructed to compare with the theoretical data. Each possible initial reaction path caused by the activated vibration mode is analyzed to reveal the triggering mechanism of the initial chemical bond breaking for nitromethane. The results show that the intermolecular double-proton transfer is the first reaction step which is triggered by the intermolecular hydrogen bond reinforced with the active C–H stretching vibration (near 3000 cm−1), which is not reported before. The scanning energy barrier related to the increase of C–H bond lengths is about 60 kcal/mol, which is consistent with the calculated activation barrier (59.8 kcal/mol) of intermolecular proton transfer reaction. This proves the reliability of our conclusions. The results provide a new interpretation for the long-debated initial decomposition mechanism of nitromethane.

Electric field-driven up-and-down motion of the flexible tail of Al13+ cluster system-a nano-scale flipper.

Dynamic metal nanoclusters have become a hot area of research in the field of nanoscience and nanotechnology due to their potential applications in micro devices. One such dynamic cluster is a quasi-planar ground state (GS) Al13+ cluster which exhibits anelectric field driven up and down flipping motion of the flexible tail which oscillates with respect to the mean plane. A Car-Parrinello molecular dynamics (CPMD) simulation has been carried out to understand the nature of dynamics of the cluster. CPMD simulation study reveals that the flexible tail region of the Al13+ isomeric system (two ground states M1, M2 and a transition state TS connecting them) can be engaged in a systematic up down flipping motion by the application of a transverse electric field. A saw tooth electric field of amplitude 5.19V/nm is sufficient to induce the up-and-down flipping oscillation of the cluster, which has an average oscillation frequency of around 20 THz. AIM, NICS and AdNDP analyses also have been carried out to understand the fluxional nature of the cluster from the electronic structural perspective. Electronic structural analysis of selected optimized intermediate states in the presence of transverse electric field has also been analyzed to correlate the electronic structure with the dynamic nature of the cluster. Single-point energies of all intermediate states between two minima of Al13+ clusters connected through a transition state cluster. Optimized geometries of Al13+ clusters in the presence of electric field of different strengths have been carried out by using the Gaussian 03 package. 6-311 + G(d) basis set and B3LYP hybrid density functional have been utilized for these studies. To establish the flipping motion, Car-Parrinello molecular dynamics (CPMD) has been performed using the cp.x module of the Quantum ESPRESSO 6.3.0 program package using the Perdew-Burke-Ernzerhof (PBE) functional, plane-wave basis set and ultrasoft pseudopotentials. ORTEP-3 and POV ray-3.7 software packages have been used for visualization and graphics generation. Atoms in molecule (AIM), Adaptive Natural Density Partitioning (AdNDP) analysis have been carried out using Multiwfn 3.7 program package.

Charge carrier dynamics in conducting polymer PEDOT using abinitio molecular dynamics simulations.

As conducting polymers become increasingly important in electronic devices, understanding their charge transport is essential for material and device development. Various semi-empirical approaches have been used to describe temporal charge carrier dynamics in these materials, but there have yet to be any theoretical approaches utilizing abinitio molecular dynamics. In this work, we develop a computational technique based on abinitio Car-Parrinello molecular dynamics to trace charge carrier temporal motion in archetypical conducting polymer poly(3,4-ethylenedioxythiophene) (PEDOT). Particularly, we analyze charge dynamics in a single PEDOT chain and in two coupled chains with different degrees of coupling and study the effect of temperature. In our model we first initiate a positively charged polaron (compensated by a negative counterion) at one end of the chain, and subsequently displace the counterion to the other end of the chain and trace polaron dynamics in the system by monitoring bond length alternation in the PEDOT backbone and charge density distribution. We find that at low temperature (T = 1K) the polaron distortion gradually disappears from its initial location and reappears near the new position of the counterion. At the room temperature (T = 300K), we find that the distortions induced by polaron, and atomic vibrations are of the same magnitude, which makes tracking the polaron distortion challenging because it is hidden behind the temperature-induced vibrations. The novel approach developed in this work can be used to study polaron mobility along and between the chains, investigate charge transport in highly doped polymers, and explore other flexible polymers, including n-doped ones.

Enabling High-Performance All-Solid-State Batteries via Guest Wrench in Zeolite Strategy.

All-solid-state batteries with a high energy density and safety are desirable candidates for next-generation energy storage applications. However, conventional solid electrolytes for all-solid-state batteries encounter limitations such as poor ionic conduction, interfacial compatibility, instability, and high cost. Herein, taking advantage of the ingenious capability of zeolite to incorporate functional guests in its void space, we present an innovative ionic activation strategy based on the "guest wrench" mechanism, by introducing a pair of cation and anion of LiTFSI-based guest species (GS) into the supercage of the LiX zeolite, to fabricate a zeolite membrane (ZM)-based solid electrolyte (GS-ZM) with high Li ionic conduction and interfacial compatibility. The restriction of zeolite frameworks toward the framework-associated Li ions is significantly reduced through the dynamic coordination of Li ions with the "oxygen wrench" of TFSI- at room temperature as shown by experiments and Car-Parrinello molecular dynamics simulations. Consequently, the GS-ZM shows an ∼100% increase in ionic conductivity compared with ZM and an outstanding Li+ transference number of 0.97. Remarkably, leveraging the superior ionic conduction of GS-ZM with the favorable interface structure between GS-ZM and electrodes, the assembled all-solid-state Li-ion and Li-air batteries based on GS-ZM exhibit the best-level electrochemical performance much superior to batteries based on liquid electrolytes: a capacity retention of 99.3% after 800 cycles at 1 C for all-solid-state Li-ion batteries and a cycle life of 909 cycles at 500 mA g-1 for all-solid-state Li-air batteries. The mechanistic discovery of a "guest wrench" in zeolite will significantly enhance the adaptability of zeolite-based electrolytes in a variety of all-solid-state energy storage systems with high performance, high safety, and low cost.

Ab Initio Study of Glycine Formation in the Condensed Phase: Carbon Monoxide, Formaldimine, and Water Are Enough

Glycine is considered to be crucial in the formation of proteins and prebiotic substances. Nevertheless, the mechanism of spontaneous glycine formation under prebiotic Earth conditions or within the interstellar medium (ISM) remains a topic of debate, given the changing geochemical environment over Earth’s history and the difficulty of detecting it within the ISM. Yet it is believed that its formation could be possible in interstellar water-rich ice. In this study, using molecular dynamics (MD) simulations at the ab initio level of theory enhanced with modern free energy calculations, we modeled the chemical reaction between carbon monoxide, formaldimine, and water to produce glycine. We estimated under what conditions, in condensed phase at 50, 70, 100, and 300 K, glycine is formed. We also explored the effect of different electric fields on this process. Our results show that glycine could be formed with energy barriers as low as 0.5 kcal mol−1 at 50 K. We discuss whether this reaction could be a suitable candidate for explaining the mechanism of glycine formation under conditions that resemble various astrophysical environments, such as planets, exoplanets, and Earth. This study is relevant to finding a consensus among various proposals for glycine formation. Moreover, it highlights the importance of metadynamics and Car–Parrinello MD methods as tools in finding unknown complex, multistep reaction mechanism pathways, possibly important to the astronomical phenomena.

Competitive adsorption between sodium citrate and naphthenic acids on alumina surfaces: Experimental and computational study

The adsorption of naphthenic acids (NAs) on clay surfaces induces an increase in surface hydrophobicity, which is detrimental to crude oil recovery and froth quality. Sodium citrate (Na3Cit), with more functional groups than NAs, is expected to compete with NAs for adsorption on mineral surfaces. The adsorption of NAs on alumina surfaces was quantified in the presence and absence of Na3Cit at pH 8 ± 0.1. An increase in hydrophobicity of the alumina surface was observed after the adsorption of NAs without Na3Cit. However, adding Na3Cit significantly reduced NAs adsorption and protected the surface wettability. The viscoelasticity and electrokinetic properties of adsorbed materials indicate Na3Cit was preferentially covered on alumina surfaces and repelled NAs by its negative charges. The competition between Na3Cit and NAs for adsorption sites was explained using ab initio Car-Parrinello molecular dynamics (CPMD) simulation and periodic plane-wave density functional theory (DFT) to calculate inner-sphere adsorption and outer-sphere adsorption, respectively. Favourable adsorption of sodium citrate than NAs representatives on aluminum (hydr)oxide cluster and the basal and edge planes of gibbsite were observed from the perspective of free energy barrier and adsorption energy, respectively.

Unraveling the Ag+ ion coordination and solvation thermodynamics in the 1-butyl-3-methylimidazolium tetrafluoroborate ionic liquid

The solvation of the Ag+ ion in the 1-butyl-3-methylimidazolium tetrafluoroborate ([C4mim][BF4]) ionic liquid (IL) has been studied by means of experimental and theoretical methods with the aim of elucidating the cation coordination structure and thermodynamic properties. Car-Parrinello molecular dynamics (CPMD) simulations showed that the Ag+ ion is coordinated by an average number of four [BF4]− anions in a pseudo-tetrahedral geometry. A high configurational disorder of the first solvation sphere is found, where the anions can be found both in mono- and bidentate coordination mode around the Ag+ ion. Also, a solvational equilibrium is observed as due to [BF4]− anion dissociation along the trajectory. The analysis of X-ray absorption spectroscopy data confirmed the picture provided by the CPMD simulation. Classical molecular dynamics simulations were carried out to obtain the single-ion solvation thermodynamic parameters. The negative water → IL free energy of transfer suggests that the Ag+ ion is more favorably solvated in the [C4mim][BF4] IL than in water. This behavior is due to a balance between the enthalpic and entropic contributions, which allows to find a rationale to the strong solvation capabilities of BF4-based ILs towards Ag+.