Source author record

Davide Donadio

Davide Donadio appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

21works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

UV-Visible Absorption Spectra of Solvated Molecules by Quantum Chemical Machine Learning

Predicting UV-visible absorption spectra is essential to understanding photochemical processes and designing energy materials. Quantum chemical methods can deliver accurate calculations of UV-visible absorption spectra, but they are computationally expensive, especially for large systems or when one computes line shapes from thermal averages. Here, we present an approach to predicting UV-visible absorption spectra of solvated aromatic molecules by quantum chemistry (QC) and machine learning (ML). We show that a ML model, trained on the high-level QC calculation of the excitation energy of a set of aromatic molecules, can accurately predict the line shape of the lowest-energy UV-visible absorption band of several related molecules with less than 0.1 eV deviation with respect to reference experimental spectra. Applying linear decomposition analysis on the excitation energies, we unveil that our ML models probe vertical excitations of these aromatic molecules primarily by learning the atomic environment of their phenyl rings, which align with the physical origin of the $π\rightarrowπ^\star$ electronic transition. Our study provides an effective workflow that combines ML with quantum chemical methods to accelerate the calculations of UV-visible absorption spectra for various molecular systems.

preprint2021arXiv

Mode Localization and Suppressed Heat Transport in Amorphous Alloys

Glasses usually represent the lower limit for the thermal conductivity of solids, but a fundamental understanding of lattice heat transport in amorphous materials can provide design rules to beat such a limit. Here we investigate the role of mass disorder in glasses by studying amorphous silicon-germanium alloy (a-Si$_{1-x}$Ge$_x$) over the full range of atomic concentration from $x=0$ to $x=1$, using molecular dynamics and the quasi-harmonic Green-Kubo lattice dynamics formalism. We find that the thermal conductivity of a-Si$_{1-x}$Ge$_x$ as a function of $x$ exhibits a smoother U-shape than in crystalline mass-disordered alloys. The main contribution to the initial drop of thermal conductivity at low Ge concentration stems from the localization of otherwise extended modes that make up the lowest 8\% of the population by frequency. Contributions from intermediate frequency modes are decreased more gradually with increasing Ge to reach a broad minimum thermal conductivity between concentrations of Ge from $x=0.25$ to $0.75$.

preprint2020arXiv

Carbon Dioxide, Bicarbonate and Carbonate Ions in Aqueous Solutions at Deep Earth Conditions

We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids at thermodynamic conditions that correspond to the Earth's deep Crust and Upper Mantle. Our first-principles molecular dynamics simulations provide mechanistic insight into the hydration shell of carbon dioxide, bicarbonate and carbonate ions, and on the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperature of 1000 K and higher our simulations can sample the chemical equilibrium of these acid/base reactions, thus allowing us to estimate the chemical composition of diluted carbon dioxide and (bi)carbonate ions as a function of acidity and thermodynamic conditions. We find that, especially at the highest temperature, the acidity of the solution is essential to determine the stability domain of carbon dioxide, bicarbonate and carbonate ions.

preprint2020arXiv

Transferability of neural network potentials for varying stoichiometry: phonons and thermal conductivity of Mn$_x$Ge$_y$ compounds

Germanium manganese compounds exhibit a variety of stable and metastable phases with different stoichiometry. These materials entail interesting electronic, magnetic and thermal properties both in their bulk form and as heterostructures. Here we develop and validate a transferable machine learning potential, based on the high-dimensional neural network formalism, to enable the study of Mn$_x$Ge$_y$ materials over a wide range of compositions. We show that a neural network potential fitted on a minimal training set reproduces successfully the structural and vibrational properties and the thermal conductivity of systems with different local chemical environments, and it can be used to predict phononic effects in nanoscale heterostructures.

preprint2019arXiv

Anisotropic In-Plane Phonon Transport in Ultrathin Silicon Membranes Guided by Nano-Surface-Resonators

Anisotropic phonon transport along different lattice directions of two-dimensional (2D) materials has been observed, however, the effect decreases with increasing the thickness beyond a few atomic layers. Here we establish a novel mechanism to induce anisotropic phonon transport in quasi-2D materials with isotropic symmetry. The phonon propagation is guided by resonance hybridization with surface nanostructures. We demonstrate that the thermal conductivity of 3 nm-thick silicon membrane with surface nanofins is greater by $\sim50\%$ parallel to the fins than that perpendicular to the fins.

preprint2019arXiv

Influence of Boundaries and Thermostatting on Nonequilibrium Molecular Dynamics Simulations of Heat Conduction in Solids

Nonequilibrium molecular dynamics (NEMD) has been extensively used to study thermal transport at various length scales in many materials. In this method, two local thermostats at different temperatures are used to generate a nonequilibrium steady state with a constant heat flux. Conventionally, the thermal conductivity of a finite system is calculated as the ratio between the heat flux and the temperature gradient extracted from the linear part of the temperature profile away from the local thermostats. Here we show that, with a proper choice of the thermostat, the nonlinear part of the temperature profile should actually not be excluded in thermal transport calculations. We compare NEMD results against those from the atomistic Green's function method in the ballistic regime, and those from the homogeneous nonequilibrium molecular dynamics method in the ballistic-to-diffusive regime. These comparisons suggest that in all the transport regimes, one should directly calculate the thermal conductance from the temperature difference between the heat source and sink and, if needed, convert it to the thermal conductivity by multiplying it with the system length. Furthermore, we find that the Langevin thermostat outperforms the Nosé-Hoover (chain) thermostat in NEMD simulations because of its stochastic and local nature. We show that this is particularly important for studying asymmetric carbon-based nanostructures, for which the Nosé-Hoover thermostat can produce artifacts leading to unphysical thermal rectification. Our findings are important to obtain correct results from molecular dynamics simulations of nanoscale heat transport as the accuracy of the interatomic potentials is rapidly improving.

preprint2016arXiv

Microscopic Mechanism and Kinetics of Ice Formation at Complex Interfaces: Zooming in on Kaolinite

Most ice in nature forms thanks to impurities which boost the exceedingly low nucleation rate of pure supercooled water. However, the microscopic details of ice nucleation on these substances remain largely unknown. Here, we have unraveled the molecular mechanism and the kinetics of ice formation on kaolinite, a clay mineral playing a key role in climate science. We find that the formation of ice at strong supercooling in the presence of this clay is twenty orders of magnitude faster than homogeneous freezing. The critical nucleus is substantially smaller than that found for homogeneous nucleation and, in contrast to the predictions of classical nucleation theory (CNT), it has a strong 2D character. Nonetheless, we show that CNT describes correctly the formation of ice at this complex interface. Kaolinite also promotes the exclusive nucleation of hexagonal ice, as opposed to homogeneous freezing where a mixture of cubic and hexagonal polytypes is observed.

preprint2015arXiv

Force and heat current formulas for many-body potentials in molecular dynamics simulation with applications to thermal conductivity calculations

We derive expressions of interatomic force and heat current for many-body potentials such as the Tersoff, the Brenner, and the Stillinger-Weber potential used extensively in molecular dynamics simulations of covalently bonded materials. Although these potentials have a many-body nature, a pairwise force expression that follows Newton's third law can be found without referring to any partition of the potential. Based on this force formula, a stress applicable for periodic systems can be unambiguously defined. The force formula can then be used to derive the heat current formulas using a natural potential partitioning. Our heat current formulation is found to be equivalent to most of the seemingly different heat current formulas used in the literature, but to deviate from the stress-based formula derived from two-body potential. We validate our formulation numerically on various systems descried by the Tersoff potential, namely three-dimensional silicon and diamond, two-dimensional graphene, and quasi-one-dimensional carbon nanotube. The effects of cell size and time used in the simulation are examined.

preprint2015arXiv

From classical to quantum and back: Hamiltonian coupling of classical and Path Integral models of atoms

In computer simulations, quantum delocalization of atomic nuclei can be modeled making use of the Path Integral (PI) formulation of quantum statistical mechanics. This approach, however, comes with a large computational cost. By restricting the PI modeling to a small region of space, this cost can be significantly reduced. In the present work we derive a Hamiltonian formulation for a bottom-up, theoretically solid simulation protocol that allows molecules to change their resolution from quantum-mechanical to classical and vice versa on the fly, while freely diffusing across the system. This approach renders possible simulations of quantum systems at constant chemical potential. The validity of the proposed scheme is demonstrated by means of simulations of low temperature parahydrogen. Potential future applications include simulations of biomolecules, membranes, and interfaces.

preprint2015arXiv

Mechanical Tuning of Thermal Transport in a Molecular Junction

Understanding and controlling heat transport in molecular junctions would provide new routes to design nanoscale coupled electronic and phononic devices. Using first principles full quantum calculations, we tune thermal conductance of a molecular junction by mechanically compressing and extending a short alkane chain connected to graphene leads. We find that the thermal conductance of the compressed junction drops by half in comparison to the extended junction, making it possible to turn on and off the heat current. The low conductance of the off state does not vary by further approaching the leads and stems from the suppression of the transmission of the in--plane transverse and longitudinal channels. Furthermore, we show that misalignment of the leads does not reduce the conductance ratio. These results also contribute to the general understanding of thermal transport in molecular junctions.

preprint2014arXiv

Length-dependent thermal conductivity in suspended single-layer graphene

Graphene exhibits extraordinary electronic and mechanical properties, and extremely high thermal conductivity. Being a very stable atomically thick membrane that can be suspended between two leads, graphene provides a perfect test platform for studying thermal conductivity in two-dimensional systems, which is of primary importance for phonon transport in low-dimensional materials. Here we report experimental measurements and non-equilibrium molecular dynamics simulations of thermal conduction in suspended single layer graphene as a function of both temperature and sample length. Interestingly and in contrast to bulk materials, when temperature at 300K, thermal conductivity keeps increasing and remains logarithmic divergence with sample length even for sample lengths much larger than the average phonon mean free path. This result is a consequence of the two-dimensional nature of phonons in graphene and provides fundamental understanding into thermal transport in two-dimensional materials.

preprint2014arXiv

Statistical Mechanics of Hamiltonian Adaptive Resolution Simulations

The Adaptive Resolution Scheme (AdResS) is a hybrid scheme that allows one to treat a molecular system with different levels of resolution depending on the location of the molecules. The construction of a Hamiltonian based on the this idea (H-AdResS) allows one to formulate the usual tools of ensembles and statistical mechanics. We present a number of exact and approximate results that provide a statistical mechanics foundation for this simulation method. We also present simulation results that illustrate the theory.

preprint2013arXiv

Divergence of the Thermal Conductivity in Uniaxially Strained Graphene

We investigate the effect of strain and isotopic disorder on thermal transport in suspended graphene by equilibrium molecular dynamics simulations. We show that the thermal conductivity of unstrained graphene, calculated from the fluctuations of the heat current at equilibrium is finite and converges with size at finite temperature. In contrast, the thermal conductivity of strained graphene diverges logarithmically with the size of the models, when strain exceeds a relatively large threshold value of 2%. An analysis of phonon populations and lifetimes explains the divergence of the thermal conductivity as a consequence of changes in the occupation of low-frequency out-of-plane phonons and an increase in their lifetimes due to strain.

preprint2013arXiv

Monte Carlo adaptive resolution simulation of multicomponent molecular liquids

Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The non-matching properties of high- and low-resolution models, however, lead to thermodynamic imbalances between the system's subdomains. Such inhomogeneities can be healed by appropriate compensation forces, whose calculation requires nontrivial iterative procedures. In this work we employ the recently developed Hamiltonian Adaptive Resolution Simulation method to perform Monte Carlo simulations of a binary mixture, and propose an efficient scheme, based on Kirkwood Thermodynamic Integration, to regulate the thermodynamic balance of multi-component systems.

preprint2013arXiv

Thermal conductivity of one-, two- and three-dimensional carbon

Carbon atoms can form structures in one, two, and three dimensions due to its unique chemical versatility. In terms of thermal conductivity, carbon polymorphs cover a wide range from very low values with amorphous carbon to very high values with diamond, carbon nanotubes and graphene. Schwarzites are a class of three-dimensional fully covalent sp$^2$-bonded carbon polymorphs, with the same local chemical environment as graphene and carbon nanotubes, but negative Gaussian curvature. We calculate the thermal conductivity of a (10,0) carbon nanotube, graphene, and two schwarzites with different curvature, by molecular dynamics simulations based on the Tersoff empirical potential. We find that schwarzites present a thermal conductivity two orders of magnitude smaller than nanotubes and graphene. The reason for such large difference is explained by anharmonic lattice dynamics calculations, which show that phonon group velocities and mean free paths are much smaller in schwarzites than in nanotubes and graphene. Their reduced thermal conductivity, in addition to tunable electronic properties, indicate that schwarzites could pave the way towards all-carbon thermoelectric technology with high conversion efficiency.

preprint2012arXiv

Atomistic simulations of heat transport in real-scale silicon nanowire devices

Utilizing atomistic lattice dynamics and scattering theory, we study thermal transport in nanodevices made of 10 nm thick silicon nanowires, from 10 to 100 nm long, sandwiched between two bulk reservoirs. We find that thermal transport in devices differs significantly from that of suspended extended nanowires, due to phonon scattering at the contact interfaces. We show that thermal conductance and the phonon transport regime can be tuned from ballistic to diffusive by varying the surface roughness of the nanowires and their length. In devices containing short crystalline wires phonon tunneling occurs and enhances the conductance beyond that of single contacts.

preprint2012arXiv

Autocatalytic and cooperatively-stabilized dissociation of water on a stepped platinum surface

Water-metal interfaces are ubiquitous and play a key role in many chemical processes, from catalysis to corrosion. Whereas water adlayers on atomically flat transition metal surfaces have been investigated in depth, little is known about the chemistry of water on stepped surfaces, commonly occurring in realistic situations. Using first-principles simulations we study the adsorption of water on a stepped platinum surface. We find that water adsorbs preferentially at the step edge, forming linear clusters or chains, stabilized by the cooperative effect of chemical bonds with the substrate and hydrogen bonds. In contrast with flat Pt, at steps water molecules dissociate forming mixed hydroxyl/water structures, through an autocatalytic mechanism promoted by hydrogen bonding. Nuclear quantum effects contribute to stabilize partially dissociated cluster and chains. Together with the recently demonstrated attitude of water chains adsorbed on stepped Pt surfaces to transfer protons via thermally activated hopping, these findings candidate these systems as viable proton wires.

preprint2012arXiv

Hamiltonian adaptive resolution simulation for molecular liquids

Adaptive resolution schemes allow the simulation of a molecular fluid treating simultaneously different subregions of the system at different levels of resolution. In this work we present a new scheme formulated in terms of a global Hamiltonian. Within this approach equilibrium states corresponding to well defined statistical ensembles can be generated making use of all standard Molecular Dynamics or Monte Carlo methods. Models at different resolutions can thus be coupled, and thermodynamic equilibrium can be modulated keeping each region at desired pressure or density without disrupting the Hamiltonian framework.

preprint2011arXiv

Atomistic calculation of the thermal conductance of large scale bulk-nanowire junctions

We have developed an efficient scalable kernel method for thermal transport in open systems, with which we have computed the thermal conductance of a junction between bulk silicon and silicon nanowires with diameter up to 10 nm. We have devised scaling laws for transmission and reflection spectra, which allow us to predict the thermal resistance of bulk-nanowire interfaces with larger cross sections than those achievable with atomistic simulations. Our results indicate the characteristic size beyond which atomistic systems can be treated accurately by mesoscopic theories.

preprint2010arXiv

Ab initio investigation of the melting line of nitrogen at high pressure

Understanding the behavior of molecular systems under pressure is a fundamental problem in condensed matter physics. In the case of nitrogen, the determination of the phase diagram and in particular of the melting line, are largely open problems. Two independent experiments have reported the presence of a maximum in the nitrogen melting curve, below 90 GPa, however the position and the interpretation of the origin of such maximum differ. By means of ab initio molecular dynamics simulations based on density functional theory and thermodynamic integration techniques, we have determined the phase diagram of nitrogen in the range between 20 and 100 GPa. We find a maximum in the melting line, related to a transformation in the liquid, from molecular N_2 to polymeric nitrogen accompanied by an insulator-to-metal transition.

preprint2009arXiv

Theoretical investigation of methane under pressure

We present computer simulations of liquid and solid phases of condensed methane at pressures below 25 GPa, between 150 and 300 K, where no appreciable molecular dissociation occurs. We used molecular dynamics (MD) and metadynamics techniques, and empirical potentials in the rigid molecule approximation, whose validity was confirmed a posteriori by carrying out it ab initio MD simulations for selected pressure and temperature conditions. Our results for the melting line are in satisfactory agreement with existing measurements. We find that the fcc crystal transforms into a hcp structure with 4 molecules per unit cell (B phase) at about 10 GPa and 150 K, and that the B phase transforms into a monoclinic high pressure phase above 20 GPa. Our results for solid/solid phase transitions are consistent with those of Raman studies but the phase boundaries estimated in our calculations are at higher pressure than those inferred from spectroscopic data.