Physics-Informed Neural Networks for Cardiac Activation Mapping
Not Available
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SEARCH · Engineering Papers
Results for “Mathematical Physics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
PUMA-V - Polyhedral User Mapping Assistant and Visualizer Final Technical Report
Explore the source record for details and available documents.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)
Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.
97 MATHEMATICS AND COMPUTING↗
What can solve the strong CP problem?
Three possible strategies have been advocated to solve the strong CP problem. The first is the axion, a dynamical mechanism that relaxes any initial value of the CP violating angle $\overline{θ}$ to zero. The second is the imposition of new symmetries that are believed to set $\overline{θ}$ to zero in the UV. The third is the acceptance of the fine tuning of parameters. We argue that the latter two solutions do not solve the strong CP problem. The θ term of QCD is not a parameter — it does not exist in the Hamiltonian. Rather, it is a property of the quantum state that our universe finds itself in, arising from the fact that there are CP violating states of a CP preserving Hamiltonian. It is not eliminated by imposing parity as a symmetry since the underlying theory is already parity symmetric and that does not preclude the existence of CP violating states. Moreover, since the value of θ realized in our universe is a consequence of measurement, it is inherently random and cannot be fine tuned by choice of parameters. Rather any fine tuning would require a tuning between parameters in the theory and the random outcome of measurement. Our results considerably strengthen the case for the existence of the axion and axion dark matter. The confusion around θ arises from the fact that unlike classical mechanics, the Hamiltonian and Lagrangian are not equivalent in quantum mechanics. The Hamiltonian defines the differential time evolution, whereas the Lagrangian is a solution to this evolution. Consequently, initial conditions could in principle appear in the Lagrangian but not in the Hamiltonian. This results in aspects of the initial condition such as θ misleadingly appearing in the Lagrangian as parameters. We comment on the similarity between the θ vacua and the violations of the constraint equations of classical gauge theories in quantum mechanics.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
CICART Center For Integrated Computation And Analysis Of Reconnection And Turbulence
CICART is a partnership between the University of New Hampshire (UNH) and Dartmouth College. CICART addresses two important science needs of the DoE: the basic understanding of magnetic reconnection and turbulence that strongly impacts the performance of fusion plasmas, and the development of new mathematical and computational tools that enable the modeling and control of these phenomena. The principal participants of CICART constitute an interdisciplinary group, drawn from the communities of applied mathematics, astrophysics, computational physics, fluid dynamics, and fusion physics. It is a main premise of CICART that fundamental aspects of magnetic reconnection and turbulence in fusion devices, smaller-scale laboratory experiments, and space and astrophysical plasmas can be viewed from a common perspective, and that progress in understanding in any of these interconnected fields is likely to lead to progress in others. The establishment of CICART has strongly impacted the education and research mission of a new Program in Integrated Applied Mathematics in the College of Engineering and Applied Sciences at UNH by enabling the recruitment of a tenure-track faculty member, supported equally by UNH and CICART, and the establishment of an IBM-UNH Computing Alliance. The proposed areas of research in magnetic reconnection and turbulence in astrophysical, space, and laboratory plasmas include the following topics: (A) Reconnection and secondary instabilities in large high-Lundquist-number plasmas, (B) Particle acceleration in the presence of multiple magnetic islands, (C) Gyrokinetic reconnection: comparison with fluid and particle-in-cell models, (D) Imbalanced turbulence, (E) Ion heating, and (F) Turbulence in laboratory (including fusion-relevant) experiments. These theoretical studies make active use of three high-performance computer simulation codes: (1) The Magnetic Reconnection Code, based on extended two-fluid (or Hall MHD) equations, in an Adaptive Mesh Refinement (AMR) framework, (2) the Particle Simulation Code, a fully electromagnetic 3D Particle-In-Cell (PIC) code that includes a collision operator, and (3) GS2, an Eulerian, electromagnetic, kinetic code that is widely used in the fusion program, and simulates the nonlinear gyrokinetic equations, together with a self-consistent set of Maxwell’s equations.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
A Measurement of Atmospheric Circular Polarization with POLARBEAR
At millimeter wavelengths, the atmospheric emission is circularly polarized owing to the Zeeman splitting of molecular oxygen by the Earth's magnetic field. We report a measurement of the signal in the 150 GHz band using 3 yr of observational data with the POLARBEAR project. Nonidealities of a continuously rotating half-wave plate (HWP) partially convert circularly polarized light to linearly polarized light. While POLARBEAR detectors are sensitive to linear polarization, this effect makes them sensitive to circular polarization. Although this was not the intended use, we utilized this conversion to measure circular polarization. We reconstruct the azimuthal gradient of the circular polarization signal and measure its dependency from the scanning direction and the detector bandpass. We compare the signal with a simulation based on atmospheric emission theory, the detector bandpass, and the HWP leakage spectrum model. We find the ratio of the observed azimuthal slope to the simulated slope is 0.92 ± 0.01(stat) ± 0.07(sys). This ratio corresponds to a brightness temperature of 3.8 mK at the effective band center of 121.8 GHz and bandwidth of 3.5 GHz estimated from representative detector bandpass and the spectrum of Zeeman emission. This result validates our understanding of the instrument and reinforces the feasibility of measuring the circular polarization using the imperfection of continuously rotating HWP. Continuously rotating HWP is popular in ongoing and future cosmic microwave background experiments to modulate the polarized signal. This work shows a method for signal extraction and leakage subtraction that can help measure circular polarization in such experiments.
79 ASTRONOMY AND ASTROPHYSICS↗
Radiation dose reduction in medical x-ray CT via Fourier-based iterative reconstruction
Purpose: A Fourier-based iterative reconstruction technique, termed Equally Sloped Tomography (EST), is developed in conjunction with advanced mathematical regularization to investigate radiation dose reduction in x-ray CT. The method is experimentally implemented on fan-beam CT and evaluated as a function of imaging dose on a series of image quality phantoms and anonymous pediatric patient data sets. Numerical simulation experiments are also performed to explore the extension of EST to helical cone-beam geometry. Methods: EST is a Fourier based iterative algorithm, which iterates back and forth between real and Fourier space utilizing the algebraically exact pseudopolar fast Fourier transform (PPFFT). In each iteration, physical constraints and mathematical regularization are applied in real space, while the measured data are enforced in Fourier space. The algorithm is automatically terminated when a proposed termination criterion is met. Experimentally, fan-beam projections were acquired by the Siemens z-flying focal spot technology, and subsequently interleaved and rebinned to a pseudopolar grid. Image quality phantoms were scanned at systematically varied mAs settings, reconstructed by EST and conventional reconstruction methods such as filtered back projection (FBP), and quantified using metrics including resolution, signal-to-noise ratios (SNRs), and contrast-to-noise ratios (CNRs). Pediatric data sets were reconstructed at their original acquisition settings and additionally simulated to lower dose settings for comparison and evaluation of the potential for radiation dose reduction. Numerical experiments were conducted to quantify EST and other iterative methods in terms of image quality and computation time. The extension of EST to helical cone-beam CT was implemented by using the advanced single-slice rebinning (ASSR) method. Results: Based on the phantom and pediatric patient fan-beam CT data, it is demonstrated that EST reconstructions with the lowest scanner flux setting of 39 mAs produce comparable image quality, resolution, and contrast relative to FBP with the 140 mAs flux setting. Compared to the algebraic reconstruction technique and the expectation maximization statistical reconstruction algorithm, a significant reduction in computation time is achieved with EST. Finally, numerical experiments on helical cone-beam CT data suggest that the combination of EST and ASSR produces reconstructions with higher image quality and lower noise than the Feldkamp Davis and Kress (FDK) method and the conventional ASSR approach. Conclusions: A Fourier-based iterative method has been applied to the reconstruction of fan-bean CT data with reduced x-ray fluence. This method incorporates advantageous features in both real and Fourier space iterative schemes: using a fast and algebraically exact method to calculate forward projection, enforcing the measured data in Fourier space, and applying physical constraints and flexible regularization in real space. Our results suggest that EST can be utilized for radiation dose reduction in x-ray CT via the readily implementable technique of lowering mAs settings. Numerical experiments further indicate that EST requires less computation time than several other iterative algorithms and can, in principle, be extended to helical cone-beam geometry in combination with the ASSR method.
60 APPLIED LIFE SCIENCES↗
Broken conformal window
We show that near the edges of the conformal window of supersymmetric SU(Nc) QCD, perturbed by Anomaly Mediated Supersymmetry Breaking (AMSB), chiral symmetry can be broken depending on the initial conditions of the RG flow. We do so by perturbatively expanding around Banks-Zaks fixed points and taking advantage of Seiberg duality. Interpolating between the edges of the conformal window, we predict that non-supersymmetric QCD breaks chiral symmetry up to N f ≤ 3N c − 1, while we cannot say anything definitive for N f ≥ 3N c at this moment.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Structural optimization - Challenges and opportunities
A review of developments in structural optimization techniques and their interface with growing computer capabilities is presented. Structural design steps comprise functional definition of an object, an evaluation phase wherein external influences are quantified, selection of the design concept, material, object geometry, and the internal layout, and quantification of the physical characteristics. Optimization of a fully stressed design is facilitated by use of nonlinear mathematical programming which permits automated definition of the physics of a problem. Design iterations terminate when convergence is acquired between mathematical and physical criteria. A constrained minimum algorithm has been formulated using an Augmented Lagrangian approach and a generalized reduced gradient to obtain fast convergence. Various approximation techniques are mentioned. The synergistic application of all the methods surveyed requires multidisciplinary teamwork during a design effort.
Sobieszczanski-Sobieski, J.↗
Structural optimization: Challenges and opportunities
A review of developments in structural optimization techniques and their interface with growing computer capabilities is presented. Structural design steps comprise functional definition of an object, an evaluation phase wherein external influences are quantified, selection of the design concept, material, object geometry, and the internal layout, and quantification of the physical characteristics. Optimization of a fully stressed design is facilitated by use of nonlinear mathematical programming which permits automated definition of the physics of a problem. Design iterations terminate when convergence is acquired between mathematical and physical criteria. A constrained minimum algorithm has been formulated using an Augmented Lagrangian approach and a generalized reduced gradient to obtain fast convergence. Various approximation techniques are mentioned. The synergistic application of all the methods surveyed requires multidisciplinary teamwork during a design effort.
Sobieszczanski-Sobieski, J.↗
Physical scales in the Wigner-Boltzmann equation
The Wigner-Boltzmann equation provides the Wigner single particle theory with interactions with bosonic degrees of freedom associated with harmonic oscillators, such as phonons in solids. Quantum evolution is an interplay of two transport modes, corresponding to the common coherent particle-potential processes, or to the decoherence causing scattering due to the oscillators. Which evolution mode will dominate depends on the scales of the involved physical quantities. A dimensionless formulation of the Wigner-Boltzmann equation is obtained, where these scales appear as dimensionless strength parameters. A notion called scaling theorem is derived, linking the strength parameters to the coupling with the oscillators. It is shown that an increase of this coupling is equivalent to a reduction of both the strength of the electric potential, and the coherence length. Secondly, the existence of classes of physically different, but mathematically equivalent setups of the Wigner-Boltzmann evolution is demonstrated. - Highlights: Black-Right-Pointing-Pointer Dimensionless parameters determine the ratio of quantum or classical WB evolution. Black-Right-Pointing-Pointer The scaling theorem evaluates the decoherence effect due to scattering. Black-Right-Pointing-Pointer Evolution processes are grouped into classes of equivalence.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Geometric uncertainty relation for mixed quantum states
In this paper we use symplectic reduction in an Uhlmann bundle to construct a principal fiber bundle over a general space of unitarily equivalent mixed quantum states. The bundle, which generalizes the Hopf bundle for pure states, gives in a canonical way rise to a Riemannian metric and a symplectic structure on the base space. With these we derive a geometric uncertainty relation for observables acting on quantum systems in mixed states. We also give a geometric proof of the classical Robertson-Schrödinger uncertainty relation, and we compare the two. They turn out not to be equivalent, because of the multiple dimensions of the gauge group for general mixed states. We give examples of observables for which the geometric relation provides a stronger estimate than that of Robertson and Schrödinger, and vice versa.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
SASHIMI-SIDM: semi-analytical subhalo modelling for self-interacting dark matter at sub-galactic scales
We combine the semi-analytical structure formation model, SASHIMI, which predicts subhalo populations in collisionless, cold dark matter (CDM), with a parametric model that maps CDM halos to self-interacting dark matter (SIDM) halos. The resulting model, SASHIMI-SIDM, generates SIDM subhalo populations down to sub-galactic mass scales, for an arbitrary input cross section, in minutes. We show that SASHIMI-SIDM agrees with SIDM subhalo populations from high-resolution cosmological zoom-in simulations in resolved regimes. Crucially, we predict that the fraction of core-collapsed subhalos peaks at a mass scale determined by the input SIDM cross section and decreases toward higher halo masses, consistent with the predictions of gravothermal models and cosmological simulations. For the first time, we also show that the core-collapsed fraction decreases toward lower halo masses. While the dependence of the collapse time on mass and concentration implies such behaviour, our semi-analytical approach allows us to quantify and illustrate this trend clearly across the full mass spectrum of subhalos, including for subhalo masses below the resolution limit of any current cosmological SIDM simulation. As a proof of principle, we apply SASHIMI-SIDM to predict the boost to the local dark matter density and annihilation rate from core-collapsed SIDM subhalos, which can be enhanced relative to CDM by an order of magnitude for viable SIDM models. Thus, SASHIMI-SIDM provides an efficient and reliable tool for scanning SIDM parameter space and testing it with astrophysical observations. The code is publicly available at https://github.com/shinichiroando/sashimi-si.
cosmological simulations↗
Multiple axions save high-scale inflation
Many models of dark matter QCD axion require inflation at a scale H inf ≲ 10 6 GeV ,which precludes a detectable tensor mode fluctuation. This constraint arises because the domain wall problem forces the Peccei-Quinn symmetry to be broken during inflation and the axions to be produced by the misalignment mechanism. Here, we point out that theories with multiple axions can evade this constraint, thereby allowing high-scale inflation with detectable tensor mode. The only requirement is a condition on the anomaly coefficients that ensures a unique minimum of the axion potential without fine-tuning or the introduction of small parameters.
Kondo, Dan [University of Tokyo (Japan). Kavli Ins↗
Near-SUSY to non-SUSY crossover
Gauge theories can be solved exactly slightly away from the supersymmetric (SUSY) limit softly broken by anomaly mediation when the size of SUSY breaking is much smaller than the dynamical scale (𝑚 ≪ Λ). We show empirical evidence that the near-SUSY limit is continuously connected to the non-SUSY limit (𝑚 ≫ Λ) in SU(𝑁 𝑐 ) gauge theories with 𝑁 𝑓 quarks in the fundamental representation. The evidence includes the behavior of quark bilinear condensate and gluon condensates, light hadron spectra, and consistency with the large 𝑁 𝑐 limit including the Witten-Veneziano relations between the topological susceptibility and 𝑚 𝜂′ and 𝑚 𝜋 . In addition, we present new predictions when 𝑁 𝑓 /𝑁 𝑐 ≳ 𝑂(1).
Effective field theory↗
Uncertainty quantification and propagation in nuclear density functional theory
Nuclear density functional theory (DFT) is one of the main theoretical tools used to study the properties of heavy and superheavy elements, or to describe the structure of nuclei far from stability. While on-going eff orts seek to better root nuclear DFT in the theory of nuclear forces, energy functionals remain semi-phenomenological constructions that depend on a set of parameters adjusted to experimental data in fi nite nuclei. In this study, we review recent eff orts to quantify the related uncertainties, and propagate them to model predictions. In particular, we cover the topics of parameter estimation for inverse problems, statistical analysis of model uncertainties and Bayesian inference methods. Illustrative examples are taken from the literature.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Design of quantum optical experiments with logic artificial intelligence
Logic Artificial Intelligence (AI) is a subfield of AI where variables can take two defined arguments, True or False, and are arranged in clauses that follow the rules of formal logic. Several problems that span from physical systems to mathematical conjectures can be encoded into these clauses and solved by checking their satisfiability (SAT). In contrast to machine learning approaches where the results can be approximations or local minima, Logic AI delivers formal and mathematically exact solutions to those problems. In this work, we propose the use of logic AI for the design of optical quantum experiments. We show how to map into a SAT problem the experimental preparation of an arbitrary quantum state and propose a logic-based algorithm, called Klaus, to find an interpretable representation of the photonic setup that generates it. We compare the performance of Klaus with the state-of-the-art algorithm for this purpose based on continuous optimization. We also combine both logic and numeric strategies to find that the use of logic AI significantly improves the resolution of this problem, paving the path to developing more formal-based approaches in the context of quantum physics experiments.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Laboratory Directed Research and Development FY2010 Annual Report
A premier applied-science laboratory, Lawrence Livermore National Laboratory (LLNL) has at its core a primary national security mission - to ensure the safety, security, and reliability of the nation's nuclear weapons stockpile without nuclear testing, and to prevent and counter the spread and use of weapons of mass destruction: nuclear, chemical, and biological. The Laboratory uses the scientific and engineering expertise and facilities developed for its primary mission to pursue advanced technologies to meet other important national security needs - homeland defense, military operations, and missile defense, for example - that evolve in response to emerging threats. For broader national needs, LLNL executes programs in energy security, climate change and long-term energy needs, environmental assessment and management, bioscience and technology to improve human health, and for breakthroughs in fundamental science and technology. With this multidisciplinary expertise, the Laboratory serves as a science and technology resource to the U.S. government and as a partner with industry and academia. This annual report discusses the following topics: (1) Advanced Sensors and Instrumentation; (2) Biological Sciences; (3) Chemistry; (4) Earth and Space Sciences; (5) Energy Supply and Use; (6) Engineering and Manufacturing Processes; (7) Materials Science and Technology; Mathematics and Computing Science; (8) Nuclear Science and Engineering; and (9) Physics.
36 MATERIALS SCIENCE↗