Nature of partial sigma bond

This study investigates the formation of partial sigma (σ) covalent bonds in experimentally synthesizable biradicals formed from hydrogenated and fluorinated C8, C20, and C60 cage structures, by assessing their stability, geometry, and bonding character in singlet and triplet states using restricted B3LYP-D3/6-31+G(d,p) theory, natural bond orbital (NBO) analysis, and complete active space self-consistent field (CASSCF) method. The results show that these partial σCC bonds have Wiberg bond orders of 0.38 to 0.48 and bond lengths ranging from 2.62 Å to 5.93 Å. Cage size influences the characteristics of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), with electrons favoring more antibonding orbitals in smaller cages where electrons reside more on the exterior of the cage and favoring bonding orbitals in larger ones where electrons are more in the interior. Fluorination enhances electron density on bonding orbitals. The analysis further clarified that the differentiation between antibonding and bonding features of HOMOs and LUMOs extends beyond merely electron transfer from s- to p-atomic orbitals, also noting possible interactions of the same symmetry repel. The study also introduces hyperconjugation from α-position CH bonds as a factor in stabilizing partial σ-bond formation. The results also caution against the use of broken symmetry methodology in unrestricted SCF wavefunctions for biradicals, such as those in this study as it may cause large spin contamination and thus errors in the calculated electronic properties results.

High-throughput dataset of impurity adsorption on common catalysts in biomass upgrading applications

An extensive dataset consisting of adsorption energies of pernicious impurities present in biomass upgrading processes on common catalysts and support materials has been generated. This work aims to inform catalyst and process development for the conversion of biomass-derived feedstocks to fuels and chemicals. A high-throughput workflow was developed to execute density functional theory calculations for a diverse set of atomic (Al, B, Ca, Cl, Fe, K, Mg, Mn, N, Na, P, S, Si, Zn) and molecular (COS, H2S, HCl, HCN, K2O, KCl, NH3) species on 35 unique surfaces for transition-metal (Ag, Au, Co, Cu, Fe, Ir, Ni, Pd, Pt, Re, Rh, Ru) and metal-oxide (Al2O3, MgO, anatase-TiO2, rutile-TiO2, ZnO, ZrO2) catalysts and supports. Approximately 3,000 unique adsorption geometries and corresponding adsorption energies were obtained.

Ultrahigh Stability and Operation Performance in Bi-doped GeTe/Sb2Te3 Superlattices Achieved by Tailoring Bonding and Structural Properties

Changes in bond types and the reversible switching process between metavalent and covalent bonds are related to the operating mechanism of the phase-change (PC) behavior. Thus, controlling the bonding characteristics is the key to improving the PC memory performance. In this study, we have controlled the bonding characteristics of GeTe/Sb2Te3 superlattices (SLs) via bismuth (Bi) doping. The incorporation of Bi into the GeTe sublayers tailors the metavalent bond. We observed significant improvement in device reliability, set speed, and power consumption induced upon increasing Bi incorporation. The introduction of Bi was found to suppress the change in density between the SET and RESET states, resulting in a significant increase in device reliability. The reduction in Peierls distortion, leading to a more octahedral-like atomic arrangement, intensifies electron-phonon coupling with increased bond polarizability, which are responsible for the fast set speed and low power consumption. This study demonstrates how the structural and thermodynamic changes in phase change materials alter phase change characteristics due to systematic changes of bonding and provides an important methodology for the development of PC devices.

Lithium Orbital Hybridization Chemistry to Stimulate Oxygen Redox with Reversible Phase Evolution in Sodium-Layered Oxide Cathodes

Searching for high energy-density electrode materials for sodium ion batteries has revealed Na-deficient intercalation compounds with lattice oxygen redox as promising high-capacity cathodes. However, anionic redox reactions commonly encountered poor electrochemical reversibility and unfavorable structural transformations during dynamic (de)sodiation processes. To address this issue, we employed lithium orbital hybridization chemistry to create Na-O-Li configuration in a prototype P2-layered Na43/60Li1/20Mg7/60Cu1/6Mn2/3O2 (P2-NaLMCM') cathode material. That Li+ ions, having low electronegativity, reside in the transition metal slabs serves to stimulate unhybridized O 2p orbitals to facilitate the stable capacity contribution of oxygen redox at high state of charge. The prismatic-type structure evolving to an intergrowth structure of the Z phase at high charging state could be simultaneously alleviated by reducing the electrostatic repulsion of O-O layers. As a consequence, P2-NaLMCM' delivers a high specific capacity of 183.8 mAh g-1 at 0.05 C and good cycling stability with a capacity retention of 80.2% over 200 cycles within the voltage range of 2.0-4.5 V. Our findings provide new insights into both tailoring oxygen redox chemistry and stabilizing dynamic structural evolution for high-energy battery cathode materials.

An automated protocol to construct flexibility parameters for classical forcefields: applications to metal-organic frameworks

In this work, forcefield flexibility parameters were constructed and validated for more than 100 metal-organic frameworks (MOFs). We used atom typing to identify bond types, angle types, and dihedral types associated with bond stretches, angle bends, dihedral torsions, and other flexibility interactions. Our work used Manz's angle-bending and dihedral-torsion model potentials. For a crystal structure containing N atoms in its unit cell, the number of independent flexibility interactions is 3(N atoms - 1). Because the number of bonds, angles, and dihedrals is normally much larger than 3(N atoms - 1), these internal coordinates are redundant. To reduce (but not eliminate) this redundancy, our protocol prunes dihedral types in a way that preserves symmetry equivalency. Next, each dihedral type is classified as non-rotatable, hindered, rotatable, or linear. We introduce a smart selection method that identifies which particular torsion modes are important for each rotatable dihedral type. Then, we computed the force constants for all flexibility interactions together via LASSO regression (i.e., regularized linear least-squares fitting) of the training dataset. LASSO automatically identifies and removes unimportant forcefield interactions. For each MOF, the reference dataset was quantum-mechanically-computed in VASP via DFT with dispersion and included: (i) finite-displacement calculations along every independent atom translation mode, (ii) geometries randomly sampled via ab initio molecular dynamics (AIMD), (iii) the optimized ground-state geometry using experimental lattice parameters, and (iv) rigid torsion scans for each rotatable dihedral type. After training, the flexibility model was validated across geometries that were not part of the training dataset. For each MOF, we computed the goodness of fit (R-squared value) and the root-mean-squared error (RMSE) separately for the training and validation datasets. We compared flexibility models with and without bond-bond cross terms. Even without cross terms, the model yielded R-squared values of 0.910 (avg across all MOFs) ± 0.018 (st. dev.) for atom-in-material forces in the validation datasets. Our SAVESTEPS protocol should find widespread applications to parameterize flexible forcefields for material datasets. We performed molecular dynamics simulations using these flexibility parameters to compute heat capacities and thermal expansion coefficients for two MOFs.

Engineered Nanoporous Frameworks for Adsorption Cooling Applications

The energy demand for traditional vapor-compressed technology for space cooling continues to soar year after year due to global warming and the increasing human population's need to improve living and working conditions. Thus, there is a growing demand for eco-friendly technologies that use sustainable or waste energy resources. This review discusses the properties of various refrigerants used for adsorption cooling applications followed by a brief discussion on the thermodynamic cycle. Next, sorbents traditionally used for cooling are reviewed to emphasize the need for advanced capture materials with superior properties to improve refrigerant sorption. The remainder of the review focus on studies using engineered nanoporous frameworks (ENFs) with various refrigerants for adsorption cooling applications. The effects of the various factors that play a role in ENF-refrigerant pair selection, including pore structure/dimension/shape, morphology, open-metal sites, pore chemistry and possible presence of defects, are reviewed. Next, in-depth insights into the sorbent-refrigerant interaction, and pore filling mechanism gained through a combination of characterization techniques and computational modeling are discussed. Finally, we outline the challenges and opportunities related to using ENFs for adsorption cooling applications and provide our views on the future of this technology.

PACMAN: A Robust Partial Atomic Charge Predicter for Nanoporous Materials Based on Crystal Graph Convolution Networks

We report a fast and easy method (PACMAN) to assign partial atomic charges on metal-organic framework (MOF) and covalent-organic framework (COF) crystal structures based on graph convolution networks (GCNs) trained on >1.8 million high-fidelity partial atomic charge data obtained from the Quantum Metal-Organic Framework (QMOF) database. The developed model shows outstanding performance, achieving a mean absolute error (MAE) of 0.0055 e (test set performance) while maintaining consistency with DDEC6, Bader, and CM5 charges across diverse chemistry and topologies of MOFs and COFs. We find that the new method accurately assigns partial atomic charges for ion-containing nanoporous materials, which has not been possible in previous machine learning (ML) models. Grand canonical Monte Carlo (GCMC) simulation results for CO2 and N2 uptakes and the Widom particle insertion calculation for Henry's law constant of water results based on PACMAN and the original DDEC6 charges show excellent agreements compared to other ML models reported in the literature. The runtime analysis of the new method demonstrates that the partial atomic charges of MOF and COF structures with up to 500 atoms can be obtained in less than 10 s. An easy-to-use web interface has been developed to facilitate the adoption of the developed model.

Pt-based graphene quantum dots for water dissociation

It is widely recognized that Pt nanostructures exhibit favorable catalytic properties for several important technological reactions. Furthermore, selecting an appropriate support has the potential to enhance the catalytic activity of these materials. In this study, we investigate Pt nanoparticles deposited on quantum dots using quantum chemical calculations. We explore the utilization of low-dimensional carbonaceous support by employing graphene quantum dots (GQDs), which offer abundant active sites, such as edges, and diverse conformations. This provides excellent tuning possibilities for both chemical and physical properties. Our goal is to gather information on the alterations in electronic properties, charge redistribution and reactivity of platinum particles on GQD, also analyzing their potential role as catalysts in the water dissociation reaction. Based on thermodynamic and kinetic considerations, our calculations suggest that a Pt3nanoparticle adsorbed on the edge of the GQD exhibits favorable energetics, leading to a promising catalytic material.

Radiation-Induced Defects in Uranyl Trinitrate Solids

Actinides are inherently radioactive; thus, ionizing radiation is emitted by these elements can have profound effects on its surrounding chemical environment through the formation of free radical species. While previous work has noted that the presence of free radicals in the system impacts the redox state of the actinides, there is little atomistic understanding of how these metal cations interact with free radicals. Herein, we explore the effects of radiation (UV and γ) on three U(VI) trinitrate complexes, M[UO2(NO3)3] (where M=K+, Rb+, Cs+), and their respective nitrate salts in the solid state via electron paramagnetic resonance (EPR) and Raman spectroscopy paired with Density Functional Theory (DFT) methods. We find that the alkali salts form nitrate radicals under UV and γ irradiation, but also note the presence of additional degradation products. M[UO2(NO3)3] solids also form nitrate radicals and additional DFT calculations indicate the species corresponds to a change from the bidentate bound nitrate anion into a monodentate NO3 • radical. Computational studies also highlight the need to include the second sphere coordination environment around the [UO2(NO3)3]0,1 species to gain agreement between the experimental and predicted EPR signatures.

Quantifying the Uncertainty of Force Field Selection on Adsorption Predictions in MOFs

Comparisons between simulated and experimental adsorption isotherms in MOFs are fraught with challenges. On the experimental side, there is significant variation between isotherms measured on the same system, with a significant percentage (∼20%) of published data being considered outliers. On the simulation side, force fields are often chosen "off-the-shelf" with little or no validation. The effect of this choice on the reliability of simulated adsorption predictions has not yet been rigorously quantified. In this work, we fill this gap by systematically quantifying the uncertainty arising from force field selection on adsorption isotherm predictions. We choose methane adsorption, where electrostatic interactions are negligible, to independently study the effect of the framework Lennard-Jones parameters on a series of prototypical materials that represent the most widely studied MOF "families". Using this information, we compute an adsorption "consensus isotherm" from simulations, including a quantification of uncertainty, and compare it against a manually curated set of experimental data from the literature. By considering many experimental isotherms measured by different groups and eliminating outliers in the data using statistical analysis, we conduct a rigorous comparison that avoids the pitfalls of the standard approach of comparing simulation predictions to a single experimental data set. Our results show that (1) the uncertainty in simulated isotherms can be as large as 15% and (2) standard force fields can provide reliable predictions for some systems but can fail dramatically for others, highlighting systematic shortcomings in those models. Based on this, we offer recommendations for future simulation studies of adsorption, including high-throughput computational screening of MOFs.

Quantitative Analogue Simulation of Planar Molecules

Synthetic quantum systems provide a pathway for exploring the physics of complex quantum matter in a programmable fashion. This approach becomes particularly advantageous when it comes to systems that are thermodynamically unfavorable. By sculpting the potential landscape of Cu(111) surfaces with carbon monoxide quantum corrals in a cryogenic scanning tunneling microscope, we created analogue simulators of planar organic molecules, including antiaromatic and non-Kekulé species that are generally reactive or unstable. Spectroscopic imaging of such synthetic molecules reveals close replications of molecular orbitals obtained from ab initio calculations of the organic molecules. We further illustrate the quantitative nature of such analogue simulators by faithful extraction of bond orders and global aromaticity indices, which are otherwise technically daunting using real molecules. Our approach therefore sets the stage for new research frontiers pertaining to the quantum physics and chemistry of designer nanostructures.

Supramolecular compounds assembled from the heteroleptic tetrahedral complex [{CpMo(CO)2}2(μ,η2-AsSb)] and metal salts

The reaction of the tetrahedral complex [{CpMo(CO)2}2(μ,η2-AsSb)] with CuI and AgI salts is presented which gives unprecedented neutral and cationic supramolecular aggregates featuring mixed As/Sb-donor molecules as ligands/linkers between metal ions.

pyCOFBuilder: A Python Package for Automated Creation of Covalent Organic Framework Models Based on the Reticular Approach

Covalent organic frameworks (COFs) have gained significant popularity in recent years due to their unique ability to provide a high surface area and customizable pore geometry and chemistry, making them an ideal choice for a wide range of applications. However, exploring COFs experimentally can be arduous and time-consuming due to their immense number of potential structures. As a result, computational high-throughput studies have become an attractive option. Nevertheless, generating COF structures can also be a challenging and time-consuming task. To address this challenge, here, we introduce the pyCOFBuilder, an open-source Python package designed to facilitate the generation of COF structures for computational studies. The pyCOFBuilder software provides an easy-to-use set of functionalities to generate COF structures following the reticular approach. In this paper, we describe the implementation, main features, and capabilities of the pyCOFBuilder, demonstrating its utility for generating COF structures with varying topologies and chemical properties. pyCOFBuilder is freely available on GitHub at https://github.com/lipelopesoliveira/pyCOFBuilder.

Ultralow Glassy Thermal Conductivity and Controllable, Promising Thermoelectric Properties in Crystalline o-CsCu5S3

We thoroughly investigated the anharmonic lattice dynamics and microscopic mechanisms of the thermal and electronic transport characteristics in orthorhombic o-CsCu5S3 at the atomic level. Taking into account the phonon energy shifts and the wave-like tunneling phonon channel, we predict an ultralow κL of 0.42 w/mK at 300 K with an extremely weak temperature dependence following ∼T-0.33. These findings agree well with experimental values along with the parallel to the Bridgman growth direction. The κL in o-CsCu5S3 is suppressed down to the amorphous limit, primarily due to the unconventional Cu-S bonding induced by the p-d hybridization antibonding state coupled with the stochastic oscillation of Cs atoms. The nonstandard temperature dependence of κL can be traced back to the critical or dominant role of wave-like tunneling of phonon contributions in thermal transport. Moreover, the p-d hybridization of Cu(3)-S bonding results in the formation of a valence band with "pudding-mold" and high-degeneracy valleys, ensuring highly efficient electron transport characteristics. By properly adjusting the carrier concentration, excellent thermoelectric performance is achieved with a maximum thermoelectric conversion efficiency of 18.4% observed at 800 K in p-type o-CsCu5S3. Our work not only elucidates the anomalous electronic and thermal transport behavior in the copper-based chalcogenide o-CsCu5S3 but also provides insights for manipulating its thermal and electronic properties for potential thermoelectric applications.

Chemically tuned intermediate band states in atomically thin CuxGeSe/SnS quantum material for photovoltaic applications

A new generation of quantum material derived from intercalating zerovalent atoms such as Cu into the intrinsic van der Waals gap at the interface of atomically thin two-dimensional GeSe/SnS heterostructure is designed, and their optoelectronic features are explored for next-generation photovoltaic applications. Advanced ab initio modeling reveals that many-body effects induce intermediate band (IB) states, with subband gaps (~0.78 and 1.26 electron volts) ideal for next-generation solar devices, which promise efficiency greater than the Shockley-Queisser limit of ~32%. The charge carriers across the heterojunction are both energetically and spontaneously spatially confined, reducing nonradiative recombination and boosting quantum efficiency. Using this IB material in a solar cell prototype enhances absorption and carrier generation in the near-infrared to visible light range. Tuning the active layer's thickness increases optical activity at wavelengths greater than 600 nm, achieving ~190% external quantum efficiency over a broad solar wavelength range, underscoring its potential in advanced photovoltaic technology.

When Polymorphism in Metal-Organic Frameworks Enables Water Sorption Profile Tunability for Enhancing Heat Allocation and Water Harvesting Performance

The development of thermally driven water-sorption-based technologies relies on high-performing water vapor adsorbents. Here, polymorphism in Al-metal-organic frameworks is disclosed as a new strategy to tune the hydrophilicity of MOFs. This involves the formation of MOFs built from chains of either trans- or cis- µ-OH-connected corner-sharing AlO4(OH)2 octahedra. Specifically, [Al(OH)(muc)] or MIP-211, is made of trans, trans-muconate linkers, and cis-µ-OH-connected corner-sharing AlO4(OH)2 octahedra giving a 3D network with sinusoidal channels. The polymorph MIL-53-muc has a tiny change in the chain structure that results in a shift of the step position of the water isotherm from P/P0 ≈ 0.5 in MIL-53-muc, to P/P0 ≈ 0.3 in MIP-211. Solid-state NMR and Grand Canonical Monte Carlo reveal that the adsorption occurs initially between two hydroxyl groups of the chains, favored by the cis-positioning in MIP-211, resulting in a more hydrophilic behavior. Finally, theoretical evaluations show that MIP-211 would allow achieving a coefficient of performance for cooling (COPc) of 0.63 with an ultralow driving temperature of 60 °C, outperforming benchmark sorbents for small temperature lifts. Combined with its high stability, easy regeneration, huge water uptake capacity, green synthesis, MIP-211 is among the best adsorbents for adsorption-driven air conditioning and water harvesting from the air.

First-principles study of lithium aluminosilicate glass scintillators

Radiation sensors are an important enabling technology in several fields, such as medicine, scientific research, energy, defence, meteorology, and homeland security. Glass-based scintillators have been in use for more than 50 years and offer many benefits, including their ability to respond to different types of radiation, and to be readily formed into various shapes. There is, however, the prospect to develop new and improved glass scintillators, with low self-absorption, low refractive indices, and high radiative recombination rates. To investigate the factors limiting the improvement of glass scintillator properties, this work provides insight from atomic scale simulations of the cerium-doped lithium aluminosilicate (SiO2-Al2O3-MgO-Li2O-Ce2O3) glass scintillator system. Three glass compositions were studied using molecular dynamics and density functional theory to investigate the effect of the ratio (with RAl/M = [0.1, 0.8 and 1.2]) on the structural and electronic properties. For a ratio RAl/M > 1, it has been shown that glasses with increased polymerization allow for more effective incorporation of Ce3+ cations. The structural analysis also showed that the bond order of Al-O can be affected in the presence of a lithium-rich environment. Electronic density of states and Bader charge analysis indicate a decline in the population of localized trapping states with increasing RAl/M. This suggests a higher probability of radiative recombination which can increase the photon yield of these scintillators. These findings provide valuable guidance for optimizing Li-glasses in neutron detection systems by highlighting the intricate challenges.

Synthesis, characterization, and density functional theory investigation of (CH6N3)2[NpO2Cl3] and Rb[NpO2Cl2(H2O)] chain structures

The actinyl tetrachloro complex [An(V/VI)O2Cl4]2-/3- tends to form discrete molecular units in both solution and solid state materials, but related aquachloro complexes have been observed as both discrete coordination compounds and 1-D chain topologies. Subtle differences in the inner sphere coordination significantly influence the formation of structural topologies in the actinyl chloride system, but the exact reasoning for these variations has not been delineated. In the current study, we present the synthesis, structural characterization, and vibrational analysis of two 1-D neptunyl(V) chain compounds: (CH6N3)2[NpO2Cl3] (Np-Gua) and Rb[NpO2Cl2(H2O)] (Np-Rb). Bonding and non-covalent interactions (NCIs) in the systems were evaluated using periodic Density Functional Theory (DFT) to link these properties to related phases. We observed ∼6.5% and ∼3.9% weakening of NpO bonds in Np-Gua and Np-Rb compared to the reference Cs3[NpO2Cl4]. NCI analysis distinguished specific assembly modes, where Np-Gua was connected via hydrogen bonding (N-H⋯Cleq and N-H⋯Oyl) and Np-Rb contained both cation interactions (Rb+⋯Oyl and Rb+⋯Cleq) and hydrogen bonding (Oeq-H⋯Oyl) networks. Thermodynamically viable formation pathways for both compounds were explored using DFT methodology. The [NpO2Cl4](aq)3- and [NpO2Cl3(H2O)](aq)2- substructures were identified as precursors to Np-Gua and [NpO2Cl3(H2O)](aq)2- and [NpO2Cl2(H2O)2](aq)- were isolated as the primary building units of Np-Rb. Finally, we utilized DFT to analyze the vibrational modes for Np-Gua and Np-Rb, where we found evidence of the NpO bond weakening within the Np(V) chain structures compared to [NpO2Cl4]3-.

Understanding Carbon Monoxide Binding and Interactions in M-MOF-74 (M = Mg, Mn, Ni, Zn)

Small molecules may adsorb strongly in metal-organic frameworks (MOFs) through interactions with under-coordinated open metal sites (OMS) that often exist within these structures. Among adsorbates, CO is attractive to study both for its relevance in energy-related applications and for its ability to engage in both σ-donation and π-backbonding interactions with the OMS in MOFs. Concomitant with strong adsorption, structural changes arise due to modifications of the electronic structure of both the adsorbate and adsorbent. These structural changes affect the separation performance of materials, and accurately capturing these changes and the resulting energetics is critical for accurate predictive modeling of adsorption. Traditional approaches to modeling using classical force fields typically do not capture or account for changes at the electronic level. To characterize the structural and energetic effects of the local structural changes, we employed density functional theory (DFT) to study CO adsorption in M-MOF-74s. M-MOF-74s feature OMS at which CO is known to adsorb strongly and can be synthesized with a variety of divalent metal cations with distinct performance in adsorption. We considered M-MOF-74s with a range of metals of varied d-band occupations (Mg (3d0), Mn (3d5), Ni (3d8), and Zn (3d10)) with various structural constraints ranging from geometrically constrained adsorbent and adsorbate ions to fully optimized geometries to deconvolute the relative contributions of various structural effects to the adsorption energetics and binding distances observed. Our data indicate that the most significant structural changes during adsorption correlate with the greatest π-backbonding behaviors and commensurately result in a sizable binding energy change observed for CO adsorption. The insights built from this work are relevant to two longstanding research challenges within the MOF community: rational design of materials for separations and the design of force fields capable of accurately modeling adsorption.

Purely single-bonded spiral nitrogen chains stabilized by trivalent lanthanum ions

Inspired by the single-bonded nitrogen chains stabilized by tetravalent cerium, pentavalent tantalum, and hexavalent tungsten atoms, we explored the possibility of single-bonded nitrogen polymorphs stabilized by trivalent lanthanum ions. To achieve this, we utilized the crystal structure search method on the phase diagram of binary La-N compounds. We identified three novel thermodynamically stable phases, the C2/c LaN3, P-1 LaN4, and P-1 LaN8. Among them, the C2/c phase with infinite helical poly-N6 chains becomes thermodynamically stable above 50 GPa. Each nitrogen atom in the poly-N6 chain acquires one extra electron, and the spiral chain is purely single-bonded. The C2/c phase has an indirect band gap of ∼1.6 eV at 60 GPa. Notably, the band gap exhibits non-monotonic behavior, decreases first and then increases with increasing pressure. This abnormal behavior is attributed to the significant bonding of two La-N bonds at around 35 GPa. Phonon spectrum calculations and AIMD simulations have confirmed that the C2/c phase can be quenched to ambient conditions with slight distortion, and it exhibits excellent detonation properties. Additionally, we also discovered armchair-like nitrogen chains in LaN4 and the armchair and zigzag-like mixed nitrogen chains in LaN8. These results provide valuable insights into the electronic and bonding properties of nitrides under high pressure and may have important implications for the design and development of novel functional materials.

Two-Dimensional Noble Metal Chalcogenides in the Frustrated Snub-Square Lattice

We study two-dimensional noble metal chalcogenides, with compositions {Cu, Ag, Au}2{S, Se, Te}, crystallizing in a snub-square lattice. This is a semiregular two-dimensional tesselation formed by triangles and squares that exhibits geometrical frustration. We use for comparison a square lattice, from which the snub-square tiling can be derived by a simple rotation of the squares. The monolayer snub-square chalcogenides are very close to thermodynamic stability, with the most stable system (Ag2Se) a mere 7 meV/atom above the convex hull of stability. All compounds studied in the square and snub-square lattice are semiconductors, with band gaps ranging from 0.1 to more than 2.5 eV. Excitonic effects are strong, with an exciton binding energy of around 0.3 eV. We propose the Cu (001) surface as a possible substrate to synthesize Cu2Se, although many other metal and semiconducting surfaces can be found with very good lattice matching.

Band edges positions prediction of the of Ag nanocluster-decorated titania surfaces and their relationship to NO and NO2 interaction from first-principles calculations

Metal nanoclusters deposited on oxides have been widely used in photocatalysis playing an important role in the design of model catalysts with applications in heterogeneous catalysis. In particular, we are interested in the potential activity of these cluster-supported systems for the removal of nitrogen oxides either by possible catalytic reduction and/or by their adsorption. In this work, using first-principles methods, we evaluate the main characteristics of Agn (n = 1-4) nanoclusters isolated and deposited on anatase TiO2(101) and rutile TiO2(110) surfaces. Our results indicate that they are preferably adsorbed on rutile surface. The different formation energy at each surface can be explained using a Bader charge analysis. Particularly for Ag4 the lowest formation energy is obtained for tetrahedral geometry while the isolated Ag4 geometry is planar. Small silver deposits placed superficially on titania surfaces modify its electronic structures and improve the conduction band edges positions for possible NO reduction. Band edges positions with respect to the vacuum potential have been studied. The comparison of the conduction band minimum with the reduction potentials of NO/N2O and N2O/N2 shows that they are higher, being Ag3 on rutile and Ag1, Ag2 and Ag4P on anatase better for NO reduction. To complete the analysis, the calculation of work function, energy gap, ionization energy and electron affinity are relevant since they allow the location of semiconductor band edges at point of zero charge. Finally, the adsorption of nitrogen oxides is studied where the NO2 adsorption is favored over NO.

Graph Transformer with Convolution Parallel Networks for Predicting Single and Binary Component Adsorption Performance of Metal-Organic Frameworks

Metal-organic frameworks (MOFs) are considered one of the most important materials for carbon capture and storage (CCS) due to the advantages of porosity, multifunction, diverse structure, and controllable chemical composition. With the continuous development of artificial intelligence (AI) technology, more and more machine learning models are used to identify MOFs with high performance within a massive search space. However, current works have yet to form a model that uses graph-structured data only, which can predict the adsorption properties of single and binary components. In this work, we proposed and developed a graph transformer, combined with convolution parallel networks, called GC-Trans. The model can accurately and efficiently predict the adsorption performance of MOFs under the single- and binary-component adsorption conditions using only the features of the crystal diagram as inputs. By extracting and fusing local and global feature information, the model has stronger expression and generalization abilities. Thus, we used it to screen the ARC-MOF database and analyze the MOF structures that meet the target requirements. Additionally, to demonstrate the transferability of the model, we applied transfer learning methods to predict the CO2/CH4 separations and CH4 uptake, both of which showed good predictive performance.

Three-Dimensional Noncovalent Interaction Network within [NpO2Cl4]2- Coordination Compounds: Influence on Thermochemical and Vibrational Properties

Noncovalent interactions (NCIs) can influence the stability and chemical properties of pentavalent and hexavalent actinyl (AnO2+/2+) compounds. In this work, the impact of NCIs (actinyl-hydrogen and actinyl-cation interactions) on the enthalpy of formation (ΔHf) and vibrational features was evaluated using Np(VI) tetrachloro compounds as the model system. We calculated the ΔHf values of these solid-state compounds through density functional theory+ thermodynamics (DFT+ T) and validated the results against experimental ΔHf values obtained through isothermal acid calorimetry. Three structural descriptors were evaluated to develop predictors for ΔHf, finding a strong link between ΔHf and hydrogen bond energy (EHtotal) for neptunyl-hydrogen interactions and total electrostatic attraction energy (Eelectrostatictotal) for neptunyl-cation interactions. Finally, we used Raman spectroscopy together with bond order analysis to probe Np=O bond perturbation due to NCIs. Our results showed a strong correlation between the degree of NCIs by axial oxygen and red-shifting of Np=O symmetrical stretch (ν1) wavenumbers and quantitatively demonstrated that NCIs can weaken the Np=O bond. These properties were then compared to those of related U(VI) and Np(V) phases to evaluate the effects of subtle differences in the NCIs and overall properties. In general, the outcomes of our study demonstrated the role of NCIs in stabilizing actinyl solid materials, which consequently governs their thermochemical behaviors and vibrational signatures.

Local environmental engineering for highly stable single-atom Pt1/CeO2catalysts: first-principles insights

Single-atom Pt1/CeO2catalysts may cope with the high cost and durability issues of fuel cell electrocatalysts. In the present study, the stability and underlying interaction mechanisms of the Pt1/CeO2system are systematically investigated using first-principles calculations. The Pt adsorption energy on CeO2surfaces can be divided into chemical interaction and surface deformation parts. The interaction energy, mainly associated with the local chemical environment, i.e. the number of Pt-O bonds, plays a major role in Pt1/CeO2stability. When forming a Pt-4O configuration, the catalytic system has the highest stability and Pt is oxidized to Pt2+. An electronic metal-support interaction mechanism is proposed for understanding Pt1/CeO2stability. In addition, our calculations show that the Pt1/CeO2(100) system is dynamically stable, and the external O environment can promote the further oxidation of Pt to Ptn+(2 ≤n< 4). The present study provides useful guidance for the experimental development of highly stable and efficient electrocatalysts for fuel cell applications.

Discriminative Mechanical and Thermal Response of the H-N Bonds for the Energetic LLM-105 Molecular Assembly

Molecular interactions in energetic materials form the key not only to the "structure stability, energy storage, ignition, and detonation" dynamics but also to the sensitivity to the loading of perturbation and the power intensity of radiation for the energetic substance, with the nature of the interactions remaining elusive. With the aid of perturbative Raman spectroscopy and the pressure-resolved density functional theory, we uncovered that the H-N bond of the intermolecular O:H-N bonds for LLM-105 shares the same negative compressibility and thermal expansivity of the H-O bond for the coupling O:H-O bond of water [Phys. Rep. 2023, 998, 1-68]. In contrast, the dangling H-N bond vibrating at a 3440 cm-1 high frequency does otherwise due to the absence of coupling interaction and the undercoordination-driven bond contraction. These findings should deepen our insight into interactions involving electron lone pairs and offer an efficient means for discriminating the performance of individual bonds.

Synthesis, Characterization, and Density Functional Theory Investigation of the Solid-State [UO2Cl4(H2O)]2- Complex

A significant number of solid-state [UO2Cl4]2- coordination compounds have been synthesized and structurally characterized. Yet, despite their purposive relative abundance in aqueous solutions, characterization of aquachlorouranium(VI) complexes remain rare. In the current study, a solid-state uranyl aqua chloro complex ((C4H12N2)2[UO2Cl4(H2O)]Cl2) was synthesized using piperazinium as a charge-balancing ligand, and the structure was determined using single-crystal X-ray diffraction. Using periodic density functional theory, the electronic structure of the [UO2Cl4(H2O)]2- complex was compared to [UO2Cl4]2- to uncover the strengthening of the U═O bond in [UO2Cl4(H2O)]2-. Changes in the strength of the U═O bond were validated further with Raman and IR spectroscopy, where uranyl symmetrical (ν1) and asymmetrical (ν3) stretches were blue-shifted compared to the reference [UO2Cl4]2- complex. Furthermore, the formation energy of the solid-state (C4H12N2)2[UO2Cl4(H2O)]Cl2 complex was calculated to be -287.60 ± 1.75 kJ mol-1 using isothermal acid calorimetry. The demonstrated higher stability relative to the related [UO2Cl4]2- complex was related to the relative stoichiometry of the counterions.

Guiding Principles for the Rational Design of Hybrid Materials: Use of DFT Methodology for Evaluating Non-Covalent Interactions in a Uranyl Tetrahalide Model System

Together with the synthesis and experimental characterization of 14 hybrid materials containing [UO2 X4 ]2- (X=Cl- and Br- ) and organic cations, we report on novel methods for determining correlation trends in their formation enthalpy (ΔHf ) and observed vibrational signatures. ΔHf values were analyzed through isothermal acid calorimetry and a Density Functional Theory+Thermodynamics (DFT+T) approach with results showing good agreement between theory and experiment. Three factors (packing efficiency, cation protonation enthalpy, and hydrogen bonding energy [E H , norm total ${{E}_{H,{\rm { norm}}}^{{\rm { total}}}}$ ]) were assessed as descriptors for trends in ΔHf . Results demonstrated a strong correlation betweenE H , norm total ${E_{{\rm{H}},{\rm{norm}}}^{{\rm{total}}} }$ and ΔHf , highlighting the importance of hydrogen bonding networks in determining the relative stability of solid-state hybrid materials. Lastly, we investigate how hydrogen bonding networks affect the vibrational characteristics of uranyl solid-state materials using experimental Raman and IR spectroscopy and theoretical bond orders and find that hydrogen bonding can red-shift U≡O stretching modes. Overall, the tightly integrated experimental and theoretical studies presented here bridge the trends in macroscopic thermodynamic energies and spectroscopic features with molecular-level details of the geometry and electronic structure. This modeling framework forms a basis for exploring 3D hydrogen bonding as a tunable design feature in the pursuit of supramolecular materials by rational design.

The Influence of Ni Addition in the Mechanism of CO2 Electroreduction on Cu Crystals-Mechanistic Insight from DFT Simulations

We present a DFT analysis of the role of the Cu-Ni synergistic effect for the CO₂ reduction to C₂H₄, in comparison to the pure Cu catalyst. The analysis is focused on the thermodynamic stability of reactive intermediates along the proposed pathway of C₂ species formation. We have observed that the potential needed for the reaction decreases with the addition of Ni in the investigated model. In addition, we have observed the differences in the preferred pathway based on the significant differences in stability of the reactive intermediates depending on th Cu:Ni ratio. The results suggest that despite the fact the Cu surface is always exposed, and it is the only one that is able to directly interact with the intermediates, the presence of the Ni in the underlying sections of the crystal is significant enough to change the mechanism of the reaction.

Grain Size-Dependent Thermal Expansion of Nanocrystalline Metals

In the present work, we have used classical molecular dynamics and quantum mechanical density functional theory modeling to investigate the grain size-dependent thermal expansion coefficient (CTE) of nanocrystalline Cu. We find that the CTE increases by up to 20% with a gradually decreasing grain size. This behavior emerges as a result of the increased population of occupied anti-bonding states and bond order variation in the grain boundary regions, which contribute to the reduced resistance against thermally-induced bond stretching and dictate the thermal expansion behavior in the small grain size limit. As a part of the present work, we have established a procedure to produce ab initio thermal expansion maps that can be used for the prediction of the grain size-dependent CTE. This can serve as a modeling tool, e.g., to explore the impact of grain boundary impurity segregation on the CTE.

Theoretical investigation of CO2 capture in the MIL-88 series: effects of organic linker modification

CO2 capture is a crucial strategy to mitigate global warming and protect a sustainable environment. Metal-organic frameworks with large surface area, high flexibility, and reversible adsorption and desorption of gases are good candidates for CO2 capture. Among the synthesized metal-organic frameworks, the MIL-88 series has attracted our attention due to their excellent stability. However, a systematic investigation of CO2 capture in the MIL-88 series with different organic linkers is not available. Therefore, we clarified the topic via two sections: (1) elucidate physical insights into the CO2@MIL-88 interaction by van der Waals-dispersion correction density functional theory calculations, and (2) quantitatively study the CO2 capture capacity by grand canonical Monte Carlo simulations. We found that the 1πg, 2σu/1πu, and 2σg peaks of the CO2 molecule and the C and O p orbitals of the MIL-88 series are the predominant contributors to the CO2@MIL-88 interaction. The MIL-88 series, i.e., MIL-88A, B, C, and D, has the same metal oxide node but different organic linkers: fumarate (MIL-88A), 1,4-benzene-dicarboxylate (MIL-88B), 2,6-naphthalene-dicarboxylate (MIL-88C), and 4,4'-biphenyl-dicarboxylate (MIL-88D). The results exhibited that fumarate should be the best replacement for both the gravimetric and volumetric CO2 uptakes. We also pointed out a proportional relationship between the capture capacities with electronic properties and other parameters.

In silico characterization of nanoparticles

Nanoparticles (NPs) make for intriguing heterogeneous catalysts due to their large active surface area and excellent and often size-dependent catalytic properties that emerge from a multitude of chemically different surface reaction sites. NP catalysts are, in principle, also highly tunable: even small changes to the NP size or surface facet composition, doping with heteroatoms, or changes of the supporting material can significantly alter their physicochemical properties. Because synthesis of size- and shape-controlled NP catalysts is challenging, the ability to computationally predict the most favorable NP structures for a catalytic reaction of interest is an in-demand skill that can help accelerate and streamline the material optimization process. Fundamentally, simulations of NP model systems present unique challenges to computational scientists. Not only must considerable methodological hurdles be overcome in performing calculations with hundreds to thousands of atoms while retaining appropriate accuracy to be able to probe the desired properties. Also, the data generated by simulations of NPs are typically more complex than data from simulations of, for example, single crystal surface models, and therefore often require different data analysis strategies. To this end, the present work aims to review analytical methods and data analysis strategies that have proven useful in extracting thermodynamic trends from NP simulations.

The Stiffest and Strongest Predicted Material: C2 N Atomic Chains Approach the Theoretical Limits

Though linear atomic chains exhibit extreme properties, it is presently unclear how these properties can be maximized by the choice of elemental composition. Considering that boron, carbon, and nitrogen can form high modulus and high strength atomic chains, here an algorithm is developed to construct 143 possible atomic chains of these elements with 6 or fewer atoms in the primitive cell and explore their stabilities and mechanical properties by first-principles calculations. It is found that the gravimetric modulus (1032 GPa g^-1 cm³ ) and strength (108 GPa g^-1 cm³ ) of the C₂ N chain significantly exceed those of any known material, including the previously stiffest predicted material (C chain, 945 GPa g^-1 cm³ ) and the previously strongest predicted material (BC chain, 92 GPa g^-1 cm³ ), and also approach the theoretical limits of gravimetric modulus (1036 GPa g^-1 cm³ ) and strength (130 GPa g^-1 cm³ ). Mechanistic analyses demonstrate that the higher gravimetric modulus and strength of the C₂ N chain, compared with the C and BC chains, originate from its short, stiff chemical bonding and the abnormal decrease in bond length alternation under tension. The likely ease of fabrication and potential synthesis routes for C₂ N chains are discussed.

CRAFTED: An exploratory database of simulated adsorption isotherms of metal-organic frameworks

Grand Canonical Monte Carlo is an important method for performing molecular-level simulations and assisting the study and development of nanoporous materials for gas capture applications. These simulations are based on the use of force fields and partial charges to model the interaction between the adsorbent molecules and the solid framework. The choice of the force field parameters and partial charges can significantly impact the results obtained, however, there are very few databases available to support a comprehensive impact evaluation. Here, we present a database of simulations of CO₂ and N₂ adsorption isotherms on 690 metal-organic frameworks taken from the CoRE MOF 2014 database. We performed simulations with two force fields (UFF and DREIDING), six partial charge schemes (no charges, Qeq, EQeq, MPNN, PACMOF, and DDEC), and three temperatures (273, 298, 323 K). The resulting isotherms compose the Charge-dependent, Reproducible, Accessible, Forcefield-dependent, and Temperature-dependent Exploratory Database (CRAFTED) of adsorption isotherms.

Understanding the Reactivity of Supported Late Transition Metals on a Bare Anatase (101) Surface: A Periodic Conceptual DFT Investigation

The rapidly growing interest for new heterogeneous catalytic systems providing high atomic efficiency along with high stability and reactivity triggered an impressive progress in the field of single-atom catalysis. Nevertheless, unravelling the factors governing the interaction strength between the support and the adsorbed metal atoms remains a major challenge. Based on periodic density functional theory (DFT) calculations, this paper provides insight into the adsorption of single late transition metals on a defect-free anatase surface. The obtained adsorption energies fluctuate, with the exception of Pd, between -3.11 and -3.80 eV and are indicative of a strong interaction. Depending on the considered transition metal, we could attribute the strength of this interaction with the support to i) an electron transfer towards anatase (Ru, Rh, Ni), ii) s-d orbital hybridisation effects (Pt), or iii) a synergistic effect between both factors (Fe, Co, Os, Ir). The driving forces behind the adsorption were also found to be strongly related to Klechkowsky's rule for orbital filling. In contrast, the deviating behaviour of Pd is most likely associated with the lower dissociation enthalpy of the Pd-O bond. Additionally, the reactivity of these systems was evaluated using the Fermi weighted density of states approach. The resulting softness values can be clearly related to the electron configuration of the catalytic systems as well as with the net charge on the transition metal. Finally, these indices were used to construct a model that predicts the adsorption strength of CO on these anatase-supported d-metal atoms. The values obtained from this regression model show, within a 95 % probability interval, a correlation of 84 % with the explicitly calculated CO adsorption energies.

Structure and size-dependent vibrational and thermal properties of Ni clusters: A systematic ab initio approach

There is scarce information on the vibrational and thermal properties of small Ni clusters. Here, the outcomes of ab initio spin-polarized density functional theory calculations on the size and geometry effects upon the vibrational and thermal properties of Ni_n (n = 13 and 55) clusters, are discussed. For theses clusters a comparison is presented between the closed shell symmetric octahedral (O_h) and the icosahedral (I_h) geometries. The results indicate that the I_h isomers are lower in energy. Besides, ab initio molecular dynamics runs at T = 300K show that Ni₁₃ and Ni₅₅ clusters transform from their initial O_h geometries towards the corresponding I_h ones. For Ni₁₃, we also consider the lowest energy less symmetric layered 1-3-6-3 structure, and the cuboid, recently observed experimentally for Pt₁₃, which is competitive in energy but is unstable, as phonon analysis reveals. We calculate their vibrational density of states (νDOS) and heat capacity, and compare with the Ni FCC bulk counterpart. The characteristic features of the νDOS curves of these clusters are interpreted in terms of the clusters' sizes, the interatomic distance contractions, the bond order values as well as the internal pressure and strains of the clusters. We find that the softest possible frequency of the clusters is size and structure-dependent, being the smallest for the O_h ones. We identify mostly shear, tangential type displacements involving mainly surface atoms for the lowest frequency of the spectra of both I_h and O_h isomers. For the maximum frequencies of these clusters the central atom shows anti-phase movements against groups of nearest neighbor atoms. An excess of heat capacity at low temperatures with respect to the bulk is found, while at high temperatures a constant limiting value, close but lower to the Dulong and Petit value, is determined.

Carbon Dioxide Reduction on Transition Metal Dichalcogenides with Ni and Cu Edge Doping: A Density-Functional Theory Study

Transition-metal dichalcogenides (TMDs) have promising properties for their use as catalysts of CO₂ reduction to methane via the Sabatier reaction. In this article we use density-functional theory calculations to gain insight into the energetics of this reaction for Mo/W-based and S/Se-based TMDs with non-, Ni- and Cu-doping. We show that sulfur-based TMDs with Ni/Cu doping exhibit better indicators for catalytic performance of the CO₂ reduction reaction than non-doped and doped TMDs without active sites. In addition, the role of the transition metal was found to a much smaller influence in the reaction than the role of the chalcogen and dopant atoms, which influence the bonding strength and type, respectively.

Giant Tunability of Charge Transport in 2D Inorganic Molecular Crystals by Pressure Engineering

The unique intermolecular van der Waals force in emerging two-dimensional inorganic molecular crystals (2DIMCs) endows them with highly tunable structures and properties upon applying external stimuli. Using high pressure to modulate the intermolecular bonding, here we reveal the highly tunable charge transport behavior in 2DIMCs for the first time, from an insulator to a semiconductor. As pressure increases, 2D α-Sb₂ O₃ molecular crystal undergoes three isostructural transitions, and the intermolecular bonding enhances gradually, which results in a considerably decreased band gap by 25 % and a greatly enhanced charge transport. Impressively, the in situ resistivity measurement of the α-Sb₂ O₃ flake shows a sharp drop by 5 orders of magnitude in 0-3.2 GPa. This work sheds new light on the manipulation of charge transport in 2DIMCs and is of great significance for promoting the fundamental understanding and potential applications of 2DIMCs in advanced modern technologies.

Diantimony Complexes [CpR 2 Mo2 (CO)4 (μ,η2 -Sb2 )] (CpR =C5 H5 , C5 H4 t Bu) as Unexpected Ligands Stabilizing Silver(I)n (n=1-4) Monomers, Dimers and Chains

Synthesis and reactivity of transition metal compounds bearing "naked" pnictogen atoms is an active research area with remarkable bonding patterns observed in the formed compounds. Within this field, intense investigations on the coordination behavior of complexes possessing Pn and Asn (2≤n≤5) moieties have been conducted. However, studies on heavier analogues have been ignored so far due to arduous challenges related to low yields and moderate air stability. Herein, we present the first in-depth study addressing the reactivity of organometallic complexes containing Sb-donor atoms with several AgI salts. These reactions afforded twelve unprecedented aggregates as monomers, dimers as well as three- and four-membered chains of AgI ions claimed in the literature to be inaccessible. Interatomic distances as well as computational evidence obtained with help of several different methods suggest the presence of Ag⋅⋅⋅Ag interactions in all complexes containing more than one AgI ion.

Plasmonic Ag/Cu/PEG nanofluids prepared when solids meet liquids in the gas phase

Since the time of Faraday's experiments, the optical response of plasmonic nanofluids has been tailored by the shape, size, concentration, and material of nanoparticles (NPs), or by mixing different types of NPs. To date, water-based liquids have been the most extensively investigated host media, while polymers, such as poly(ethylene glycol) (PEG), have frequently been added to introduce repulsive steric interactions and protect NPs from agglomeration. Here, we introduce an inverse system of non-aqueous nanofluids, in which Ag and Cu NPs are dispersed in PEG (400 g mol^-1), with no solvents or chemicals involved. Our single-step approach comprises the synthesis of metal NPs in the gas phase using sputtering-based gas aggregation cluster sources, gas flow transport of NPs, and their deposition (optionally simultaneous) on the PEG surface. Using computational fluid dynamics simulations, we show that NPs diffuse into PEG at an average velocity of the diffusion front of the order of μm s^-1, which is sufficient for efficient loading of the entire polymer bulk. We synthesize yellow Ag/PEG, green Cu/PEG, and blue Ag/Cu/PEG nanofluids, in which the color is given by the position of the plasmon resonance. NPs are prone to partial agglomeration and sedimentation, with a slower kinetics for Cu. Density functional theory calculations combined with UV-vis data and zeta-potential measurements prove that the surface oxidation to Cu₂O and stronger electrostatic repulsion are responsible for the higher stability of Cu NPs. Adopting the De Gennes formalism, we estimate that PEG molecules adsorb on the NP surface in mushroom coordination, with the thickness of the adsorbed layer L < 1.4 nm, grafting density σ < 0.20, and the average distance between the grafted chains D > 0.8 nm. Such values provide sufficient steric barriers to retard, but not completely prevent, agglomeration. Overall, our approach offers an excellent platform for fundamental research on non-aqueous nanofluids, with metal-polymer and metal-metal interactions unperturbed by the presence of solvents or chemical residues.

Unexpected polarization properties of sub-nanosized magnesium clusters

The isotropic electrostatic polarizability (IEP) of sub-nanosized magnesium clusters Mg2-Mg32 was studied in an extensive set comprising 1237 structurally unique isomers. These isomers were found in the course of the global search for the potential energy surface minima of the magnesium clusters at the BP86/6-31G(d) level. The calculation of the polarizability at the same DFT level reveals an unexpected property of the IEP: the linear correlation between the polarizability of the most favorable isomers (and only them) and the cluster nuclearity n. Moreover, for each n, the most stable cluster isomer demonstrates nearly minimal IEP value among all found isomers of a given nuclearity. Surprisingly, these observed features are independent of the cluster structures which are quite different. We hypothesize that the energetic favorability of a cluster structure is related to their low polarizability. Apparently, the atoms forming the cluster tend to arrange themselves in such a way as to provide the most compact distribution of the cluster electron density. A possible explanation of the observed trends, their significance for cluster structure prediction, and the practical applications are discussed.

Discovery of Clustered-P1 Borophene and Its Application as the Lightest High-Performance Transistor

The two-dimensional network of boron atoms (borophene) has attracted attention for its ultralow molar mass and remarkable polymorphism. Synthesized polymorphs of borophene (striped, β₁₂, χ₃, and honeycomb), so far, are all found to be metallic. Employing a genetic algorithm-based structure searching technique, here we discover an allotrope, clustered-P1, which is located very close to the global energy minimum. Clustered-P1 exhibits a bulk silicon-like band gap (1.08 eV) with symmetric effective masses (∼0.2 m₀) for electrons and holes along the transport direction. Phonon dispersion and beyond room-temperature ab initio molecular dynamics studies further confirm its excellent dynamic and structural stability. Since two-dimensional semiconductors are promising silicon alternatives for complementary metal-oxide semiconductor (CMOS) technology extension, we further investigate the characteristics of clustered-P1-based transistors using self-consistent quantum transport models for channel lengths of 10-3 nm. The performance of these devices has been found to be balanced for p- and n-type transistors and meets the requirements of the International Roadmap for Devices and Systems (IRDS). Our study may aid in the experimental realization of the lightest high-performance transistor.

Conversion between Metavalent and Covalent Bond in Metastable Superlattices Composed of 2D and 3D Sublayers

Reversible conversion over multimillion times in bond types between metavalent and covalent bonds becomes one of the most promising bases for universal memory. As the conversions have been found in metastable states, an extended category of crystal structures from stable states via redistribution of vacancies, research on kinetic behavior of the vacancies is highly in demand. However, it remains lacking due to difficulties with experimental analysis. Herein, the direct observation of the evolution of chemical states of vacancies clarifies the behavior by combining analysis on charge density distribution, electrical conductivity, and crystal structures. Site-switching of vacancies of Sb₂Te₃ gradually occurs with diverged energy barriers owing to their own activation code: the accumulation of vacancies triggers spontaneous gliding along atomic planes to relieve electrostatic repulsion. Studies on the behavior can be further applied to multiphase superlattices composed of Sb₂Te₃ (2D) and GeTe (3D) sublayers, which represent superior memory performances, but their operating mechanisms were still under debate due to their complexity. The site-switching is favorable (suppressed) when Te-Te bonds are formed as physisorption (chemisorption) over the interface between Sb₂Te₃ (2D) and GeTe (3D) sublayers driven by configurational entropic gain (electrostatic enthalpic loss). Depending on the type of interfaces between sublayers, phases of the superlattices are classified into metastable and stable states, where the conversion could only be achieved in the metastable state. From this comprehensive understanding on the operating mechanism via kinetic behaviors of vacancies and the metastability, further studies toward vacancy engineering are expected in versatile materials.

Predicting Frequency from the External Chemical Environment: OH Vibrations on Hydrated and Hydroxylated Surfaces

Robust correlation curves are essential to decipher structural information from IR-vibrational spectra. However, for surface-adsorbed water and hydroxides, few such correlations have been presented in the literature. In this paper, OH vibrational frequencies are correlated against 12 structural descriptors representing the quantum mechanical or geometrical environment, focusing on those external to the vibrating molecule. A nonbiased fitting procedure based on Gaussian process regression (GPR) was used alongside simple analytical functional forms. The training data consist of 217 structurally unique OH groups from 38 water/metal oxide interface systems for MgO, CaO and CeO₂, all optimized at the DFT level, and the fully anharmonic and uncoupled OH vibrational signatures were calculated. Among our results, we find the following: (i) The intermolecular R(H···O) hydrogen bond distance is particularly strong, indicating the primary cause of the frequency shift. (ii) Similarly, the electric field along the H-bond vector is also a good descriptor. (iii) Highly detailed machine learning descriptors (ACSF, SOAP) are less intuitive but were found to be more capable descriptors. (iv) Combinations of geometric and QM descriptors give the best predictions, supplying complementary information.

Ab initio calculations of CO2 adsorption on β-C2S(100) and M3-C3S(001) surfaces: An exploration of early CO2 sequestration pathways

Investigating CO₂ sequestration in cement-based materials is significant for achieving carbon neutrality in the cement and concrete industries. The early CO₂ sequestration pathways on cement-based materials are fundamental for CO₂ sequestration, which is not clear. Towards this, the adsorption behavior of CO₂ on β-C₂S(100) and M3-C₃S(001) was investigated at the atomic level using density functional theory calculations, which were then compared with water adsorption results. The molecular adsorption configurations of CO₂ on both β-C₂S(100) and M3-C₃S(001) were tilted from their initial configurations due to the influence of surface Ca and O atoms. The CO₂ adsorption energy on M3-C₃S(001) and β-C₂S(100) were -0.458 eV and -0.426 eV, respectively, indicating adsorption on M3-C₃S(001) was more energetically favorable. After CO₂ adsorption, electrons were transferred from the surface to the CO₂ molecule. Furthermore, the Ca-O bond orders of β-C₂S(100) and M3-C₃S(001) after CO₂ adsorption were maximally decreased by 2.79% and 6.99%, respectively. A more significant adsorption influence on surfaces was found for H₂O, with more negative adsorption energy, more evident electron transfer, and a greater decrease in bond order. The CO₂ adsorption on β-C₂S(100) and M3-C₃S(001) were still spontaneous at 298 K and 1 atm. This study provides important theoretical insights into early CO₂ sequestration at the atomic level, which has practical implications for the design of efficient CO₂ sequestration technologies.

Water structure, dynamics and reactivity on a TiO2-nanoparticle surface: new insights from ab initio molecular dynamics

Water structure, dynamics and reactivity at the surface of a small TiO₂-nanoparticle fully immersed in water was investigated by an ab initio molecular dynamics simulation. Several modes of water binding were identified by assigning each atom to an atom type, representing a distinct chemical environment in the ab initio ensemble, and then computing radial distribution functions between the atom types. Surface reactivity was investigated by monitoring how populations of atom types change during the simulation. In order to acquire further insight, electron densities for a set of representative system snapshots were analyzed using an atoms-in-molecules approach. Our results reveal that water dissociation, where a water molecule splits at a bridging oxygen site to form a hydroxyl group and a protonated oxygen bridge, can occur by a mechanism involving transfer of a proton over several water molecules. The hydroxyl group and protonated oxygen bridge formed in the process persist (on a 10 ps time scale) and the hydroxyl group undergoes exchange using a mechanism similar to the one responsible for water dissociation. Rotational and translational dynamics of water molecules around the nanoparticle were analyzed in terms of reorientational time correlation functions and mean square displacement. While reorientation of water O-H vectors decreases quickly in the proximity of the nanoparticle surface, translational diffusion slows down more gradually. Our results give new insight into water structure, dynamics and reactivity on TiO₂-nanoparticle surfaces and suggest that water dissociation on curved TiO₂-nanoparticle surfaces can occur via more complex mechanisms than those previously identified for flat defect-free surfaces.