skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A shortcut to the thermodynamic limit for quantum many-body calculations of metals
Abstract Computationally efficient and accurate quantum mechanical approximations to solve the many-electron Schrödinger equation are crucial for computational materials science. Methods such as coupled cluster theory show potential for widespread adoption if computational cost bottlenecks can be removed. For example, extremely densek-point grids are required to model long-range electronic correlation effects, particularly for metals. Although these grids can be made more effective by averaging calculations over an offset (or twist angle), the resultant cost in time for coupled cluster theory is prohibitive. We show here that a single special twist angle can be found using the transition structure factor, which provides the same benefit as twist averaging with one or two orders of magnitude reduction in computational time. We demonstrate that this not only works for metal systems but also is applicable to a broader range of materials, including insulators and semiconductors.  more » « less
Award ID(s):
2045046
PAR ID:
10361023
Author(s) / Creator(s):
; ; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Computational Science
Volume:
1
Issue:
12
ISSN:
2662-8457
Page Range / eLocation ID:
p. 801-808
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Moiré lattices formed in twisted van der Waals bilayers provide a unique, tunable platform to realize coupled electron or exciton lattices unavailable before. While twist angle between the bilayer has been shown to be a critical parameter in engineering the moiré potential and enabling novel phenomena in electronic moiré systems, a systematic experimental study as a function of twist angle is still missing. Here we show that not only are moiré excitons robust in bilayers of even large twist angles, but also properties of the moiré excitons are dependant on, and controllable by, the moiré reciprocal lattice period via twist-angle tuning. From the twist-angle dependence, we furthermore obtain the effective mass of the interlayer excitons and the electron inter-layer tunneling strength, which are difficult to measure experimentally otherwise. These findings pave the way for understanding and engineering rich moiré-lattice induced phenomena in angle-twisted semiconductor van der Waals heterostructures. 
    more » « less
  2. Abstract The investigation of twisted stacked few‐layer MoS2has revealed novel electronic, optical, and vibrational properties over an extended period. For the successful integration of twisted stacked few‐layer MoS2into a wide range of applications, it is crucial to employ a noninvasive, versatile technique for characterizing the layered architecture of these complex structures. In this work, we introduce a machine learning‐assisted low‐frequency Raman spectroscopy method to characterize the twist angle of few‐layer stacked MoS2samples. A feedforward neural network (FNN) is utilized to analyze the low‐frequency breathing mode as a function of the twist angle. Moreover, using finite difference method (FDM) and density functional theory (DFT) calculations, we show that the low‐frequency Raman spectra of MoS2are mainly influenced by the effect of the nearest and second nearest layers. A new improved linear chain model (TA‐LCM) with taking the twist angle into the consideration is developed to understand the interlayer breathing modes of stacked few‐layer MoS2. This approach can be extended to other 2D materials systems and provides an intelligent way to investigate naturally stacked and twisted interlayer interactions. 
    more » « less
  3. Abstract Forming a hetero-interface is a materials-design strategy that can access an astronomically large phase space. However, the immense phase space necessitates a high-throughput approach for an optimal interface design. Here we introduce a high-throughput computational framework, InterMatch, for efficiently predicting charge transfer, strain, and superlattice structure of an interface by leveraging the databases of individual bulk materials. Specifically, the algorithm reads in the lattice vectors, density of states, and the stiffness tensors for each material in their isolated form from the Materials Project. From these bulk properties, InterMatch estimates the interfacial properties. We benchmark InterMatch predictions for the charge transfer against experimental measurements and supercell density-functional theory calculations. We then use InterMatch to predict promising interface candidates for doping transition metal dichalcogenide MoSe2. Finally, we explain experimental observation of factor of 10 variation in the supercell periodicity within a few microns in graphene/α-RuCl3by exploring low energy superlattice structures as a function of twist angle using InterMatch. We anticipate our open-source InterMatch algorithm accelerating and guiding ever-growing interfacial design efforts. Moreover, the interface database resulting from the InterMatch searches presented in this paper can be readily accessed online. 
    more » « less
  4. Abstract Reliable simulations of molecules in condensed phase require the combination of an accurate quantum mechanical method for the core region, and a realistic model to describe the interaction with the environment. Additionally, this combination should not significantly increase the computational cost of the calculation compared to the corresponding in vacuo case. In this review, we describe the combination of methods based on coupled cluster (CC) theory with polarizable classical models for the environment. We use the polarizable continuum model (PCM) of solvation to discuss the equations, but we also show how the same theoretical framework can be extended to polarizable force fields. The theory is developed within the perturbation theory energy and singles‐T density (PTES) scheme, where the environmental response is computed with the CC single excitation amplitudes as an approximation for the full one‐particle reduced density. The CC‐PTES combination provides the best compromise between accuracy and computational effort for CC calculations in condensed phase, because it includes the response of the environment to the correlation density at the same computational cost of in vacuo CC. We discuss a number of numerical applications for ground and excited state properties, based on the implementation of CC‐PTES with single and double excitations (CCSD‐PTES), which show the reliability and computational efficiency of the method in reproducing experimental or full‐CC data. This article is characterized under:Electronic Structure Theory > Ab Initio Electronic Structure MethodsElectronic Structure Theory > Combined QM/MM MethodsSoftware > Quantum Chemistry 
    more » « less
  5. Jonathan R. Whitlock (Ed.)
    IntroductionUnderstanding the neural code has been one of the central aims of neuroscience research for decades. Spikes are commonly referred to as the units of information transfer, but multi-unit activity (MUA) recordings are routinely analyzed in aggregate forms such as binned spike counts, peri-stimulus time histograms, firing rates, or population codes. Various forms of averaging also occur in the brain, from the spatial averaging of spikes within dendritic trees to their temporal averaging through synaptic dynamics. However, how these forms of averaging are related to each other or to the spatial and temporal units of information representation within the neural code has remained poorly understood. Materials and methodsIn this work we developed NeuroPixelHD, a symbolic hyperdimensional model of MUA, and used it to decode the spatial location and identity of static images shown ton= 9 mice in the Allen Institute Visual Coding—NeuroPixels dataset from large-scale MUA recordings. We parametrically varied the spatial and temporal resolutions of the MUA data provided to the model, and compared its resulting decoding accuracy. ResultsFor almost all subjects, we found 125ms temporal resolution to maximize decoding accuracy for both the spatial location of Gabor patches (81 classes for patches presented over a 9×9 grid) as well as the identity of natural images (118 classes corresponding to 118 images) across the whole brain. This optimal temporal resolution nevertheless varied greatly between different regions, followed a sensory-associate hierarchy, and was significantly modulated by the central frequency of theta-band oscillations across different regions. Spatially, the optimal resolution was at either of two mesoscale levels for almost all mice: the area level, where the spiking activity of all neurons within each brain area are combined, and the population level, where neuronal spikes within each area are combined across fast spiking (putatively inhibitory) and regular spiking (putatively excitatory) neurons, respectively. We also observed an expected interplay between optimal spatial and temporal resolutions, whereby increasing the amount of averaging across one dimension (space or time) decreases the amount of averaging that is optimal across the other dimension, and vice versa. DiscussionOur findings corroborate existing empirical practices of spatiotemporal binning and averaging in MUA data analysis, and provide a rigorous computational framework for optimizing the level of such aggregations. Our findings can also synthesize these empirical practices with existing knowledge of the various sources of biological averaging in the brain into a new theory of neural information processing in which theunit of informationvaries dynamically based on neuronal signal and noise correlations across space and time. 
    more » « less