Engineering Papers⌕ Search

Engineering topics

Jakowski, Jacek

Publications and source records attributed to Jakowski, Jacek.

Structural, electronic, and optical-absorption properties of 2D Si thin films

Recent experimental studies highlighted the potential of thin-film crystalline silicon (Si) for high-efficiency solar cells. Using density functional theory, we investigated 2D Si thin films across various orientations, thicknesses, and surface structures to elucidate their structure–property relationships. Through surface-energy calculations and Wulff construction, we determined the crystal habit of Si, which aligns with available experimental observations. Electronic-structure calculations underscored the critical role of valence saturation on surfaces in enabling semiconducting behavior in Si thin films, essential for optical applications. From optical-absorption calculations, we identified the surface index exhibiting the highest absorption coefficients for thin films Si solar cell applications.

36 MATERIALS SCIENCE↗

First-Principles Simulation of Beam-Induced Processes Underlying Atomic Manipulation in Electron Microscopes

The development of experimental methods and apparatuses capable of promoting atomically precise material manipulations holds great promise for realizing the ultimate limit of feature miniaturization in materials and devices. The ability to modify materials atom by atom is anticipated to usher in new technologies in areas as diverse as separation science, medicine, and quantum information science. Historically, scanning probe-based techniques have been the most prominent approaches in this space. However, these methods are best suited for the manipulation of surface-exposed regions of materials, as the strong perturbations required for bond scission are delivered most effectively to atoms in the near-proximity to the scanning probe. In contrast, convergent electron beams with energies tuned slightly below the threshold for inducing irreversible knock-on damage have recently been employed (within scanning transmission electron microscopy) to promote atomic-scale bond rearrangements in various beam-stable solids. Currently, however, the efficiency and selectivity of beam-induced atomic manipulation processes with focused electron beams are such that long irradiation times are required to induce a desired atomic rearrangement. With a better understanding of the underlying physics dictating the outcome of a given irradiation event, methods can be devised to improve the efficiency of these techniques so that their promise can be fully realized through widespread adoption.To this end, this Account details our recent efforts to develop and apply tractable first-principles simulation approaches for studying the response of materials to electric beam-like external electric potentials applied in real space. We briefly review the concepts and capabilities in the area of atomically precise materials manipulation and review the early demonstrations of accomplishments in this area, focusing on studies using scanned convergent electron beam probes in particular. We expound upon the depth of the challenge and identify critical shortcomings of theoretical methods that have previously been employed in the simulation of beam-induced processes. We then describe the computational methods that we have generalized from the concepts and tools most commonly applied to the study of molecular photochemistry and how our adaptations of these methods can be employed to capture the relevant dynamical phenomena for beam-induced processes ranging from the initial electron scattering to the ensuing multistate reactions. Here, we contextualize these methods within the current state of the art in this area, which has historically focused primarily on the simulation of inelastic image formation in the electron microscope for the purpose of interpreting the results of quantitative electron microscopy experiments. We demonstrate that the spatial distribution of state-specific excitation rates due to the presence of an external (probe) electric charge is inhomogeneous, such that irradiation at particular locations in materials can favor specific electronic transitions (and disallow others). In addition to the potential for excited-state reaction pathways to be accessed through the initial inelastic scattering of the tightly focused electron beam from the targeted atoms, we also identify favorable conditions for the electronically nonadiabatic evolution of the highly vibrationally excited system to open complex multistate reaction pathways. Implications of the early results for understanding the mechanisms and potential routes to improved efficiency and selectivity in beam-induced reactions are discussed. We conclude with a summary of the current state of theory and modeling capabilities in this area and provide our perspective on future directions for theoretical and experimental developments that we view as crucial to advancing the use of convergent electron beams in mode-specific, atomically precise platforms for direct-write materials modifications.

36 MATERIALS SCIENCE↗

Theoretical Examination of the Hydroxide Transport in Cobaltocenium-Containing Polyelectrolytes

Polymers incorporating cobaltocenium groups have received attention as promising components of anion-exchange membranes (AEMs), exhibiting a good balance of chemical stability and high ionic conductivity. In this work, we analyze the hydroxide diffusion in the presence of cobaltocenium cations in an aqueous environment based on the molecular dynamics of model systems confined in one dimension to mimic the AEM channels. In order to describe the proton hopping mechanism, the forces are obtained from the electronic structure computed at the density-functional tight-binding level. We find that the hydroxide diffusion depends on the channel size, modulation of the electrostatic interactions by the solvation shell, and its rearrangement ability. Hydroxide diffusion proceeds via both the vehicular and structural diffusion mechanisms with the latter playing a larger role at low diffusion coefficients. The highest diffusion coefficient is observed under moderate water densities (around half the density of liquid water) when there are enough water molecules to form the solvation shell, reducing the electrostatic interaction between ions, yet there is enough space for the water rearrangements during the proton hopping. Furthermore, the effects of cobaltocenium separation, orientation, chemical modifications, and the role of nuclear quantum effects are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gene Expression Programming for Quantum Computing

Here, we introduce QuantumGEP, a scientific computer program that uses gene expression programming (GEP) to find a quantum circuit that either (1) maps a given set of input states to a given set of output states or (2) transforms a fixed initial state to minimize a given physical quantity of the output state. QuantumGEP is a driver program that uses evendim, a generic computational engine for GEP, both of which are free and open source. We apply QuantumGEP as a powerful solver for MaxCut in graphs and for condensed matter quantum many-body Hamiltonians.

97 MATHEMATICS AND COMPUTING↗

Quantum Chemical Simulations of CO 2 and N 2 Capture in Reline, a Prototypical Deep Eutectic Solvent

Deep eutectic solvents such as reline are an emerging class of low-cost, environmentally friendly solvents with tunable properties that are potentially applicable for the capture and separation of CO 2 . Experimental measurements showed that a reline-based membrane contactor can capture and separate CO 2 via physisorption through a dissolution process with 96.7% purity from a mixed gas containing CO 2 and N 2 (50:50% molar ratio). Here, we examine the nature of the interaction of CO 2 and N 2 with reline employing quantum chemical methods. We focus on explaining the mechanism by which CO 2 and N 2 bind to reline and the reason for the high selectivity for absorption of CO 2 compared to N 2 . We analyze the dynamics, energetics, and binding motifs for CO 2 and N 2 in reline employing density functional theory, density functional tight binding, and ab initio molecular dynamics. We also investigate the effect of reline on the vibrational spectra of CO 2 and reline. Our simulations indicate that the selective capture of CO 2 from the mixture of CO 2 and N 2 is due to the interplay between attractive electrostatic and charge polarization forces with opposing entropic effects, which shift the energetic balance and make the N 2 absorption unfavorable in reline.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Symmetry in VQE: A Dual Approach for Trapped-Ion Simulations of Benzene

Understanding complex chemical systems—such as biomolecules, catalysts, and novel materials—is a central goal of quantum simulations. Near-term strategies hinge on the use of variational quantum eigensolver (VQE) algorithms combined with a suitable ansatz. However, straightforward application of many chemically-inspired ansatze yields prohibitively deep circuits. In this work, we employ several circuit optimization methods tailored for trapped-ion quantum devices to enhance the feasibility of intricate chemical simulations. The techniques aim to lessen the depth of the unitary coupled cluster with singles and doubles (uCCSD) ansatz's circuit compilation, a considerable challenge on current noisy quantum devices. Furthermore, we use symmetry-inspired classical post-selection methods to further refine the outcomes and minimize errors in energy measurements, without adding quantum overhead. Our strategies encompass optimal mapping from orbital to qubit, term reordering to minimize entangling gates, and the exploitation of molecular spin and point group symmetry to eliminate redundant parameters. The inclusion of error mitigation via post-selection based on known molecular symmetries improves the results to near milli-Hartree accuracy. These methods, when applied to a benzene molecule simulation, enabled the construction of an 8-qubit circuit with 69 two-qubit entangling operations, pushing the limits for variational quantum eigensolver (VQE) circuits executed on quantum hardware to date. 1 1 This manuscript has been authored in part by UT-Battelle, LLC, under contract DE-AC05-000R22725 with the US Department of Energy (DOE). The publisher acknowledges the US government license to provide public access under the DOE Public Access Plan(https://energy.gov/doe-public-access-plan).

Goings, Joshua↗

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Membrane Contactor Enabling Energy-Efficient CO 2 Capture from Point Sources with Deep Eutectic Solvents

We demonstrate a scalable and energy-efficient hollow fiber membrane contactor (HFMC)-based process using a green solvent for CO 2 capture. This process uses a deep eutectic solvent (DES) in an HFMC to provide close interfacial interactions and contact between the DES and CO 2 . This approach overcomes disadvantages associated with direct absorption in DES and could potentially be applied to a variety of solvent-based CO 2 capture methods. Commercial low-cost polymer hollow fiber membranes (e.g., microporous polypropylene) were evaluated for CO 2 capture with reline, a prototypical DES. Single-gas measurements showed that the DES-based polypropylene HFMC can capture and separate CO 2 while rejecting N2. From a mixed gas containing 50 mol % N 2 and 50 mol % CO ν , the DES-based HFMC separated CO 2 with a purity of 96.9 mol %. The effect of several process parameters including solvent flow rate, pressure, and temperature on the CO 2 separation performance was studied. The flux of the recovered CO 2 was 67.43 mmole/m 2 /h at a feed pressure of 4 bar. In situ Fourier transform infrared (FTIR) measurements combined with density functional theory (DFT)-based molecular dynamics simulations revealed that reline absorbs CO 2 by physical absorption without forming a new chemical compound, and CO 2 separation by reline occurs via the pressure swing mechanism. This research provides fundamental insights about physical solvent-based separation processes and a pathway toward practical deployment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Extracting Inelastic Scattering Cross Sections for Finite and Aperiodic Materials from Electronic Dynamics Simulations

Explicit time-dependent electronic structure theory methods are increasingly prevalent in the areas of condensed matter physics and quantum chemistry, with the broad-band optical absorptivity of molecular and small condensed-phase systems nowadays routinely studied with such approaches. Here, in this paper, it is demonstrated that electronic dynamics simulations can similarly be employed to study cross sections for the scattering-induced electronic excitations probed in nonresonant inelastic X-ray scattering and momentum-resolved electron energy loss spectroscopies. A method is put forth for evaluating the electronic dynamic structure factor, which involves the application of a momentum boost-type perturbation and transformation of the resulting reciprocal space density fluctuations into the frequency domain. Good agreement is first demonstrated between the dynamic structure factor extracted from these electronic dynamics simulations and the corresponding transition matrix elements from linear response theory. The method is then applied to some extended (quasi)one-dimensional systems, for which the wave vector becomes a good quantum number in the thermodynamic limit. Finally, the dispersion of many-body excitations in a series of hydrogen-terminated graphene flakes (and twisted bilayers thereof) is investigated to highlight the utility of the presented approach for capturing morphology-dependent effects in the inelastic scattering cross sections of nanostructured and/or noncrystalline materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗