Incorporating solvent effects in DFT: insights from cation exchange in faujasites.

Zeolites are versatile materials renowned for their extra-framework cation exchange capabilities, with applications spanning diverse fields, including nuclear waste treatment. While detailed experimental characterization offers valuable insight, density functional theory (DFT) proves particularly adept at investigating ion exchange in zeolites, owing to its atomic and electronic resolution. However, the prevalent occurrence of zeolitic ion exchange in aqueous environments poses a challenge to conventional DFT modeling, traditionally conducted in a vacuum. This study seeks to enhance zeolite modeling by systematically evaluating predictive differences across varying degrees of aqueous solvent inclusion. Specifically focusing on monovalent cation exchange in Na-X zeolites, we explore diverse modeling approaches. These range from simple dehydrated systems (representing bare reference states in vacuum) to more sophisticated models that incorporate aqueous solvent effects through explicit water molecules and/or a dielectric medium. Through comparative analysis of DFT and semi-empirical DFT approaches, along with their validation against experimental results, our findings underscore the necessity to concurrently consider explicit and implicit solvent effects for accurate prediction of zeolitic ionic exchange.

Friends or Foes: Fundamental Principles of Th-Organic Scaffold Chemistry Using Zr-Analogs as a Guide.

The fundamental interest in actinide chemistry, particularly for the development of thorium-based materials, is experiencing a renaissance owing to the recent and rapidly growing attention to fuel cycle reactors, radiological daughters for nuclear medicine, and efficient nuclear stockpile development. Herein, we uncover fundamental principles of thorium chemistry on the example of Th-based extended structures such as metal-organic frameworks in comparison with the discrete systems and zirconium extended analogs, demonstrating remarkable over two-and-half-year chemical stability of Th-based frameworks as a function of metal node connectivity, amount of defects, and conformational linker rigidity through comprehensive spectroscopic and crystallographic analysis as well as theoretical modeling. Despite exceptional chemical stability, we report the first example of studies focusing on the reactivity of the most chemically stable Th-based frameworks in comparison with the discrete Th-based systems such as metal-organic complexes and a cage, contrasting multicycle recyclability and selectivity (>97%) of the extended structures in comparison with the molecular compounds. Overall, the presented work not only establishes the conceptual foundation for evaluating the capabilities of Th-based materials but also represents a milestone for their multifaceted future and foreshadows their potential to shape the next era of actinide chemistry.

f-block MOFs: A Pathway to Heterometallic Transuranics.

A novel series of heterometallic f-block-frameworks including the first examples of transuranic heterometallic 238 U/239 Pu-metal-organic frameworks (MOFs) and a novel monometallic 239 Pu-analog are reported. In combination with theoretical calculations, we probed the kinetics and thermodynamics of heterometallic actinide(An)-MOF formation and reported the first value of a U-to-Th transmetallation rate. We concluded that formation of uranyl species could be a driving force for solid-state metathesis. Density of states near the Fermi edge, enthalpy of formation, band gap, proton affinity, and thermal/chemical stability were probed as a function of metal ratios. Furthermore, we achieved 97 % of the theoretical maximum capacity for An-integration. These studies shed light on fundamental aspects of actinide chemistry and also foreshadow avenues for the development of emerging classes of An-containing materials, including radioisotope thermoelectric generators or metalloradiopharmaceuticals.

Capture Instead of Release: Defect-Modulated Radionuclide Leaching Kinetics in Metal-Organic Frameworks.

Comparison of defect-controlled leaching-kinetics modulation of metal-organic frameworks (MOFs) and porous functionalized silica-based materials was performed on the example of a radionuclide and radionuclide surrogate for the first time, revealing an unprecedented readsorption phenomenon. On a series of zirconium-based MOFs as model systems, we demonstrated the ability to capture and retain >99% of the transuranic 241Am radionuclide after 1 week of storage. We report the possibility of tailoring radionuclide release kinetics in MOFs through framework defects as a function of postsynthetically installed organic ligands including cation-chelating crown ether-based linkers. Based on comprehensive analysis using spectroscopy (EXAFS, UV-vis, FTIR, and NMR), X-ray crystallography (single crystal and powder), and theoretical calculations (nine kinetics models and structure simulations), we demonstrated the synergy of radionuclide integration methods, topological restrictions, postsynthetic scaffold modification, and defect engineering. This combination is inaccessible in any other material and highlights the advantages of using well-defined frameworks for gaining fundamental knowledge necessary for the advancement of actinide-based material development, providing a pathway for addressing upcoming challenges in the nuclear waste administration sector.

Heterometallic Actinide-Containing Photoresponsive Metal-Organic Frameworks: Dynamic and Static Tuning of Electronic Properties.

Acquiring fundamental knowledge of properties of actinide-based materials is a necessary step to create new possibilities for addressing the current challenges in the nuclear energy and nuclear waste sectors. In this report, we established a photophysics-electronics correlation for actinide-containing metal-organic frameworks (An-MOFs) as a function of excitation wavelength, for the first time. A stepwise approach for dynamically modulating electronic properties was applied for the first time towards actinide-based heterometallic MOFs through integration of photochromic linkers. Optical cycling, modeling of density of states near the Fermi edge, conductivity measurements, and photoisomerization kinetics were employed to shed light on the process of tailoring optoelectronic properties of An-MOFs. Furthermore, the first photochromic MOF-based field-effect transistor, in which the field-effect response could be changed through light exposure, was constructed. As a demonstration, the change in current upon light exposure was sufficient to operate a two-LED fail-safe indicator circuit.

Heterometallic multinuclear nodes directing MOF electronic behavior.

Metal node engineering in combination with modularity, topological diversity, and porosity of metal-organic frameworks (MOFs) could advance energy and optoelectronic sectors. In this study, we focus on MOFs with multinuclear heterometallic nodes for establishing metal-property trends, i.e., connecting atomic scale changes with macroscopic material properties by utilization of inductively coupled plasma mass spectrometry, conductivity measurements, X-ray photoelectron and diffuse reflectance spectroscopies, and density functional theory calculations. The results of Bader charge analysis and studies employing the Voronoi-Dirichlet partition of crystal structures are also presented. As an example of frameworks with different nodal arrangements, we have chosen MOFs with mononuclear, binuclear, and pentanuclear nodes, primarily consisting of first-row transition metals, that are incorporated in HHTP-, BTC-, and NIP-systems, respectively (HHTP3- = triphenylene-2,3,6,7,10,11-hexaone; BTC3- = 1,3,5-benzenetricarboxylate; and NIP2- = 5-nitroisophthalate). Through probing framework electronic profiles, we demonstrate structure-property relationships, and also highlight the necessity for both comprehensive analysis of trends in metal properties, and novel avenues for preparation of heterometallic multinuclear isoreticular structures, which are critical components for on-demand tailoring of properties in heterometallic systems.

Thermodynamics and Electronic Properties of Heterometallic Multinuclear Actinide-Containing Metal-Organic Frameworks with "Structural Memory".

Thermodynamic studies of actinide-containing metal-organic frameworks (An-MOFs), reported herein for the first time, are a step toward addressing challenges related to effective nuclear waste administration. In addition to An-MOF thermochemistry, enthalpies of formation were determined for the organic linkers, 2,2'-dimethylbiphenyl-4,4'-dicarboxylic acid (H2Me2BPDC) and biphenyl-4,4'-dicarboxylic acid (H2BPDC), which are commonly used building blocks for MOF preparation. The electronic structure of the first example of An-MOF with mixed-metal AnAn'-nodes was influenced through coordination of transition metals as shown by the density of states near the Fermi edge, changes in the Tauc plot, conductivity measurements, and theoretical calculations. The "structural memory" effect (i.e., solvent-directed crystalline-amorphous-crystalline structural dynamism) was demonstrated as a function of node coordination degree, which is the number of organic linkers per metal node. Remarkable three-month water stability was reported for Th-containing frameworks herein, and the mechanism is also considered for improvement of the behavior of a U-based framework in water. Mechanistic aspects of capping linker installation were highlighted through crystallographic characterization of the intermediate, and theoretical calculations of free energies of formation (ΔGf) for U- and Th-MOFs with 10- and 12-coordinated secondary building units (SBUs) were performed to elucidate experimentally observed transformations during the installation processes. Overall, these results are the first thermochemical, electronic, and mechanistic insights for a relatively young class of actinide-containing frameworks.

Galvanic Replacement-Driven Transformations of Atomically Intermixed Bimetallic Colloidal Nanocrystals: Effects of Compositional Stoichiometry and Structural Ordering.

Galvanic replacement reactions dictated by deliberately designed nanoparticulate templates have emerged as a robust and versatile approach that controllably transforms solid monometallic nanocrystals into a diverse set of architecturally more sophisticated multimetallic hollow nanostructures. The galvanic atomic exchange at the nanoparticle/liquid interfaces induces a series of intriguing structure-transforming processes that interplay over multiple time and length scales. Using colloidal Au-Cu alloy and intermetallic nanoparticles as structurally and compositionally fine-tunable bimetallic sacrificial templates, we show that atomically intermixed bimetallic nanocrystals undergo galvanic replacement-driven structural transformations remarkably more complicated than those of their monometallic counterparts. We interpret the versatile structure-transforming behaviors of the bimetallic nanocrystals in the context of a unified mechanistic picture that rigorously interprets the interplay of three key structure-evolutionary pathways: dealloying, Kirkendall diffusion, and Ostwald ripening. By deliberately tuning the compositional stoichiometry and atomic-level structural ordering of the Au-Cu bimetallic nanocrystals, we have been able to fine-maneuver the relative rates of dealloying and Kirkendall diffusion with respect to that of Ostwald ripening through which an entire family of architecturally distinct complex nanostructures are created in a selective and controllable manner upon galvanic replacement reactions. The insights gained from our systematic comparative studies form a central knowledge framework that allows us to fully understand how multiple classic effects and processes interplay within the confinement by a colloidal nanocrystal to synergistically guide the structural transformations of complex nanostructures at both the atomic and nanoparticulate levels.

Graphene-Titanium Interfaces from Molecular Dynamics Simulations.

Unraveling the physical and chemical properties of graphene-metal contacts is a key step toward the development of graphitic electronic nanodevices. Although many studies have revealed the way that various metals interact with graphene, few have described the structure and behavior of large pieces of graphene-metal nanostructures under different conditions. Here, we present the first classical molecular dynamics study of graphene-titanium (G-Ti) structures, with and without substrates. Physical and chemical properties of equilibrium structures of G-Ti interfaces with different amounts of titanium coverage are investigated. Adhesion of Ti films on graphene is shown to be enhanced by the vacancies in graphene or the electrostatic influence of substrates. The dynamics of pristine G-Ti structures at different temperatures on planar and nonplanar substrates are investigated, and the results show that G-Ti interfaces are thermally stable, that is, not prone to any reaction toward the formation of titanium carbide.

Effect of Surface Chemistry on Water Interaction with Cu(111).

The interfacial dynamics of water in contact with bare, oxidized, and hydroxylated copper surfaces are examined using classical molecular dynamics (MD) simulations. A third-generation charge-optimized many-body (COMB3) potential is used in the MD simulations to investigate the adsorption of water molecules on Cu(111), and the results are compared to the findings of density functional theory (DFT) calculations. The adsorption energies and structures predicted by COMB3 are generally consistent with those determined with DFT. The COMB3 potential is then used to investigate the wetting behavior of water nanodroplets on Cu(111) at 20, 130, and 300 K. At room temperature, the simulations predict that the spreading rate of the base radius, R0, of a water droplet with a diameter of about 1.5 nm exhibits a spreading rate of R0 ≈ t(0.16) and a final base radius of 3.5 nm. At 20 and 130 K, water droplets are predicted to retain their structure after adsorption on Cu(111) and to undergo minimal spreading in agreement with scanning tunneling microscopy data. When the same water droplet encounters a reconstructed, oxidized Cu(111) surface, the classical MD simulations predict wetting with a spreading rate of R ≈ t(0.14) and a final base radius of 3.0 nm. Similarly, our MD simulations predict a spreading rate of R ≈ t(0.14) and a final base radius of 2.5 nm when water encounters OH-covered Cu(111). These results indicate that oxidation and hydroxylation cause a reduction in the degree of spreading and final base radius that is directly associated with a decreased spreading rate for water nanodroplets on copper.

Charge optimized many-body (COMB) potential for dynamical simulation of Ni-Al phases.

An interatomic potential for the Ni-Al system is presented within the third-generation charge optimized many-body (COMB3) formalism. The potential has been optimized for Ni3Al, or the γ' phase in Ni-based superalloys. The formation energies predicted for other Ni-Al phases are in reasonable agreement with first-principles results. The potential further predicts good mechanical properties for Ni3Al, which includes the values of the complex stacking fault (CSF) and the anti-phase boundary (APB) energies for the (1 1 1) and (1 0 0) planes. It is also used to investigate dislocation propagation across the Ni3Al (1 1 0)-Ni (1 1 0) interface, and the results are consistent with simulation results reported in the literature. The potential is further used in combination with a recent COMB3 potential for Al2O3 to investigate the Ni3Al (1 1 1)-Al2O3 (0 0 01) interface, which has not been modeled previously at the classical atomistic level due to the lack of a reactive potential to describe both Ni3Al and Al2O3 as well as interactions between them. The calculated work of adhesion for this interface is predicted to be 1.85 J m(-2), which is in agreement with available experimental data. The predicted interlayer distance is further consistent with the available first-principles results for Ni (1 1 1)-Al2O3 (0 0 0 1).

Charge optimized many-body (COMB) potential for Al2O3 materials, interfaces, and nanostructures.

This work presents the development and applications of a new empirical, variable charge potential for Al2O3 systems within the charge optimized many-body (COMB) potential framework. The potential can describe the fundamental physical properties of Al2O3, including cohesive energy, elastic constants, defect formation energies, surface energies and phonon properties of α-Al2O3 comparable to that obtained from experiments and first-principles calculations. The potential is further employed in classical molecular dynamics (MD) simulations to validate and predict the properties of the Al (1 1 1)-Al2O3 (0 0 0 1) interface, tensile properties of Al nanowires, Al2O3 nanowires, Al2O3-covered Al nanowires, and defective Al2O3 nanowires. The results demonstrate that the potential is well-suited to model heterogeneous material systems involving Al and Al2O3. Most importantly, the parameters can be seamlessly coupled with COMB3 parameters for other materials to enable MD simulations of a wide range of heterogeneous material systems.

Charge optimized many-body potential for aluminum.

An interatomic potential for Al is developed within the third generation of the charge optimized many-body (COMB3) formalism. The database used for the parameterization of the potential consists of experimental data and the results of first-principles and quantum chemical calculations. The potential exhibits reasonable agreement with cohesive energy, lattice parameters, elastic constants, bulk and shear modulus, surface energies, stacking fault energies, point defect formation energies, and the phase order of metallic Al from experiments and density functional theory. In addition, the predicted phonon dispersion is in good agreement with the experimental data and first-principles calculations. Importantly for the prediction of the mechanical behavior, the unstable stacking fault energetics along the [Formula: see text] direction on the (1 1 1) plane are similar to those obtained from first-principles calculations. The polycrsytal when strained shows responses that are physical and the overall behavior is consistent with experimental observations.

A charge optimized many-body (COMB) potential for titanium and titania.

This work proposes an empirical, variable charge potential for Ti and TiO(2) systems based on the charge-optimized many-body (COMB) potential framework. The parameters of the potential function are fit to the structural and mechanical properties of the Ti hcp phase, the TiO(2) rutile phase, and the energetics of polymorphs of both Ti and TiO(2). The relative stabilities of TiO(2) rutile surfaces are predicted and compared to the results of density functional theory (DFT) and empirical potential calculations. The transferability of the developed potential is demonstrated by determining the adsorption energy of Cu clusters of various sizes on the rutile TiO(2)(1 1 0) surface using molecular dynamics simulations. The results indicate that the adsorption energy is dependent on the number of Cu-Cu bonds and Cu-O bonds formed at the Cu/TiO(2) interface. The adsorption energies of Cu clusters on the reduced and oxidized TiO(2)(1 1 0) surfaces are also investigated, and the COMB potential predicts enhanced bonding between Cu clusters and the oxidized surface, which is consistent with both experimental observations and the results of DFT calculations for other transition metals (Au and Ag) on this oxidized surface.

A charge-optimized many-body potential for the U-UO2-O2 system.

Building on previous charge-optimized many-body (COMB) potentials for metallic α-U and gaseous O2, we have developed a new potential for UO2, which also allows the simulation of U-UO2-O2 systems. The UO2 lattice parameter, elastic constants and formation energies of stoichiometric and non-stoichiometric intrinsic defects are well reproduced. Moreover, this is the first rigid-ion potential that produces the correct deviation of the Cauchy relation, as well as the first classical interatomic potential that is able to determine the defect energies of non-stoichiometric intrinsic point defects in UO2 with an appropriate reference state. The oxygen molecule interstitial in the α-U structure is shown to decompose, with some U-O bonds approaching the natural bond length of perfect UO2. Finally, we demonstrate the capability of this COMB potential to simulate a complex system by performing a simulation of the α-U + O2 â†’ UO2 phase transformation. We also identify a possible mechanism for uranium oxidation and the orientation of the resulting fluorite UO2 structure relative to the coordinate system of orthorhombic α-U.

Phonon lifetime investigation of anharmonicity and thermal conductivity of UO2 by neutron scattering and theory.

Inelastic neutron scattering measurements of individual phonon lifetimes and dispersion at 295 and 1200 K have been used to probe anharmonicity and thermal conductivity in UO2. They show that longitudinal optic phonon modes carry the largest amount of heat, in contrast to past simulations and that the total conductivity demonstrates a quantitative correspondence between microscopic and macroscopic phonon physics. We have further performed first-principles simulations for UO2 showing semiquantitative agreement with phonon lifetimes at 295 K, but larger anharmonicity than measured at 1200 K.

Variable charge reactive potential for hydrocarbons to simulate organic-copper interactions.

A variable charge reactive empirical potential for carbon-based materials, hydrocarbons, organometallics, and their interfaces is developed within the framework of charge optimized many-body (COMB) potentials. The resulting potential contains improved expressions for the bond order and self-energy, which gives a flexible, robust, and integrated treatment of different bond types in multicomponent and multifunctional systems. It furthermore captures the dissociation and formation of the chemical bonds and appropriately and dynamically determines the associated charge transfer, thus providing a powerful method to simulate the complex chemistry of many-atom systems in changing environments. The resulting COMB potential is used in a classical molecular dynamics simulation of the room temperature, low energy deposition of ethyl radicals on the Cu (111) surface (a system with ∼5000 atoms) to demonstrate its capabilities at describing organic-metal interactions in a dynamically changing environment.

Classical interatomic potential for orthorhombic uranium.

A classical interatomic potential for uranium metal is derived within the framework of the charge optimized many body (COMB) formalism. The potential is fitted with a database obtained from experiment and density functional theory (DFT) calculations. The potential correctly predicts orthorhombic α-U to be the ground state. Good agreement with experimental values is obtained for the lattice parameters, nearest neighbor distances, and elastic constants. Molecular dynamics simulations also correctly show the anisotropy in the coefficient of thermal expansion and the temperature dependence of the nearest neighbor distances.

Effect of fluorocarbon molecules confined between sliding self-mated PTFE surfaces.

We report on the frictional response and atomic process that occur when molecular fluorocarbon molecules of varying lengths are sheared between two polytetrafluoroethylene (PTFE) surfaces. The thicknesses of the molecular layers are also varied. The approach is classical molecular dynamics simulations using a reactive bond-order potential parametrized for fluorocarbons. The results indicate that the presence of the molecules has a significant impact on the measured friction and wear of the surfaces, and that this impact depends on the nature of the fluorocarbon molecules and the thickness of the molecular film. The molecular mechanisms responsible for these differences are presented.

Structure and energetics of 180° domain walls in PbTiO3 by density functional theory.

Density functional theory at the level of the local density approximation with the projector augmented wave method is used to determine the structure of 180° domain walls in tetragonal ferroelectric PbTiO(3). In agreement with previous studies, it is found that PbO-centered {100} walls have lower energies than TiO(2)-centered {100} walls, leading to a Peierls potential barrier for wall motion along <010> of ∼36 mJ m(-2). In addition to the Ising-like polarization along the tetragonal axis, it is found that near the domain wall, there is a small polarization in the wall-normal direction away from the domain wall. These Néel-like contributions to the domain wall are analyzed in terms of the Landau-Ginzburg-Devonshire phenomenological theory for ferroelectrics. Similar characteristics are found for {110} domain walls, where OO-centered walls are energetically more favorable than the PbTiO-centered walls.