Matching Items (31)
Filtering by
Description
Monte Carlo methods often used in nuclear physics, such as auxiliary field diffusion Monte Carlo and Green's function Monte Carlo, have typically relied on phenomenological local real-space potentials containing as few derivatives as possible, such as the Argonne-Urbana family of interactions, to make sampling simple and efficient. Basis set methods such as no-core shell model or coupled-cluster techniques typically use softer non-local potentials because of their more rapid convergence with basis set size. These non-local potentials are typically defined in momentum space and are often based on effective field theory. Comparisons of the results of the two types of methods are complicated by the use of different potentials. This thesis discusses progress made in using such non-local potentials in quantum Monte Carlo calculations of light nuclei. In particular, it shows methods for evaluating the real-space, imaginary-time propagators needed to perform quantum Monte Carlo calculations using non-local potentials and universality properties of these propagators, how to formulate a good trial wave function for non-local potentials, and how to perform a "one-step" Green's function Monte Carlo calculation for non-local potentials.
ContributorsLynn, Joel E (Author) / Schmidt, Kevin E (Thesis advisor) / Alarcon, Ricardo (Committee member) / Lebed, Richard (Committee member) / Shovkovy, Igor (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2013
Description
The theory of quantum electrodynamics predicts that beta decay of the neutron into a proton, electron, and anti-neutrino should be accompanied by a continuous spectrum of photons. A recent experiment, RDK I, reported the first detection of radiative decay photons from neutron beta decay with a branching ratio of (3.09 ± 0.32) × 10-3 in the energy range of 15 keV to 340 keV. This was achieved by prompt coincident detection of an electron and photon, in delayed coincidence with a proton. The photons were detected by using a single bar of bismuth germanate scintillating crystal coupled to an avalanche photodiode. This thesis deals with the follow-up experiment, RDK II, to measure the branching ratio at the level of approximately 1% and the energy spectrum at the level of a few percent. The most significant improvement of RDK II is the use of a photon detector with about an order of magnitude greater solid angle coverage than RDK I. In addition, the detectable energy range has been extended down to approximately 250 eV and up to the endpoint energy of 782 keV. This dissertation presents an overview of the apparatus, development of a new data analysis technique for radiative decay, and results for the ratio of electron-proton-photon coincident Repg to electron-proton coincident Rep events.
ContributorsO'Neill, Benjamin (Author) / Alarcon, Ricardo (Thesis advisor) / Drucker, Jeffery (Committee member) / Lebed, Richard (Committee member) / Comfort, Joseph (Committee member) / Chamberlin, Ralph (Committee member) / Arizona State University (Publisher)
Created2012
Description
Spin-orbit interactions are important in determining nuclear structure. They lead to a shift in the energy levels in the nuclear shell model, which could explain the sequence of magic numbers in nuclei. Also in nucleon-nucleon scattering, the large nucleon polarization observed perpendicular to the plane of scattering needs to be explained by adding the spin-orbit interactions in the potential. Their effects change the equation of state and other properties of nuclear matter. Therefore, the simulation of spin-orbit interactions is necessary in nuclear matter.
The auxiliary field diffusion Monte Carlo is an effective and accurate method for calculating the ground state and low-lying exited states in nuclei and nuclear matter. It has successfully employed the Argonne v6' two-body potential to calculate the equation of state in nuclear matter, and has been applied to light nuclei with reasonable agreement with experimental results. However, the spin-orbit interactions were not included in the previous simulations, because the isospin-dependent spin-orbit potential is difficult in the quantum Monte Carlo method. This work develops a new method using extra auxiliary fields to break up the interactions between nucleons, so that the spin-orbit interaction with isospin can be included in the Hamiltonian, and ground-state energy and other properties can be obtained.
The auxiliary field diffusion Monte Carlo is an effective and accurate method for calculating the ground state and low-lying exited states in nuclei and nuclear matter. It has successfully employed the Argonne v6' two-body potential to calculate the equation of state in nuclear matter, and has been applied to light nuclei with reasonable agreement with experimental results. However, the spin-orbit interactions were not included in the previous simulations, because the isospin-dependent spin-orbit potential is difficult in the quantum Monte Carlo method. This work develops a new method using extra auxiliary fields to break up the interactions between nucleons, so that the spin-orbit interaction with isospin can be included in the Hamiltonian, and ground-state energy and other properties can be obtained.
ContributorsZhang, Jie (Author) / Schmidt, Kevin E (Thesis advisor) / Alarcon, Ricardo (Committee member) / Lebed, Richard (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2014
Description
The properties of materials depend heavily on the spatial distribution and connectivity of their constituent parts. This applies equally to materials such as diamond and glasses as it does to biomolecules that are the product of billions of years of evolution. In science, insight is often gained through simple models with characteristics that are the result of the few features that have purposely been retained. Common to all research within in this thesis is the use of network-based models to describe the properties of materials. This work begins with the description of a technique for decoupling boundary effects from intrinsic properties of nanomaterials that maps the atomic distribution of nanomaterials of diverse shape and size but common atomic geometry onto a universal curve. This is followed by an investigation of correlated density fluctuations in the large length scale limit in amorphous materials through the analysis of large continuous random network models. The difficulty of estimating this limit from finite models is overcome by the development of a technique that uses the variance in the number of atoms in finite subregions to perform the extrapolation to large length scales. The technique is applied to models of amorphous silicon and vitreous silica and compared with results from recent experiments. The latter part this work applies network-based models to biological systems. The first application models force-induced protein unfolding as crack propagation on a constraint network consisting of interactions such as hydrogen bonds that cross-link and stabilize a folded polypeptide chain. Unfolding pathways generated by the model are compared with molecular dynamics simulation and experiment for a diverse set of proteins, demonstrating that the model is able to capture not only native state behavior but also partially unfolded intermediates far from the native state. This study concludes with the extension of the latter model in the development of an efficient algorithm for predicting protein structure through the flexible fitting of atomic models to low-resolution cryo-electron microscopy data. By optimizing the fit to synthetic data through directed sampling and context-dependent constraint removal, predictions are made with accuracies within the expected variability of the native state.
ContributorsDe Graff, Adam (Author) / Thorpe, Michael F. (Thesis advisor) / Ghirlanda, Giovanna (Committee member) / Matyushov, Dmitry (Committee member) / Ozkan, Sefika B. (Committee member) / Treacy, Michael M. J. (Committee member) / Arizona State University (Publisher)
Created2011
Description
Quark matter at sufficiently high density and low temperature is expected to be a color superconductor, and may exist in the interior of neutron stars. The properties of two simplest possible color-superconducting phases, i.e., the color-flavor-locked (CFL) and two-flavor superconducting (2SC) phases, are reviewed. The effect of a magnetic field on the pairing dynamics in two-flavor color-superconducting dense quark matter is investigated. A universal form of the gap equation for an arbitrary magnetic field is derived in the weakly coupled regime of QCD at asymptotically high density, using the framework of Schwinger-Dyson equation in the improved rainbow approximation. The results for the gap in two limiting cases, weak and strong magnetic fields, are obtained and discussed. It is shown that the superconducting gap function in the weak magnetic field limit develops a directional dependence in momentum space. This property of the gap parameter is argued to be a consequence of a long-range interaction in QCD.
ContributorsYu, Lang (Author) / Shovkovy, Igor A. (Thesis advisor) / Lunardini, Cecilia (Committee member) / Schmidt, Kevin (Committee member) / Alarcon, Ricardo (Committee member) / Lebed, Richard (Committee member) / Arizona State University (Publisher)
Created2012
Description
Conformational changes in biomolecules often take place on longer timescales than are easily accessible with unbiased molecular dynamics simulations, necessitating the use of enhanced sampling techniques, such as adaptive umbrella sampling. In this technique, the conformational free energy is calculated in terms of a designated set of reaction coordinates. At the same time, estimates of this free energy are subtracted from the potential energy in order to remove free energy barriers and cause conformational changes to take place more rapidly. This dissertation presents applications of adaptive umbrella sampling to a variety of biomolecular systems. The first study investigated the effects of glycosylation in GalNAc2-MM1, an analog of glycosylated macrophage activating factor. It was found that glycosylation destabilizes the protein by increasing the solvent exposure of hydrophobic residues. The second study examined the role of bound calcium ions in promoting the isomerization of a cis peptide bond in the collagen-binding domain of Clostridium histolyticum collagenase. This study determined that the bound calcium ions reduced the barrier to the isomerization of this peptide bond as well as stabilizing the cis conformation thermodynamically, and identified some of the reasons for this. The third study represents the application of GAMUS (Gaussian mixture adaptive umbrella sampling) to on the conformational dynamics of the fluorescent dye Cy3 attached to the 5' end of DNA, and made predictions concerning the affinity of Cy3 for different base pairs, which were subsequently verified experimentally. Finally, the adaptive umbrella sampling method is extended to make use of the roll angle between adjacent base pairs as a reaction coordinate in order to examine the bending both of free DNA and of DNA bound to the archaeal protein Sac7d. It is found that when DNA bends significantly, cations from the surrounding solution congregate on the concave side, which increases the flexibility of the DNA by screening the repulsion between phosphate backbones. The flexibility of DNA on short length scales is compared to the worm-like chain model, and the contribution of cooperativity in DNA bending to protein-DNA binding is assessed.
ContributorsSpiriti, Justin Matthew (Author) / van der Vaart, Arjan (Thesis advisor) / Chizmeshya, Andrew (Thesis advisor) / Matyushov, Dmitry (Committee member) / Fromme, Petra (Committee member) / Arizona State University (Publisher)
Created2011
Description
Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into the free energy of the system, the dynamical stability of contacts and the role of correlated motions in directing the motions of the system. In this thesis, I use molecular dynamics simulations to provide molecular mechanisms that rationalize structural, thermodynamic, and mutation data on the interactions between the lac repressor headpiece and its O1 operator DNA as well as the ERK2 protein kinase. I performed molecular dynamics simulations of the lac repressor headpiece - O1 operator complex at the natural angle as well as at under- and overbent angles to assess the factors that determine the natural DNA bending angle. I find both energetic and entropic factors contribute to recognition of the natural angle. At the natural angle the energy of the system is minimized by optimization of protein-DNA contacts and the entropy of the system is maximized by release of water from the protein-DNA interface and decorrelation of protein motions. To identify the mechanism by which mutations lead to auto-activation of ERK2, I performed a series of molecular dynamics simulations of ERK1/2 in various stages of activation as well as the constitutively active Q103A, I84A, L73P and R65S ERK2 mutants. My simulations indicate the importance of domain closure for auto-activation and activity regulation. My results enable me to predict two loss-of-function mutants of ERK2, G83A and Q64C, that have been confirmed in experiments by collaborators. One of the powerful capabilities of MD simulations in biochemistry is the ability to find low free energy pathways that connect and explain disparate structural data on biomolecular systems. An extention of the targeted molecular dynamics technique using constraints on internal coordinates will be presented and evaluated. The method gives good results for the alanine dipeptide, but breaks down when applied to study conformational changes in GroEL and adenylate kinase.
ContributorsBarr, Daniel Alan (Author) / van der Vaart, Arjan (Thesis advisor) / Matyushov, Dmitry (Committee member) / Wolf, George (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2011
Description
In a typical living cell, millions to billions of proteins—nanomachines that fluctuate and cycle among many conformational states—convert available free energy into mechanochemical work. A fundamental goal of biophysics is to ascertain how 3D protein structures encode specific functions, such as catalyzing chemical reactions or transporting nutrients into a cell. Protein dynamics span femtosecond timescales (i.e., covalent bond oscillations) to large conformational transition timescales in, and beyond, the millisecond regime (e.g., glucose transport across a phospholipid bilayer). Actual transition events are fast but rare, occurring orders of magnitude faster than typical metastable equilibrium waiting times. Equilibrium molecular dynamics (EqMD) can capture atomistic detail and solute-solvent interactions, but even microseconds of sampling attainable nowadays still falls orders of magnitude short of transition timescales, especially for large systems, rendering observations of such "rare events" difficult or effectively impossible.
Advanced path-sampling methods exploit reduced physical models or biasing to produce plausible transitions while balancing accuracy and efficiency, but quantifying their accuracy relative to other numerical and experimental data has been challenging. Indeed, new horizons in elucidating protein function necessitate that present methodologies be revised to more seamlessly and quantitatively integrate a spectrum of methods, both numerical and experimental. In this dissertation, experimental and computational methods are put into perspective using the enzyme adenylate kinase (AdK) as an illustrative example. We introduce Path Similarity Analysis (PSA)—an integrative computational framework developed to quantify transition path similarity. PSA not only reliably distinguished AdK transitions by the originating method, but also traced pathway differences between two methods back to charge-charge interactions (neglected by the stereochemical model, but not the all-atom force field) in several conserved salt bridges. Cryo-electron microscopy maps of the transporter Bor1p are directly incorporated into EqMD simulations using MD flexible fitting to produce viable structural models and infer a plausible transport mechanism. Conforming to the theme of integration, a short compendium of an exploratory project—developing a hybrid atomistic-continuum method—is presented, including initial results and a novel fluctuating hydrodynamics model and corresponding numerical code.
Advanced path-sampling methods exploit reduced physical models or biasing to produce plausible transitions while balancing accuracy and efficiency, but quantifying their accuracy relative to other numerical and experimental data has been challenging. Indeed, new horizons in elucidating protein function necessitate that present methodologies be revised to more seamlessly and quantitatively integrate a spectrum of methods, both numerical and experimental. In this dissertation, experimental and computational methods are put into perspective using the enzyme adenylate kinase (AdK) as an illustrative example. We introduce Path Similarity Analysis (PSA)—an integrative computational framework developed to quantify transition path similarity. PSA not only reliably distinguished AdK transitions by the originating method, but also traced pathway differences between two methods back to charge-charge interactions (neglected by the stereochemical model, but not the all-atom force field) in several conserved salt bridges. Cryo-electron microscopy maps of the transporter Bor1p are directly incorporated into EqMD simulations using MD flexible fitting to produce viable structural models and infer a plausible transport mechanism. Conforming to the theme of integration, a short compendium of an exploratory project—developing a hybrid atomistic-continuum method—is presented, including initial results and a novel fluctuating hydrodynamics model and corresponding numerical code.
ContributorsSeyler, Sean L (Author) / Beckstein, Oliver (Thesis advisor) / Chamberlin, Ralph (Committee member) / Matyushov, Dmitry (Committee member) / Thorpe, Michael F (Committee member) / Vaiana, Sara (Committee member) / Arizona State University (Publisher)
Created2017
Description
A series of experiments using a polarized beam incident on a polarized frozen spin target
(FROST) was conducted at Jefferson Lab in 2010. Results presented here were taken
during the second running period with the FROST target using the CEBAF Large Acceptance
Spectrometer (CLAS) detector at Jefferson Lab, which used transversely-polarized
protons in a butanol target and a circularly-polarized incident tagged photon beam with
energies between 0.62 and 2.93 GeV. Data are presented for the F and T polarization observables
for h meson photoproduction on the proton from W = 1.55 GeV to 1.80 GeV.
The data presented here will improve the world database and refine theoretical approaches
of nucleon structure.
(FROST) was conducted at Jefferson Lab in 2010. Results presented here were taken
during the second running period with the FROST target using the CEBAF Large Acceptance
Spectrometer (CLAS) detector at Jefferson Lab, which used transversely-polarized
protons in a butanol target and a circularly-polarized incident tagged photon beam with
energies between 0.62 and 2.93 GeV. Data are presented for the F and T polarization observables
for h meson photoproduction on the proton from W = 1.55 GeV to 1.80 GeV.
The data presented here will improve the world database and refine theoretical approaches
of nucleon structure.
ContributorsTucker, Ross (Author) / Ritchie, Barry (Thesis advisor) / Dugger, Michael (Committee member) / Alarcon, Ricardo (Committee member) / Lebed, Richard (Committee member) / Arizona State University (Publisher)
Created2016
Description
In this thesis, I present the study of nucleon structure from distinct perspectives. I start by elaborating the motivations behind the endeavors and then introducing the key concept, namely the generalized parton distribution functions (GPDs), which serves as the frame- work describing hadronic particles in terms of their fundamental constituents. The second chapter is then devoted to a detailed phenomenological study of the Virtual Compton Scattering (VCS) process, where a more comprehensive parametrization is suggested. In the third chapter, the renormalization kernels that enters the QCD evolution equations at twist- 4 accuracy are computed in terms of Feynman diagrams in momentum space, which can be viewed as an extension of the work by Bukhvostov, Frolov, Lipatov, and Kuraev (BKLK). The results can be used for determining the QCD background interaction for future precision measurements.
ContributorsJi, Yao, Ph. D (Author) / Belitsky, Andrei (Thesis advisor) / Lebed, Richard (Committee member) / Schmidt, Kevin E (Committee member) / Vachaspati, Tanmay (Committee member) / Arizona State University (Publisher)
Created2016