Inorganic Polysulfides in Solution: Structural Properties and Conformational Isomerism.

We present a comprehensive theoretical examination of the structural properties of dianionic polysulfides [Sn]2- (n = 2-6), their conjugated monoacids [HSn]- (n = 2-6), and a selection of 1e--oxidized radical anions [Sn]•- (n = 2-4), in aqueous and dimethyl sulfoxide (DMSO) solutions. We investigated the structures and stabilities of various conformational isomers within these families of compounds by employing Quantum Mechanics-Molecular Mechanics (QM-MM) Molecular Dynamics (MD) simulations. The explicit inclusion of solvent molecules in the calculations revealed stable conformational structures that were previously unreported and might have appreciable concentrations in real systems. The interconversions between the isomeric structures proceed on the order of hundreds of picoseconds and are energetically similar to the isomerization processes in substituted cyclohexanes. We also conducted a detailed analysis of the stability of different isomers of the radical anion [S4]•- in solution. Our findings highlight the significant influence of the solvent on the isomerizations, a result that could be particularly relevant for enhancing the performance of metal-sulfur batteries.

Modeling the electroluminescence of atomic wires from quantum dynamics simulations.

Static and time-dependent quantum-mechanical approaches have been employed in the literature to characterize the physics of light-emitting molecules and nanostructures. However, the electromagnetic emission induced by an input current has remained beyond the realm of molecular simulations. This is the challenge addressed here with the help of an equation of motion for the density matrix coupled to a photon bath based on a Redfield formulation. This equation is evolved within the framework of the driven-Liouville von Neumann approach, which incorporates open boundaries by introducing an applied bias and a circulating current. The dissipated electromagnetic power can be computed in this context from the time derivative of the energy. This scheme is applied in combination with a self-consistent tight-binding Hamiltonian to investigate the effects of bias and molecular size on the electroluminescence of metallic and semiconducting chains. For the latter, a complex interplay between bias and molecular length is observed: there is an optimal number of atoms that maximizes the emitted power at high voltages but not at low ones. This unanticipated behavior can be understood in terms of the band bending produced along the semiconducting chain, a phenomenon that is captured by the self-consistency of the method. A simple analytical model is proposed that explains the main features revealed by the simulations. The methodology, applied here at a self-consistent tight-binding level but extendable to more sophisticated Hamiltonians such as density functional tight binding and time dependent density functional theory, promises to be helpful for quantifying the power and quantum efficiency of nanoscale electroluminescent devices.

Nanobubble Stability and Formation on Solid-Liquid Interfaces in Open Environments.

Are surface nanobubbles transient or thermodynamically stable structures? This question remained controversial until recently, when the stability of gas nanobubbles at solid-liquid interfaces was demonstrated from thermodynamic arguments in closed systems, establishing that bubbles with radii of hundreds of nanometers can be stable at modest supersaturations if the gas amount is finite. Here we develop a grand-canonical description of bubble formation that predicts that nanobubbles can nucleate and remain thermodynamically stable in open boundaries at high supersaturations when pinned to hydrophobic supports as small as a few nanometers. While larger bubbles can also be stable at lower supersaturations, the corresponding barriers are orders of magnitude above kT, meaning that their formation cannot proceed via heterogeneous nucleation on a uniform solid interface but must follow some alternative path. Moreover, we conclude that a source of growth-limiting mechanism, such as pinning or gas availability, is necessary for the thermodynamic stabilization of surface bubbles.

Fluorescence in quantum dynamics: Accurate spectra require post-mean-field approaches.

Real time modeling of fluorescence with vibronic resolution entails the representation of the light-matter interaction coupled to a quantum-mechanical description of the phonons and is therefore a challenging problem. In this work, taking advantage of the difference in timescales characterizing internal conversion and radiative relaxation-which allows us to decouple these two phenomena by sequentially modeling one after the other-we simulate the electron dynamics of fluorescence through a master equation derived from the Redfield formalism. Moreover, we explore the use of a recent semiclassical dissipative equation of motion [C. M. Bustamante et al., Phys. Rev. Lett. 126, 087401 (2021)], termed coherent electron electric-field dynamics (CEED), to describe the radiative stage. By comparing the results with those from the full quantum-electrodynamics treatment, we find that the semiclassical model does not reproduce the right amplitudes in the emission spectra when the radiative process involves the de-excitation to a manifold of closely lying states. We argue that this flaw is inherent to any mean-field approach and is the case with CEED. This effect is critical for the study of light-matter interaction, and this work is, to our knowledge, the first one to report this problem. We note that CEED reproduces the correct frequencies in agreement with quantum electrodynamics. This is a major asset of the semiclassical model, since the emission peak positions will be predicted correctly without any prior assumption about the nature of the molecular Hamiltonian. This is not so for the quantum electrodynamics approach, where access to the spectral information relies on knowledge of the Hamiltonian eigenvalues.

Water Self-Dissociation is Insensitive to Nanoscale Environments.

Nanoconfinement effects on water dissociation and reactivity remain controversial, despite their importance to understand the aqueous chemistry at interfaces, pores, or aerosols. The pKw in confined environments has been assessed from experiments and simulations in a few specific cases, leading to dissimilar conclusions. Here, with the use of carefully designed ab initio simulations, we demonstrate that the energetics of bulk water dissociation is conserved intact to unexpectedly small length-scales, down to aggregates of only a dozen molecules or pores of widths below 2 nm. The reason is that most of the free-energy involved in water autoionization comes from breaking the O-H covalent bond, which has a comparable barrier in the bulk liquid, in a small droplet of nanometer size, or in a nanopore in the absence of strong interfacial interactions. Thus, dissociation free-energy profiles in nanoscopic aggregates or in 2D slabs of 1 nm width reproduce the behavior corresponding to the bulk liquid, regardless of whether the corresponding nanophase is delimited by a solid or a gas interface. The present work provides a definite and fundamental description of the mechanism and thermodynamics of water dissociation at different scales with broader implications on reactivity and self-ionization at the air-liquid interface.

Dissipative Equation of Motion for Electromagnetic Radiation in Quantum Dynamics.

The dynamical description of the radiative decay of an electronically excited state in realistic many-particle systems is an unresolved challenge. In the present investigation electromagnetic radiation of the charge density is approximated as the power dissipated by a classical dipole, to cast the emission in closed form as a unitary single-electron theory. This results in a formalism of unprecedented efficiency, critical for ab initio modeling, which exhibits at the same time remarkable properties: it quantitatively predicts decay rates, natural broadening, and absorption intensities. Exquisitely accurate excitation lifetimes are obtained from time-dependent DFT simulations for C^{2+}, B^{+}, and Be, of 0.565, 0.831, and 1.97 ns, respectively, in accord with experimental values of 0.57±0.02, 0.86±0.07, and 1.77-2.5 ns. Hence, the present development expands the frontiers of quantum dynamics, bringing within reach first-principles simulations of a wealth of photophysical phenomena, from fluorescence to time-resolved spectroscopies.

Structure and dynamics of nanoconfined water and aqueous solutions.

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Doping and coupling strength in molecular conductors: polyacetylene as a case study.

The doping mechanisms responsible for elevating the currents up to eleven orders of magnitude in semiconducting polymer films are today well characterized. Doping can also improve the performance of nanoscale devices or single molecule conductors, but the mechanism in this case appears to be different, with theoretical studies suggesting that the dopant affects the electronic properties of the junctions. In the present report, multiscale time-dependent DFT transport simulations help clarify the way in which n-type doping can raise the current flowing through a polymer chain connected to a pair of electrodes, with the focus on polyacetylene. In particular, our multiscale methodology offers control over the magnitude of the chemical coupling between the molecule and the electrodes, which allows us to analyze the effect of doping in low and strong coupling regimes. Interestingly, our results establish that the impact of dopants is the highest in weakly coupled devices, while their presence tends to be irrelevant in low-resistance junctions. Our calculations point out that both the equalization of the frontier orbitals with the Fermi level and a small gap between the HOMO and the LUMO must result from doping in order to observe any significant increase of the currents.

Spectroscopy in Complex Environments from QM-MM Simulations.

The applications of multiscale quantum-classical (QM-MM) approaches have shown an extraordinary expansion and diversification in the last couple of decades. A great proportion of these efforts have been devoted to interpreting and reproducing spectroscopic experiments in a variety of complex environments such as solutions, interfaces, and biological systems. Today, QM-MM-based computational spectroscopy methods constitute accomplished tools with refined predictive power. The present review summarizes the advances that have been made in QM-MM approaches to UV-visible, Raman, IR, NMR, electron paramagnetic resonance, and Mössbauer spectroscopies, providing in every case an introductory discussion of the corresponding methodological background. A representative number of applications are presented to illustrate the historical evolution and the state of the art of this field, highlighting the advantages and limitations of the available methodologies. Finally, we present our view of the perspectives and open challenges in the field.

A simple approximation to the electron-phonon interaction in population dynamics.

The modeling of coupled electron-ion dynamics including a quantum description of the nuclear degrees of freedom has remained a costly and technically difficult practice. The kinetic model for electron-phonon interaction provides an efficient approach to this problem, for systems evolving with low amplitude fluctuations, in a quasi-stationary state. In this work, we propose an extension of the kinetic model to include the effect of coherences, which are absent in the original approach. The new scheme, referred to as Liouville-von Neumann + Kinetic Equation (or LvN + KE), is implemented here in the context of a tight-binding Hamiltonian and employed to model the broadening, caused by the nuclear vibrations, of the electronic absorption bands of an atomic wire. The results, which show close agreement with the predictions given by Fermi's golden rule (FGR), serve as a validation of the methodology. Thereafter, the method is applied to the electron-phonon interaction in transport simulations, adopting to this end the driven Liouville-von Neumann equation to model open quantum boundaries. In this case, the LvN + KE model qualitatively captures the Joule heating effect and Ohm's law. It, however, exhibits numerical discrepancies with respect to the results based on FGR, attributable to the fact that the quasi-stationary state is defined taking into consideration the eigenstates of the closed system rather than those of the open boundary system. The simplicity and numerical efficiency of this approach and its ability to capture the essential physics of the electron-phonon coupling make it an attractive route to first-principles electron-ion dynamics.

Mechanisms of Nucleation and Stationary States of Electrochemically Generated Nanobubbles.

Gas evolving reactions are ubiquitous in the operation of electrochemical devices. Recent studies of individual gas bubbles on nanoelectrodes have resulted in unprecedented control and insights on their formation. The experiments, however, lack the spatial resolution to elucidate the molecular pathway of nucleation of nanobubbles and their stationary size and shape. Here we use molecular simulations with an algorithm that mimics the electrochemical formation of gas, to investigate the mechanisms of nucleation of gas bubbles on nanoelectrodes, and characterize their stationary states. The simulations reproduce the experimental currents in the induction and stationary stages, and indicate that surface nanobubbles nucleate through a classical mechanism. We identify three distinct regimes for bubble nucleation, depending on the binding free energy per area of bubble to the electrode, Δγbind. If Δγbind is negative, the nucleation is heterogeneous and the nanobubble remains bound to the electrode, resulting in a low-current stationary state. For very negative Δγ, the bubble fully wets the electrode, forming a one-layer-thick micropancake that nucleates without supersaturation. On the other hand, when Δγbind > 0 the nanobubble nucleates homogeneously close to the electrode, but never attaches to it. We conclude that all surface nanobubbles must nucleate heterogeneously. The simulations reveal that the size and contact angle of stationary nanobubbles increase with the reaction driving force, although their residual current is invariant. The myriad of driven nonequilibrium stationary states with the same rate of production of gas, but distinct bubble properties, suggests that these dissipative systems have attractors that control the stationary current.

Halide-Mediated Modification of Magnetism and Electronic Structure of α-Co(II) Hydroxides: Synthesis, Characterization, and DFT+U Simulations.

The present study introduces a comprehensive exploration in terms of physicochemical characterization and calculations based on density functional theory with Hubbard's correction (DFT+U) of the whole family of α-Co(II) hydroxyhalide (F, Cl, Br, I). These samples were synthesized at room temperature by employing a one-pot approach based on the epoxide route. A thorough characterization (powder X-ray diffraction, X-ray photoelectron spectroscopy, thermogravimetric analysis/mass spectroscopy, and magnetic and conductivity measurements) corroborated by simulation is presented that analyzes the structural, magnetic, and electronic aspects. Beyond the inherent tendency of intercalated anions to modify the interlayer distance, the halide's nature has a marked effect on several aspects. Such as the modulation of the CoOh to CoTd ratio, as well as the inherent tendency towards dehydration and irreversible decomposition. Whereas the magnetic behavior is strongly correlated with the CoTd amount reflected in the presence of glassy behavior with high magnetic disorder, the electrical properties depend mainly on the nature of the halide. The computed electronic structures suggest that the CoTd molar fraction exerts a minor effect on the inherent conductivity of the phases. However, the band gap of the solid turns out to be significantly dependent on the nature of the incorporated halide, governed by ligand to metal charge transfer, which minimizes the gap as the anionic radius becomes larger. Conductivity measurements of pressed pellets confirm this trend. To the best of our knowledge, this is the first report on the magnetic and electrical properties of α-Co(II) hydroxyhalides validated with in silico descriptions, opening the gate for the rational design of layered hydroxylated phases with tunable electrical, optical, and magnetic properties.

Multiscale approach to electron transport dynamics.

Molecular simulations of transport dynamics in nanostructures usually require the implementation of open quantum boundary conditions. This can be instrumented in different frameworks including Green's functions, absorbing potentials, or the driven Liouville von Neumann equation, among others. In any case, the application of these approaches involves the use of large electrodes that introduce a high computational demand when dealing with first-principles calculations. Here, we propose a hybrid scheme where the electrodes are described at a semiempirical, tight binding level, coupled to a molecule or device represented with density functional theory (DFT). This strategy allows us to use massive electrodes at a negligible computational cost, preserving the accuracy of the DFT method in the modeling of the transport properties, provided that the electronic structure of every lead is properly defined to behave as a conducting fermionic reservoir. We study the nature of the multiscale coupling and validate the methodology through the computation of the tunneling decay constant in polyacetylene and of quantum interference effects in an aromatic ring. The present implementation is applied both in microcanonical and grand-canonical frameworks, in the last case using the Driven Liouville von Neumann equation, discussing the advantages of one or the other. Finally, this multiscale scheme is employed to investigate the role of an electric field applied normally to transport in the conductance of polyacetylene. It is shown that the magnitude and the incidence angle of the applied field have a considerable effect on the electron flow, hence constituting an interesting tool for current control in nanocircuits.

Interplay of Coordination Environment and Magnetic Behavior of Layered Co(II) Hydroxichlorides: A DFT+U Study.

In this work we present a systematic computational study of the structural and magnetic properties of a layered family of Co(II) hydroxichlorides, obeying to the general formula Co(OH)2- xCl x(H2O) y. This solid contains both octahedral and tetrahedral cobalt ions, displaying a complex magnetic order arising from the particular coupling between the two kinds of metallic centers. Here, supercells representing concentrations of 12, 20, and 40% of tetrahedral sites were modeled consistently with the compositions reported experimentally. Our simulations show that the two types of cobalt ions tend to couple antiferromagnetically, giving rise to a net magnetic moment slightly out of the plane of the layers. The band gap reaches its minimum value of 1.4 eV for the most diluted fraction of tetrahedral Co(II) sites, going up to 2.2 eV when the content is 40%. Moreover, our results suggest that the presence of interlayer water stabilizes the material and at the same time strongly modifies the electronic environment of tetrahedral Co(II), leading to a further drop of the band gap. To our knowledge, this is the first theoretical investigation of this material.

Stability and Vapor Pressure of Aqueous Aggregates and Aerosols Containing a Monovalent Ion.

The incidence of charged particles on the nucleation and the stability of aqueous aggregates and aerosols was reported more than a century ago. Many studies have been conducted ever since to characterize the stability, structure, and nucleation barrier of ion-water droplets. Most of these studies have focused on the free-energy surface as a function of cluster size, with an emphasis on the role of ionic charge and radius. This knowledge is fundamental to go beyond the rudimentary ion-induced classical nucleation theory. In the present article, we address this problem from a different perspective, by computing the vapor pressures of (H2O)nLi+ and (H2O)nCl- aggregates using molecular simulations. Our calculations shed light on the structure, the critical size, the range of stability, and the role of ion-water interactions in aqueous clusters. Moreover, they allow one to assess the accuracy of the classical thermodynamic model, highlighting its strengths and weaknesses.

Electron transport in real time from first-principles.

While the vast majority of calculations reported on molecular conductance have been based on the static non-equilibrium Green's function formalism combined with density functional theory (DFT), in recent years a few time-dependent approaches to transport have started to emerge. Among these, the driven Liouville-von Neumann equation [C. G. Sánchez et al., J. Chem. Phys. 124, 214708 (2006)] is a simple and appealing route relying on a tunable rate parameter, which has been explored in the context of semi-empirical methods. In the present study, we adapt this formulation to a density functional theory framework and analyze its performance. In particular, it is implemented in an efficient all-electron DFT code with Gaussian basis functions, suitable for quantum-dynamics simulations of large molecular systems. At variance with the case of the tight-binding calculations reported in the literature, we find that now the initial perturbation to drive the system out of equilibrium plays a fundamental role in the stability of the electron dynamics. The equation of motion used in previous tight-binding implementations with massive electrodes has to be modified to produce a stable and unidirectional current during time propagation in time-dependent DFT simulations using much smaller leads. Moreover, we propose a procedure to get rid of the dependence of the current-voltage curves on the rate parameter. This method is employed to obtain the current-voltage characteristic of saturated and unsaturated hydrocarbons of different lengths, with very promising prospects.

The magnetic structure of ß-cobalt hydroxide and the effect of spin-orientation.

Synchrotron X-ray and neutron diffraction experiments at various temperatures, down to 3 K, along with ab initio calculations, are carried out to elucidate the magnetic order of layered ß-cobalt-hydroxide. This combination of techniques allows for the unambiguous assignment of the magnetic structure of this material. Our results confirm that below the Néel temperature high-spin cobalt centers are ferromagnetically coupled within a layer, and antiferromagnetically coupled across layers (magnetic propagation vector k = (0,0,½)), in agreement with the indirect interpretation based on magnetic susceptibility measurements. A paramagnetic/antiferromagnetic transition is observed at around 15 K. Moreover, the thermal expansion behavior along the c-lattice direction, perpendicular to the layers, shows an inflection slightly above this temperature, at around 30 K. The neutron diffraction patterns and the non-collinear DFT+U calculations indicate that the magnetization forms an angle of about 35° with the cobalt planes. In particular, for an isolated ferromagnetic layer, the electronic structure calculations reveal sharp cusps on the potential energy surface when the spins point parallel or perpendicular to the planes, suggesting that the ferromagnetic superexchange mechanism is strongly sensitive to the orientation of the magnetic moment.

Nitrosodisulfide [S2NO]- (perthionitrite) is a true intermediate during the "cross-talk" of nitrosyl and sulfide.

Nitrosodisulfide S2NO- is a controversial intermediate in the reactions of S-nitrosothiols with HS- that produce NO and HNO. QM-MM molecular dynamics simulations combined with TD-DFT analysis contribute to a clear identification of S2NO- in water, acetone and acetonitrile, accounting for the UV-Vis signatures and broadening the mechanistic picture of N/S signaling in biochemistry.

Hydrogen-Bond Heterogeneity Boosts Hydrophobicity of Solid Interfaces.

Experimental and theoretical studies suggest that the hydrophobicity of chemically heterogeneous surfaces may present important nonlinearities as a function of composition. In this article, this issue is systematically explored using molecular simulations. The hydrophobicity is characterized by computing the contact angle of water on flat interfaces and the desorption pressure of water from cylindrical nanopores. The studied interfaces are binary mixtures of hydrophilic and hydrophobic sites, with and without the ability to form hydrogen bonds with water, intercalated at different scales. Water is described with the mW coarse-grained potential, where hydrogen-bonds are modeled in the absence of explicit hydrogen atoms, via a three-body term that favors tetrahedral coordination. We found that the combination of particles exhibiting the same kind of coordination with water gives rise to a linear dependence of contact angle with respect to composition, in agreement with the Cassie model. However, when only the hydrophilic component can form hydrogen bonds, unprecedented deviations from linearity are observed, increasing the contact angle and the vapor pressure above their values in the purely hydrophobic interface. In particular, the maximum enhancement is seen when a 35% of hydrogen bonding molecules is randomly scattered on a hydrophobic background. This effect is very sensitive to the heterogeneity length-scale, being significantly attenuated when the hydrophilic domains reach a size of 2 nm. The observed behavior may be qualitatively rationalized via a simple modification of the Cassie model, by assuming a different microrugosity for hydrogen bonding and non-hydrogen bonding interfaces.

Vapor pressure of water nanodroplets.

Classical thermodynamics is assumed to be valid up to a certain length-scale, below which the discontinuous nature of matter becomes manifest. In particular, this must be the case for the description of the vapor pressure based on the Kelvin equation. However, the legitimacy of this equation in the nanoscopic regime can not be simply established, because the determination of the vapor pressure of very small droplets poses a challenge both for experiments and simulations. In this article we make use of a grand canonical screening approach recently proposed to compute the vapor pressures of finite systems from molecular dynamics simulations. This scheme is applied to water droplets, to show that the applicability of the Kelvin equation extends to unexpectedly small lengths, of only 1 nm, where the inhomogeneities in the density of matter occur within spatial lengths of the same order of magnitude as the size of the object. While in principle this appears to violate the main assumptions underlying thermodynamics, the density profiles reveal, however, that structures of this size are still homogeneous in the nanosecond time-scale. Only when the inhomogeneity in the density persists through the temporal average, as it is the case for clusters of 40 particles or less, do the macroscopic thermodynamics and the molecular descriptions depart from each other.