Search results for " computational"
showing 10 items of 661 documents
A weakly random Universe?
2010
The cosmic microwave background (CMB) radiation is characterized by well-established scales, the 2.7 K temperature of the Planckian spectrum and the $10^{-5}$ amplitude of the temperature anisotropy. These features were instrumental in indicating the hot and equilibrium phases of the early history of the Universe and its large scale isotropy, respectively. We now reveal one more intrinsic scale in CMB properties. We introduce a method developed originally by Kolmogorov, that quantifies a degree of randomness (chaos) in a set of numbers, such as measurements of the CMB temperature in some region. Considering CMB as a composition of random and regular signals, we solve the inverse problem of …
Sensitivity to Initial Conditions in an Extended Activator--Inhibitor Model for the Formation of Patterns
2018
Despite simplicity, the synchronous cellular automaton [D.A. Young, Math. Biosci. 72, 51 (1984)] enables reconstructing basic features of patterns of skin. Our extended model allows studying the formatting of patterns and their temporal evolution also on the favourable and hostile environments. As a result, the impact of different types of an environment is accounted for the dynamics of patterns formation. The process is based on two diffusible morphogens, the short-range activator and the long-range inhibitor, produced by differentiated cells (DCs) represented as black pixels. For a neutral environment, the extended model reduces to the original one. However, even the reduced model is stat…
Locust: C++ software for simulation of RF detection
2019
The Locust simulation package is a new C++ software tool developed to simulate the measurement of time-varying electromagnetic fields using RF detection techniques. Modularity and flexibility allow for arbitrary input signals, while concurrently supporting tight integration with physics-based simulations as input. External signals driven by the Kassiopeia particle tracking package are discussed, demonstrating conditional feedback between Locust and Kassiopeia during software execution. An application of the simulation to the Project 8 experiment is described. Locust is publicly available at https://github.com/project8/locust_mc.
Numerical viscosity in simulations of the two-dimensional Kelvin-Helmholtz instability
2020
The Kelvin-Helmholtz instability serves as a simple, well-defined setup for assessing the accuracy of different numerical methods for solving the equations of hydrodynamics. We use it to extend our previous analysis of the convergence and the numerical dissipation in models of the propagation of waves and in the tearing-mode instability in magnetohydrodynamic models. To this end, we perform two-dimensional simulations with and without explicit physical viscosity at different resolutions. A comparison of the growth of the modes excited by our initial perturbations allows us to estimate the effective numerical viscosity of two spatial reconstruction schemes (fifth-order monotonicity preservin…
Transition to turbulence in toroidal pipes
2011
AbstractIncompressible flow in toroidal pipes of circular cross-section was investigated by three-dimensional, time-dependent numerical simulations using a finite volume method. The computational domain included a whole torus and was discretized by up to ${\ensuremath{\sim} }11. 4\ensuremath{\times} 1{0}^{6} $ nodes. Two curvatures $\delta $ (radius of the cross-section/radius of the torus), namely 0.3 and 0.1, were examined; a streamwise forcing term was imposed, and its magnitude was made to vary so that the bulk Reynolds number ranged between ${\ensuremath{\sim} }3500$ and ${\ensuremath{\sim} }14\hspace{0.167em} 700$. The results were processed by different techniques in order to confirm…
Calculation of the two-loop heavy-flavor contribution to Bhabha scattering
2008
We describe in detail the calculation of the two-loop corrections to the QED Bhabha scattering cross section due to the vacuum polarization by heavy fermions. Our approach eliminates one mass scale from the most challenging part of the calculation and allows us to obtain the corrections in a closed analytical form. The result is valid for arbitrary values of the heavy fermion mass and the Mandelstam invariants, as long as s,t,u >> m_e^2.
Brownian dynamics simulations with hard-body interactions: Spherical particles
2012
A novel approach to account for hard-body interactions in (overdamped) Brownian dynamics simulations is proposed for systems with non-vanishing force fields. The scheme exploits the analytically known transition probability for a Brownian particle on a one-dimensional half-line. The motion of a Brownian particle is decomposed into a component that is affected by hard-body interactions and into components that are unaffected. The hard-body interactions are incorporated by replacing the affected component of motion by the evolution on a half-line. It is discussed under which circumstances this approach is justified. In particular, the algorithm is developed and formulated for systems with spa…
Dynamical coexistence in moderately polydisperse hard-sphere glasses
2020
We perform extensive numerical simulations of a paradigmatic model glass former, the hard-sphere fluid with 10% polydispersity. We sample from the ensemble of trajectories with fixed observation time, whereby single trajectories are generated by event-driven molecular dynamics. We show that these trajectories can be characterized in terms of the local structure, and we find a dynamical-structural (active-inactive) phase transition between two dynamical phases: one dominated by liquidlike trajectories with a low degree of local order and one dominated by glassylike trajectories with a high degree of local order. We show that both phases coexist and are separated by a spatiotemporal interface…
Simulation of many-qubit quantum computation with matrix product states
2006
Matrix product states provide a natural entanglement basis to represent a quantum register and operate quantum gates on it. This scheme can be materialized to simulate a quantum adiabatic algorithm solving hard instances of a NP-Complete problem. Errors inherent to truncations of the exact action of interacting gates are controlled by the size of the matrices in the representation. The property of finding the right solution for an instance and the expected value of the energy are found to be remarkably robust against these errors. As a symbolic example, we simulate the algorithm solving a 100-qubit hard instance, that is, finding the correct product state out of ~ 10^30 possibilities. Accum…
Domain-wall excitations in the two-dimensional Ising spin glass
2018
The Ising spin glass in two dimensions exhibits rich behavior with subtle differences in the scaling for different coupling distributions. We use recently developed mappings to graph-theoretic problems together with highly efficient implementations of combinatorial optimization algorithms to determine exact ground states for systems on square lattices with up to $10\,000\times 10\,000$ spins. While these mappings only work for planar graphs, for example for systems with periodic boundary conditions in at most one direction, we suggest here an iterative windowing technique that allows one to determine ground states for fully periodic samples up to sizes similar to those for the open-periodic…