Temporal difference learning with linear function approximation is a popular method to obtain a low-dimensional approximation of the value function of a policy in a Markov Decision Process. We give a new interpretation of this method in terms of a splitting of the gradient of an appropriately chosen function. As a consequence of this interpretation, convergence proofs for gradient descent can be applied almost verbatim to temporal difference learning. Beyond giving a new, fuller explanation of why temporal difference works, our interpretation also yields improved convergence times. We consider the setting with 1/T^{1/2} step-size, where previous comparable finite-time convergence time bounds for temporal difference learning had the multiplicative factor 1/(1-\gamma) in front of the bound, with γ being the discount factor. We show that a minor variation on TD learning which estimates the mean of the value function separately has a convergence time where 1/(1-\gamma) only multiplies an asymptotically negligible term.
more »
« less
Machine learning for a finite size correction in periodic coupled cluster theory calculations
We introduce a straightforward Gaussian process regression (GPR) model for the transition structure factor of metal periodic coupled cluster singles and doubles (CCSD) calculations. This is inspired by the method introduced by Liao and Gr\"uneis for interpolating over the transition structure factor to obtain a finite size correction for CCSD [J. Chem. Phys. 145, 141102 (2016)], and by our own prior work using the transition structure factor to efficiently converge CCSD for metals to the thermodynamic limit [Nat. Comput. Sci. 1, 801 (2021)]. In our CCSD-FS-GPR method to correct for finite size errors, we fit the structure factor to a 1D function in the momentum transfer, $$G$$.We then integrate over this function by projecting it onto a k-point mesh to obtain comparisons with extrapolated results. Results are shown for lithium, sodium, and the uniform electron gas.
more »
« less
- Award ID(s):
- 2045046
- PAR ID:
- 10326873
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Hydrazoic acid (HN3) is used as a case study for investigating the accuracy and precision by which a molecular structure—specifically, a semi-experimental equilibrium structure (reSE)—may be determined using current state-of-the-art methodology. The influence of the theoretical corrections for effects of vibration–rotation coupling and electron-mass distribution that are employed in the analysis is explored in detail. The small size of HN3 allowed us to deploy considerable computational resources to probe the basis-set dependence of these corrections using a series of coupled-cluster single, double, perturbative triple [CCSD(T)] calculations with cc-pCVXZ (X = D, T, Q, 5) basis sets. We extrapolated the resulting corrections to the complete basis set (CBS) limit to obtain CCSD(T)/CBS corrections, which were used in a subsequent reSE structure determination. The reSE parameters obtained using the CCSD(T)/cc-pCV5Z corrections are nearly identical to those obtained using the CCSD(T)/CBS corrections, with uncertainties in the bond distances and angles of less than 0.0006 Å and 0.08°, respectively. The previously obtained reSE structure using CCSD(T)/ANO2 agrees with that using CCSD(T)/cc-pCV5Z to within 0.000 08 Å and 0.016° for bond distances and angles, respectively, and with only 25% larger uncertainties, validating the idea that reSE structure determinations can be carried out with significantly smaller basis sets than those needed for similarly accurate, strictly ab initio determinations. Although the purely computational re structural parameters [CCSD(T)/cc-pCV6Z] fall outside of the statistical uncertainties (2σ) of the corresponding reSE structural parameters, the discrepancy is rectified by applying corrections to address the theoretical limitations of the CCSD(T)/cc-pCV6Z geometry with respect to basis set, electron correlation, relativity, and the Born–Oppenheimer approximation, thereby supporting the contention that the semi-experimental approach is both an accurate and vastly more efficient method for structure determinations than is brute-force computation.more » « less
-
Geosynthetic reinforced soil structures are widely used for earth retention and stabilization in many geotechnical and transportation applications. In the traditional design of geosynthetic reinforced soil structures, factor of safety is used to address the uncertainties. However, this approach cannot systematically consider the uncertainties and usually result in over-conservativeness and inconsistence in the design practice. In this paper, a reliability assessment framework of geosynthetic reinforced soil structure design is developed using probabilistic and numerical methods. The geosynthetic reinforced soil structures are modeled using finite element method. In the finite element method, the soil behavior is modelled using the Mohr-Coulomb soil model and a strength reduction method is used to determine the factor of safety value for a given geosynthetic reinforced soil structure. Then the reliability method is combined with the finite element method to obtain the probability curves for the geosynthetic reinforced soil structure. A case study of geosynthetic reinforced soil wall design is used to illustrate the significance of the proposed framework. The reliability assessment framework provides a useful decision making tool for informed design of geosynthetic reinforced soil structures in the face of uncertainties.more » « less
-
In density-functional theory, the exchange–correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n̄xc(r,r′), representing the probability depletion of finding an electron at r′ due to an electron at r. Accurate knowledge of n̄xc(r,r′) has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron–electron cusp using finite basis sets. Focusing on Hooke’s atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke’s atoms, the error introduced by the description of the electron–electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew–Burke–Ernzerhof and local density approximation XC hole models are assessed.more » « less
-
Sloan Siegrist, M. (Ed.)ABSTRACT The Agrobacterium growth pole ring (GPR) protein forms a hexameric ring at the growth pole (GP) that is essential for polar growth. GPR is large (2,115 amino acids) and contains 1,700 amino acids of continuous α-helices. To dissect potential GPR functional domains, we created deletions of regions with similarity to human apolipoprotein A-IV (396 amino acids), itself composed of α-helical domains. We also tested deletions of the GPR C terminus. Deletions were inducibly expressed as green fluorescent protein (GFP) fusion proteins and tested for merodiploid interference with wild-type (WT) GPR function, for partial function in cells lacking GPR, and for formation of paired fluorescent foci (indicative of hexameric rings) at the GP. Deletion of domains similar to human apolipoprotein A-IV in GPR caused defects in cell morphology when expressed in trans to WT GPR and provided only partial complementation to cells lacking GPR. Agrobacterium -specific domains A-IV-1 and A-IV-4 contain predicted coiled coil (CC) regions of 21 amino acids; deletion of CC regions produced severe defects in cell morphology in the interference assay. Mutants that produced the most severe effects on cell shape also failed to form paired polar foci. Modeling of A-IV-1 and A-IV-4 reveals significant similarity to the solved structure of human apolipoprotein A-IV. GPR C-terminal deletions profoundly blocked complementation. Finally, peptidoglycan (PG) synthesis is abnormally localized circumferentially in cells lacking GPR. The results support the hypothesis that GPR plays essential roles as an organizing center for membrane and PG synthesis during polar growth. IMPORTANCE Bacterial growth and division are extensively studied in model systems ( Escherichia coli , Bacillus subtilis , and Caulobacter crescentus ) that grow by dispersed insertion of new cell wall material along the length of the cell. An alternative growth mode—polar growth—is used by some Actinomycetales and Proteobacteria species. The latter phylum includes the family Rhizobiaceae , in which many species, including Agrobacterium tumefaciens , exhibit polar growth. Current research aims to identify growth pole (GP) factors. The Agrobacterium growth pole ring (GPR) protein is essential for polar growth and forms a striking hexameric ring structure at the GP. GPR is long (2,115 amino acids), and little is known about regions essential for structure or function. Genetic analyses demonstrate that the C terminus of GPR, and two internal regions with homology to human apolipoproteins (that sequester lipids), are essential for GPR function and localization to the GP. We hypothesize that GPR is an organizing center for membrane and cell wall synthesis during polar growth.more » « less
An official website of the United States government

