Search

Optical control of polyatomic molecules promises new opportunities in precision metrology, fundamental chemistry, quantum information, and many-body science. Contemporary experimental and theoretical efforts have mostly focused on cycling photons via excitation of a single electron localized to an alkaline earth (group 2)-like metal center. In this manuscript, we consider pathways towards optical cycling in polyatomic molecules with multi-electron degrees of freedom, which arise from two or more cycling electrons localized to $p$-block post-transition metal and metalloid (group 13, 14, and 15) centers. We characterize the electronic structure and rovibrational branching of several prototypical candidates using ab initio quantum chemical methods. Despite increased internal complexity and challenging design parameters, we find several molecules possessing quasi-closed photon cycling schemes with highly diagonal, visible and near-infrared transitions. Furthermore, we identify new heuristics for engineering optically controllable and laser-coolable polyatomic molecules with multi-electron cycling centers. Our results help elucidate the interplay between hybridization, repulsion, and ionicity in optically active species and provide a first step towards using polyatomic molecules with complex electronic structure as a resource for quantum science and measurement., Comment: 29 pages, 11 figures, 5 tables, and 6 supplemental data tables

Transition metal impurities such as nickel, copper, and iron, in solid-state materials like silicon have a significant impact on the electrical performance of integrated circuits and solar cells. To study the impact of copper impurities inside bulk silicon on the electrical properties of the material, one needs to understand the configurational space of copper atoms incorporated inside the silicon lattice. In this work, we performed ReaxFF reactive force field based molecular dynamics simulations, studying different configurations of individual and crystalline copper atoms inside bulk silicon by looking at the diffusional behavior of copper in silicon. The ReaxFF Cu/Si parameter set was developed by training against DFT data, including the energy barrier for an individual Cu-atom inside a silicon lattice. We found that the diffusion of copper atoms has a direct relationship with the temperature. Moreover, it is also shown that individual copper atoms start to clusterize inside bulk silicon at elevated temperatures. Our simulation results provide a comprehensive picture of the effects of temperature and copper concentration on the crystallization of individual copper inside silicon lattice. Finally, the stress-strain relationship of Cu/Si compounds under uniaxial tensile loading have been obtained. Our results indicate a decrease in the elastic modulus with increasing level of Cu-impurity concentration. We observe spontaneous microcracking of the Si during the stress-strain tests as a consequence of the formation of a small Cu clusters adjacent to the Si surface., Comment: 29 pages, 13 figures

The amorphous propagation of microstructures during the electrochemical charging of a battery is the main reason for the capacity decay and short circuit. The charge distribution across the micro-structure is the result of both local and global equilibrium and is non-convex problem merely due to stochastic placement of the atoms. As such, obtaining the charge equilibrium (QEq) is a critical factor, since the amount of charge determines the success rate of the bond formation for the ionic species approaching the microstructure and consequently the ultimate morphology of the electrochemical deposits. Herein we develop a computationally-affordable method for determining the charge allocation within such microstructures. The cost function and the span of the charge distribution correlates very closely with the trivial method as well as a conventional method, albeit having significantly less computational cost. The method can be used for optimization in non-convex environments, specially those of stochastic nature., Comment: 20 pages, 8 Figures. arXiv admin note: text overlap with arXiv:1908.02362, arXiv:2006.13349

In this paper we develop a feed-back control framework for the real-time minimization of microstructures grown within the rechargeable battery. Due to quickening nature of the branched evolution, we identify the critical ramified peaks in the early stages and based on the state we compute the relaxation time for the concentration in those branching fingers. The control parameter is a function of the maximum curvature (i.e. minimum radius) of the branched microstructure, where the higher rate dendritic evolution would lead to the more critical state to be controlled. The charging time is minimized for generating the most packed microstructures and obtained results correlate closely with those of considerably higher charging time periods. The developed framework could be utilized as a smart charging protocol for the safe and sustainable operation the rechargeable batteries, where the branching of the microstructures could be correlated to the sudden variation in the current/voltage., Comment: 17 pages, 14 figures.arXiv admin note: text overlap with arXiv:1908.02362

The ramified and stochastic evolution of dendritic microstructures has been a major issue on the safety and longevity of rechargeable batteries, particularly for the utilization high-energy metallic electrodes. We analytically develop criteria for the pulse characteristics leading to the effective halting of the ramified electrodeposits grown during extensive time scales beyond inter-ionic collisions. Our framework is based on the competitive interplay between diffusion and electromigration and tracks the gradient of ionic concentration throughout the entire cycle of pulse-rest as a critical measure for heterogeneous evolution. In particular, the framework incorporates the Brownian motion of the ions and investigates the role of the geometry of the electrodeposition interface. Our novel experimental observations verify the analytical developments, where the the dimension-free developments allows the application to the electrochemical systems of various scales., Comment: 24 Pages, 10 Figures, 3 Tables

Percolation phenomena are pervasive in nature, ranging from capillary flow, crack propagation, ionic transport, fluid permeation, etc. Modeling percolation in highly-branched media requires the use of numerical solutions, as problems can quickly become intractable due to the number of pathways available. This becomes even more challenging in dynamic scenarios where the generation of pathways can quickly become a combinatorial problem. In this work, we develop a new constriction percolation paradigm, using cellular automata to predict the transport of oxygen through a stochastically cracked Zr oxide layer within a coupled diffusion-reaction framework. We simulate such branching trees by generating a series porosity-controlled media. Additionally, we develop an analytical criterion based on compressive yielding for bridging the transition state in corrosion regime, where the percolation threshold has been achieved. Our model extends Dijkstras shortest path method to constriction pathways and predicts the arrival rate of oxygen ions at the oxide interface. This is a critical parameter to predict oxide growth in the so-called post-transition regime, when bulk diffusion is no longer the rate-limiting phenomenon., Comment: 16 pages, 14 figures

We analyze from a far field the evolution of an interface that separates ideal incompressible fluids of different densities and has an interfacial mass flux. We develop and apply the general matrix method to rigorously solve the boundary value problem involving the governing equations in the fluid bulk and the boundary conditions at the interface and at the outside boundaries of the domain. We find the fundamental solutions for the linearized system of equations, and analyze the interplay of interface stability with flow fields structure, by directly linking rigorous mathematical attributes to physical observables. New mechanisms are identified of the interface stabilization and destabilization. We find that interfacial dynamics is stable when it conserves the fluxes of mass, momentum and energy. The stabilization is due to inertial effects causing small oscillations of the interface velocity. In the classic Landau dynamics, the postulate of perfect constancy of the interface velocity leads to the development of the Landau-Darrieus instability. This destabilization is also associated with the imbalance of the perturbed energy at the interface, in full consistency with the classic results. We identify extreme sensitivity of the interface dynamics to the interfacial boundary conditions, including formal properties of fundamental solutions and qualitative and quantitative properties of the flow fields. This provides new opportunities for studies, diagnostics, and control of multiphase flows in a broad range of processes in nature and technology.

The recent application of electrosprays to characterize the air-water interface, along with the reports on dramatically accelerated chemical reactions in aqueous electrosprays, have sparked a broad interest. Herein, we report on complementary laboratory and in silico experiments tracking the oligomerization of isoprene, an important biogenic gas, in electrosprays and isoprene-water emulsions to differentiate the contributions of interfacial effects from those of high voltages leading to charge-separation and concentration of reactants in the electrosprays. To this end, we employed electrospray ionization mass spectrometry, proton nuclear magnetic resonance, and quantum mechanical simulations. We found that the oligomerization of isoprene in aqueous electrosprays involved minimally hydrated and highly reactive hydronium ions. Those conditions, however, are non-existent at pristine air-water interfaces and oil-water emulsions under normal temperature and pressure. Thus, electrosprays should be complemented with surface-specific platforms and theoretical methods to reliably investigate chemistries at the pristine air-water interface., Comment: 16 pages and 6 figures in the main paper; 30 pages and 14 figures in total including supporting information

The transformation pathways of carbon at high pressures are of broad interest for synthesis of novel materials and for revealing the Earth's geological history. We have applied large plastic shear on graphite in rotational anvils to form hexagonal and nanocrystalline cubic diamond at extremely low pressures of 0.4 and 0.7 GPa, which are 50 and 100 times lower than the transformation pressures under hydrostatic compression and well below the phase equilibrium. Large shearing accompanied with pressure elevation to 3 GPa also leads to formation of a new orthorhombic diamond phase. Our results demonstrate new mechanisms and new means for plastic shear-controlled material synthesis at drastically reduced pressures, enabling new technologies for material synthesis. The results indicate that the micro-diamonds found in the low pressure-temperature crust could have formed during a large shear producing event, such as tectonic rifting and continued plate collision, without the need to postulate subduction to the mantle., Comment: 15 pages, 6 figures, 3 tables

Na is known to deliver very low energy capacity for sodium intercalation compared to Lithium. In this study, we use quantum mechanics based metadynamics simulations to obtain the free energy landscape for sodium ion intercalation from Dimethyl sulfoxide (DMSO) solvent into graphite. We find that the lowest free energy minima from the metadynamics are associated with sodium solvated by 3 or 4 DMSO. The free energy minima of these states are activated by a free energy of solvation computed to be 0.17 eV, which in turn are the most thermodynamically stable. We observe weak interactions of sodium with graphite sheets during the unbiased and biased molecular dynamics simulations. Our simulations results indicate that solvent plays an important role in stabilizing the sodium intercalation into graphite through shielding of the sodium while modulating the interaction of the solvent with the graphite sheets. In order to facilitate this intercalation, we propose solvents with negatively charged groups and aromatic cores (e.g., cyclic ethers) that could allow a greater rate of anion exchange to increase Na+ mobility., Comment: 15 pages, 4 figures

Plasmonic hot carrier devices extract excited carriers from metal nanostructures before equilibration, and have the potential to surpass semiconductor light absorbers. However their efficiencies have so far remained well below theoretical limits, which necessitates quantitative prediction of carrier transport and energy loss in plasmonic structures to identify and overcome bottlenecks in carrier harvesting. Here, we present a theoretical and computational framework, Non-Equilibrium Scattering in Space and Energy (NESSE), to predict the spatial evolution of carrier energy distributions that combines the best features of phase-space (Boltzmann) and particle-based (Monte Carlo) methods. Within the NESSE framework, we bridge first-principles electronic structure predictions of plasmon decay and carrier collision integrals at the atomic scale, with electromagnetic field simulations at the nano- to mesoscale. Finally, we apply NESSE to predict spatially-resolved energy distributions of photo-excited carriers that impact the surface of experimentally realizable plasmonic nanostructures at length scales ranging from tens to several hundreds of nanometers, enabling first-principles design of hot carrier devices., Comment: 10 pages, 5 figures

The reliable production of two-dimensional crystals is essential for the development of new technologies based on 2D materials. However, current synthesis methods suffer from a variety of drawbacks, including limitations in crystal size and stability. Here, we report the fabrication of large-area, high-quality 2D tellurium (tellurene) using a substrate-free solution process. Our approach can create crystals with a process-tunable thickness, from monolayer to tens of nanometres, and with lateral sizes of up to 100 um. The chiral-chain van der Waals structure of tellurene gives rise to strong in-plane anisotropic properties and large thickness dependent shifts in Raman vibrational modes, which is not observed in other 2D layered materials. We also fabricate tellurene field-effect transistors, which exhibit air-stable performance at room temperature for over two months, on off ratios on the order of 106 and field-effect mobilities of around 700 cm2 per Vs. Furthermore, by scaling down the channel length and integrating with high-k dielectrics, transistors with a significant on-state current density of 1 A mm-1 are demonstrated.

First-principles calculations combining density-functional theory and continuum solvation models enable realistic theoretical modeling and design of electrochemical systems. When a reaction proceeds in such systems, the number of electrons in the portion of the system treated quantum mechanically changes continuously, with a balancing charge appearing in the continuum electrolyte. A grand-canonical ensemble of electrons at a chemical potential set by the electrode potential is therefore the ideal description of such systems that directly mimics the experimental condition. We present two distinct algorithms, a self-consistent field method (GC-SCF) and a direct variational free energy minimization method using auxiliary Hamiltonians (GC-AuxH), to solve the Kohn-Sham equations of electronic density-functional theory directly in the grand canonical ensemble at fixed potential. Both methods substantially improve performance compared to a sequence of conventional fixed-number calculations targeting the desired potential, with the GC-AuxH method additionally exhibiting reliable and smooth exponential convergence of the grand free energy. Finally, we apply grand-canonical DFT to the under-potential deposition of copper on platinum from chloride-containing electrolytes and show that chloride desorption, not partial copper monolayer formation, is responsible for the second voltammetric peak., Comment: 15 pages, 12 figures