RESUMO
Crystals with complicated geometry are often observed with mixed chemical occupancy among Wyckoff sites, presenting a unique challenge for accurate atomic modeling. Similar systems possessing exact occupancy on all the sites can exhibit superstructural ordering, dramatically inflating the unit cell size. In this work, a crystal graph convolutional neural network (CGCNN) is used to predict optimal atomic decorations on fixed crystalline geometries. This is achieved with a site permutation search (SPS) optimization algorithm based on Monte Carlo moves combined with simulated annealing and basin-hopping techniques. Our approach relies on the evidence that, for a given chemical composition, a CGCNN estimates the correct energetic ordering of different atomic decorations, as predicted by electronic structure calculations. This provides a suitable energy landscape that can be optimized according to site occupation, allowing the prediction of chemical decoration in crystals exhibiting mixed or disordered occupancy, or superstructural ordering. Verification of the procedure is carried out on several known compounds, including the superstructurally ordered clathrate compound Rb8Ga27Sb16 and vacancy-ordered perovskite Cs2SnI6, neither of which was previously seen during the neural network training. In addition, the critical temperature of an order-disorder phase transition in solid solution CuZn is probed with our SPS routines by sampling site configuration trajectories in the canonical ensemble. This strategy provides an accurate method for determining favorable decoration in complex crystals and analyzing site occupation at unprecedented speed and scale.
RESUMO
Three polyanionic tellurides, ABa6Cu31Te22 (A = K, Rb, Cs), were synthesized in salt flux. The isostructural tellurides crystallize in a new structure type, in the cubic Pa3 space group with a Wyckoff sequence of d10c2b1 and large unit cell volumes of over 5500 Å3. The structures feature a framework of [CuTe4] tetrahedra and [CuTe3] trigonal pyramids with disorder in the Cu sites. The polyanionic frameworks have large square antiprism and cuboctahedral voids where Ba and alkali metal cations are situated, forming [BaTe8] and [ATe12], respectively. The overall compositions are close to being charge balanced. The large [ATe12] cuboctahedra allowed for significant anisotropic displacement of the A cations, as observed from both single crystal X-ray diffraction and heat capacity studies. Alkali cations rattling together with Cu atom displacement and disorder leads to the dispersion of phonons, thus softening the lattice and subsequently reducing the thermal conductivity. Evaluations of the electronic band structure revealed the occurrence of a narrow bandgap together with the presence of a flat band near the valence band maximum, giving rise to the high thermopower. The Cs and Rb analogues show a slope change in the temperature dependence of electrical resistivity around room temperature, which is typical for semimetals or degenerate semiconductors. For the as-synthesized and unoptimized materials, high values of the thermoelectric figure-of-merit of â¼0.2 were observed at 623 K.
RESUMO
The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.
RESUMO
In the interstellar medium, diols and other prebiotic molecules adsorb onto icy mantles surrounding dust grains. Water in the ice may affect the reactivity and photoionization of these diols. Ethylene glycol (EG), 1,2-propylene glycol, and 1,3-propylene glycol clusters with water clusters were used as a proxy to study these interactions. The diol-water clusters were generated in a continuous supersonic molecular beam, photoionized by synchrotron-based vacuum ultraviolet light from the Advanced Light Source, and subsequently detected by reflectron time-of-flight mass spectrometry. The appearance energies for the detected clusters were determined from the mass spectra, collected at increasing photon energy. Clusters of both diol fragments and unfragmented diols with water were detected. The lowest energy geometry optimized conformers for the observed EG-water clusters and EG fragment-water clusters have been visualized using density functional theory (DFT), providing insight into hydrogen bonding networks and how these affect fragmentation and appearance energy. As the number of water molecules clustered around EG fragments (m/z 31 and 32) increased, the appearance energy for the cluster decreased, indicating a stabilization by water. This trend was supported by DFT calculations. Fragment clusters from 1,2-propylene glycol exhibited a similar trend, but with a smaller energy decrease, and no trend was observed from 1,3-propylene glycol. We discuss and suggest that the reactivity and photoionization of diols in the presence of water depend on the size of the diol, the location of the hydroxyl group, and the number of waters clustered around the diol.
RESUMO
The zinc-antimony phase space has been heavily investigated due to the structural complexity and abundance of high-performing thermoelectric materials. Consequentially, the desire to use zinc and antimony as framework elements to encage rattling cations and achieve phonon-glass-electron-crystal-type properties has remained an enticing goal with only two alkali metal clathrates to date, Cs8Zn18Sb28 and K58Zn122Sb207. Guided by Zintl electron-counting predictions, we explored the Ba-Zn-Pn (Pn = As, Sb) phase space proximal to the expected composition of the type-I clathrate. In situ powder X-ray diffraction studies revealed two "hidden" compounds which can only be synthesized in a narrow temperature range. The ex situ synthesis and crystal growth unveiled that instead of type-I clathrates, compositionally close but structurally different new clathrate-like compounds formed, Ba2Zn5As6 and Ba2Zn5Sb6. These materials crystallize in a unique structure, in the orthorhombic space group Pmna with the Wyckoff sequence i2h6gfe. Single-phase synthesis enabled the exploration of their transport properties. Rattling of the Ba cations in oversized cages manifested low thermal conductivity, which, coupled with the high Seebeck coefficients observed, are prerequisites for a promising thermoelectric material. Potential for further optimization of the thermoelectric performance by aliovalent doping was computationally analyzed.
RESUMO
We characterise the structural properties of the quasi-liquid layer (QLL) at two low-index ice surfaces in the presence of sodium chloride (Na+/Cl-) ions by molecular dynamics simulations. We find that the presence of a high surface density of Na+/Cl- pairs changes the surface melting behaviour from step-wise to gradual melting. The ions lead to an overall increase of the thickness and the disorder of the QLL, and to a low-temperature roughening transition of the air-ice interface. The local molecular structure of the QLL is similar to that of liquid water, and the differences between the basal and primary prismatic surface are attenuated by the presence of Na+/Cl- pairs. These changes modify the crystal growth rates of different facets and the solvation environment at the surface of sea-water ice with a potential impact on light scattering and environmental chemical reactions.
RESUMO
Predicting UV-visible absorption spectra is essential to understand photochemical processes and design 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 predict 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 π âπ* 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.
Assuntos
Aprendizado de Máquina , Espectrofotometria UltravioletaRESUMO
The first arsenic-based clathrate exhibiting superstructural ordering due to optimization of Au-As, As-As, and Ba-Au bonding is reported. Ba8Au16As30 crystallizes in a unique P21/c monoclinic clathrate structure. The synthesis, crystal and electronic structure, and transport properties are discussed.
RESUMO
Layering two-dimensional van der Waals materials provides a high degree of control over atomic placement, which could enable tailoring of vibrational spectra and heat flow at the sub-nanometer scale. Here, using spatially resolved ultrafast thermoreflectance and spectroscopy, we uncover the design rules governing cross-plane heat transport in superlattices assembled from monolayers of graphene (G) and MoS2 (M). Using a combinatorial experimental approach, we probe nine different stacking sequences, G, GG, MG, GGG, GMG, GGMG, GMGG, GMMG, and GMGMG, and identify the effects of vibrational mismatch, interlayer adhesion, and junction asymmetry on thermal transport. Pure G sequences display evidence of quasi-ballistic transport, whereas adding even a single M layer strongly disrupts heat conduction. The experimental data are described well by molecular dynamics simulations, which include thermal expansion, accounting for the effect of finite temperature on the interlayer spacing. The simulations show that an increase of â¼2.4% in the layer separation of GMGMG, relative to its value at 300 K, can lead to a doubling of the thermal resistance. Using these design rules, we experimentally demonstrate a five-layer GMGMG superlattice "thermal metamaterial" with an ultralow effective cross-plane thermal conductivity comparable to that of air.
RESUMO
Water is one of the most important, yet least understood, liquids in nature. Many anomalous properties of liquid water originate from its well-connected hydrogen bond network1, including unusually efficient vibrational energy redistribution and relaxation2. An accurate description of the ultrafast vibrational motion of water molecules is essential for understanding the nature of hydrogen bonds and many solution-phase chemical reactions. Most existing knowledge of vibrational relaxation in water is built upon ultrafast spectroscopy experiments2-7. However, these experiments cannot directly resolve the motion of the atomic positions and require difficult translation of spectral dynamics into hydrogen bond dynamics. Here, we measure the ultrafast structural response to the excitation of the OH stretching vibration in liquid water with femtosecond temporal and atomic spatial resolution using liquid ultrafast electron scattering. We observed a transient hydrogen bond contraction of roughly 0.04 Å on a timescale of 80 femtoseconds, followed by a thermalization on a timescale of approximately 1 picosecond. Molecular dynamics simulations reveal the need to treat the distribution of the shared proton in the hydrogen bond quantum mechanically to capture the structural dynamics on femtosecond timescales. Our experiment and simulations unveil the intermolecular character of the water vibration preceding the relaxation of the OH stretch.
RESUMO
Solid solutions of Yb2-xAxCdSb2 (A = Ca, Sr, Eu; x ≤ 1) are of interest for their promising thermoelectric (TE) properties. Of these solid solutions, Yb2-xCaxCdSb2 has end members with different crystal structures. Yb2CdSb2 crystallizes in the polar space group Cmc21, whereas Ca2CdSb2 crystallizes in the centrosymmetric space group Pnma. Other solid solutions, Yb2-xAxCdSb2 (A = Sr, Eu), crystallize in the polar space group for x ≤ 1, and compositions with x ≥ 1 have not been reported. Both structure types are composed of corner-sharing CdSb4 tetrahedra condensed into sheets that differ by the stacking of the layers. Single crystals of the solid solution Yb2-xCaxCdSb2 (x = 0-1) were studied to elucidate the structural transition between the Yb2CdSb2 and Ca2CdSb2 structure types. For x ≤ 1, the structures remain in the polar space group Cmc21. As the Ca content is increased, a positional disorder arises in the intralayer cation sites (Yb2/Ca2) and the Cd site, resulting in inversion of the CdSb4 tetrahedral chain. This phenomenon could be indicative of an intergrowth of the opposing space group. The TE properties of polycrystalline samples of Yb2-xCaxCdSb2 (x ≤ 1) were measured from 300 to 525 K. The lattice thermal conductivity is extremely low (0.3-0.4 W/m·K) and the Seebeck coefficients are high (100-180 µV/K) across the temperature range. First-principles calculations show a minimum in the thermal conductivity for the x = 0.3 composition, in good agreement with experimental data. The low thermal conductivity stems from the acoustic branches being confined to low frequencies and a large number of phonon scattering channels provided by the localized optical branches. The TE quality factor of the Yb1.7A0.3CdSb2 (A = Ca, Sr, Eu) series has been calculated and predicts that the A = Ca and Sr solid solutions may not improve with carrier concentration optimization but that the Eu series is worthy of additional modifications. Overall, the x = 0.3 compositions provide the highest zT because they provide the best electronic properties with the lowest thermal conductivity.
RESUMO
Anomalous heat transport in one-dimensional nanostructures, such as nanotubes and nanowires, is a widely debated problem in condensed matter and statistical physics, with contradicting pieces of evidence from experiments and simulations. Using a comprehensive modeling approach, comprised of lattice dynamics and molecular dynamics simulations, we proved that the infinite length limit of the thermal conductivity of a (10,0) single-wall carbon nanotube is finite but this limit is reached only for macroscopic lengths due to a thermal phonon mean free path of several millimeters. Our calculations showed that the extremely high thermal conductivity of this system at room temperature is dictated by quantum effects. Modal analysis showed that the divergent nature of thermal conductivity, observed in one-dimensional model systems, is suppressed in carbon nanotubes by anharmonic scattering channels provided by the flexural and optical modes with polarization in the plane orthogonal to the transport direction.
RESUMO
Some organic pollutants in snowpacks undergo faster photodegradation than in solution. One possible explanation for such effect is that their UV-visible absorption spectra are shifted toward lower energy when the molecules are adsorbed at the air-ice interface. However, such bathochromic shift is difficult to measure experimentally. Here, we employ a multiscale/multimodel approach that combines classical and first-principles molecular dynamics, quantum chemical methods, and statistical learning to compute the light absorption spectra of two phenolic molecules in different solvation environments at the relevant thermodynamic conditions. Our calculations provide an accurate estimate of the bathochromic shift of the lowest-energy UV-visible absorption band when these molecules are adsorbed at the air-ice interface, and they shed light into its molecular origin.
RESUMO
Intercalation offers a promising way to alter the physical properties of two-dimensional (2D) layered materials. Here, we investigate the electronic and vibrational properties of 2D layered MoSe2 intercalated with atomic manganese at ambient and high pressure up to 7 GPa by Raman scattering and electronic structure calculations. The behavior of optical phonons is studied experimentally with a diamond anvil cell and computationally through density functional theory calculations. Experiment and theory show excellent agreement in optical phonon behavior. The previously Raman inactive A2u mode is activated and enhanced with intercalation and pressure, and a new Raman mode appears upon decompression, indicating a possible onset of a localized structural transition, involving the bonding or trapping of the intercalant in 2D layered materials. Density functional theory calculations reveal a shift of the Fermi level into the conduction band and spin polarization in MnxMoSe2 that increases at low Mn concentrations and low pressure. Our results suggest that intercalation and pressurization of van der Waals materials may allow one to obtain dilute magnetic semiconductors with controllable properties, providing a viable route for the development of new materials for spintronic applications.
RESUMO
Snowpacks contain a wide variety of inorganic and organic compounds, including some that absorb sunlight and undergo direct photoreactions. How the rates of these reactions in, and on, ice compare to rates in water is unclear: some studies report similar rates, while others find faster rates in/on ice. Further complicating our understanding, there is conflicting evidence whether chemicals react more quickly at the air-ice interface compared to in liquid-like regions (LLRs) within the ice. To address these questions, we measured the photodegradation rate of guaiacol (2-methoxyphenol) in various sample types, including in solution, in ice, and at the air-ice interface of nature-identical snow. Compared to aqueous solution, we find modest rate constant enhancements (increases of 3- to 6-fold) in ice LLRs, and much larger enhancements (of 17- to 77-fold) at the air-ice interface of nature-identical snow. Our computational modeling suggests the absorption spectrum for guaiacol red-shifts and increases on ice surfaces, leading to more light absorption, but these changes explain only a small portion (roughly 2 to 9%) of the observed rate constant enhancements in/on ice. This indicates that increases in the quantum yield are primarily responsible for the increased photoreactivity of guaiacol on ice; relative to solution, our results suggest that the quantum yield is larger by a factor of roughly 3-6 in liquid-like regions and 12-40 at the air-ice interface.
Assuntos
Guaiacol , Gelo , Fotólise , Luz Solar , ÁguaRESUMO
We investigate the effect of pressure, temperature and acidity on the composition of water-rich carbon-bearing fluids under 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 into the pathways of the acid/base reactions that convert these carbon species into one another in aqueous solutions. At temperatures 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 CO2vs. HCO3-vs. CO32-.
RESUMO
Three novel unconventional clathrates with unprecedented III-V semiconducting frameworks have been synthesized: Cs8In27Sb19, Cs8Ga27Sb19, and Rb8Ga27Sb19. These clathrates represent the first examples of tetrel-free clathrates that are completely composed of main group elements. All title compounds crystallize in an ordered superstructure of clathrate-I in the Ia3Ì space group (No. 206; Z = 8). In the clathrate framework, a full ordering of {Ga or In} and Sb is observed by a combination of high-resolution synchrotron single-crystal and powder X-ray diffraction techniques. Density functional theory (DFT) calculations show that all three clathrates are energetically stable with relaxed lattice constants matching the experimental data. Due to the complexity of the crystal structure composed of heavy elements, the reported clathrates exhibit ultralow thermal conductivities of less than 1 W·m-1·K-1 at room temperature. All compounds are predicted and experimentally confirmed to be narrow-bandgap p-type semiconductors with high Seebeck thermopower values, up to 250 µV·K-1 at 300 K for Cs8In27Sb19. The latter compound shows carrier concentrations and mobilities, 1.42 × 1015 cm-3 and 880 cm2 ·V-1·s-1, which are on par with the values for parent binary InSb, one of the best electronic semiconductors. The high hole carrier mobility is uncommon for complex bulk materials and a highly desirable trait, opening ways to design semiconducting materials based on tunable III-V clathrates.
RESUMO
Understanding the thermal transport mechanisms in amorphous organic materials is of great importance to solve hot-spot issues in organic-electronics nanodevices. Here we studied thermal transport in two popular molecular electronic materials, N,N-dicarbazolyl-3,5-benzene (mCP) and N,N'-diphenyl-N,N'-di(3-methylphenyl)-(1,1'-biphenyl)-4,4'diamine (TPD), in the amorphous state by molecular dynamics simulations. We found that due to the softness of organic materials, the low thermal conductivity of both systems can be greatly enhanced under pressure. Notably, in such systems, the convective term of heat flux provides an important contribution to thermal transport as it cross-correlates with the Virial term in the Green-Kubo formula. Mode diffusivity calculations reveal that low-frequency modes can contribute significantly to thermal transport in both mCP and TPD. By increasing the pressure, the sound velocity and relaxation time of such low-frequency modes can be enhanced, and a part of these modes converts from diffusons to propagons. The cooperation of these three effects is responsible for the strong pressure dependence of thermal transport in amorphous organic systems. Molecular pair heat flux calculations demonstrate that heat transfer mainly happens between pairs of molecules with distances below 1.4 nm. This work paves the way for the optimization of thermal transport in amorphous organic materials widely used in opto-electronics, e.g. as OLED and OPV.
RESUMO
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 into 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.