Simulation of granular soil behaviour using the bullet physics library

A physics engine is computer software which provides a simulation of certain physical systems, such as rigid body dynamics, soft body dynamics and fluid dynamics. Physics engines were firstly developed for using in animation and gaming industry ; nevertheless, due to fast calculation speed they are attracting more and more attetion from researchers of the engineering fields. Since physics engines are capable of performing fast calculations on multibody rigid dynamic systems, soil particles can be modeled as distinct rigid bodies. However, up to date, it is not clear to what extent they perform accurately in modeling soil behaviour from a geotechnical viewpoint. To investigate this, examples of pluviation and vibration-induced desification were simulated using the physics engine called Bullet physics library. In order to create soil samples, first, randomly shaped polyhedrons, representing gravels, were generated using the Voronoi tessellation approach. Then, particles were pluviated through a funnel into a cylinder. Once the soil particles settled in a static state, the cylinder was subjected to horizontal sinusoidal vibration for a period of 20 seconds. The same procedure for sample perparation was performed in the laboratory. The results of pluviation and vibration tests weere recorded and compared to those of simulations. A good agreement has been found between the results of simulations and laboratory tests. The findings in this study reinforce the idea that physics engines can be employed as a geotechnical engineering simulation tool

Skyrme functional from a three-body pseudo-potential of second-order in gradients. Formalism for central terms

In one way or the other, all modern parametrizations of the nuclear energy density functional (EDF) do not respect the exchange symmetry associated with Pauli's principle. It has been recently shown that this practice jeopardizes multi-reference (MR) EDF calculations by contaminating the energy with spurious self-interactions that, for example, lead to finite steps or even divergences when plotting it as a function of collective coordinates. As of today, the only viable option to bypass these pathologies is to rely on EDF kernels that enforce Pauli's principle from the outset by strictly and exactly deriving from a genuine, i.e. density-independent, Hamilton operator. We wish to develop the most general Skyrme-like EDF parametrization containing linear, bilinear and trilinear terms in the density matrices with up to two gradients, under the key constraint that it derives strictly from an effective Hamilton operator. The most general three-body Skyrme-like pseudo-potential containing up to two gradient operators is constructed to generate the trilinear part. The present study is limited to central terms. Spin-orbit and tensor will be addressed in a forthcoming paper. (See paper for full abstract)Comment: 38 pages revtex, no figur

Pairing correlations beyond the mean field

We discuss dynamical pairing correlations in the context of configuration mixing of projected self-consistent mean-field states, and the origin of a divergence that might appear when such calculations are done using an energy functional in the spirit of a naive generalized density functional theory.Comment: Proceedings of the XIII Nuclear Physics Workshop ``Maria and Pierre Curie'' on ``Pairing and beyond - 50 years of the BCS model'', held at Kazimierz Dolny, Poland, September 27 - October 1, 2006. Int. J. Mod. Phys. E, in prin

Time-Reversal Violating Schiff Moment of 225Ra

We use the Skyrme-Hartree-Fock method, allowing all symmetries to be broken, to calculate the time-reversal-violating nuclear Schiff moment (which induces atomic electric dipole moments) in the octupole-deformed nucleus 225Ra. Our calculation includes several effects neglected in earlier work, including self consistency and polarization of the core by the last nucleon. We confirm that the Schiff moment is large compared to those of reflection-symmetric nuclei, though ours is generally a few times smaller than recent estimates.Comment: Typos corrected, references added, minor changesin text. Version to appear in PRC. 10 pages, 4 figure

Spatially resolved spectroscopy of Coma cluster early-type galaxies IV. Completing the dataset

The long-slit spectra obtained along the minor axis, offset major axis and diagonal axis are presented for 12 E and S0 galaxies of the Coma cluster drawn from a magnitude-limited sample studied before. The rotation curves, velocity dispersion profiles and the H_3 and H_4 coefficients of the Hermite decomposition of the line of sight velocity distribution are derived. The radial profiles of the Hbeta, Mg, and Fe line strength indices are measured too. In addition, the surface photometry of the central regions of a subsample of 4 galaxies recently obtained with Hubble Space Telescope is presented. The data will be used to construct dynamical models of the galaxies and study their stellar populations.Comment: 40 pages, 7 figures, 6 tables. Accepted for publication in ApJ

Stellar Orbits and the Interstellar Gas Temperature in Elliptical Galaxies

We draw attention to the close relationship between the anisotropy parameter beta(r) for stellar orbits in elliptical galaxies and the temperature profile T(r) of the hot interstellar gas. For nearly spherical galaxies the gas density can be accurately determined from X-ray observations and the stellar luminosity density can be accurately found from the optical surface brightness. The Jeans equation and hydrostatic equilibrium establish a connection between beta(r) and T(r) that must be consistent with the observed stellar velocity dispersion. Purely optical observations of the bright elliptical galaxy NGC 4472 indicate beta(r) < 0.35 within the effective radius. However, the X-ray gas temperature profile T(r) for NGC 4472 requires significantly larger anisotropy, beta = 0.6 - 0.7, about twice the optical value. This strong preference for radial stellar orbits must be understood in terms of the formation history of massive elliptical galaxies. Conversely, if the smaller, optically determined anisotropy is indeed correct, we are led to the important conclusion that the temperature profile T(r) of the hot interstellar gas in NGC 4472 must differ from that indicated by X-ray observations, or that the hot gas is not in hydrostatic equilibrium.Comment: 6 pages (emulateapj5) with 4 figures; accepted by The Astrophysical Journa

PT-Symmetry Quantum Electrodynamics--PTQED

The construction of $\mathcal{PT}$-symmetric quantum electrodynamics is reviewed. In particular, the massless version of the theory in 1+1 dimensions (the Schwinger model) is solved. Difficulties with unitarity of the $S$-matrix are discussed.Comment: 11 pages, 1 figure, contributed to Proceedings of 6th International Workshop on Pseudo-Hermitian Hamiltonians in Quantum Physic

Extending PT symmetry from Heisenberg algebra to E2 algebra

The E2 algebra has three elements, J, u, and v, which satisfy the commutation relations [u,J]=iv, [v,J]=-iu, [u,v]=0. We can construct the Hamiltonian H=J^2+gu, where g is a real parameter, from these elements. This Hamiltonian is Hermitian and consequently it has real eigenvalues. However, we can also construct the PT-symmetric and non-Hermitian Hamiltonian H=J^2+igu, where again g is real. As in the case of PT-symmetric Hamiltonians constructed from the elements x and p of the Heisenberg algebra, there are two regions in parameter space for this PT-symmetric Hamiltonian, a region of unbroken PT symmetry in which all the eigenvalues are real and a region of broken PT symmetry in which some of the eigenvalues are complex. The two regions are separated by a critical value of g.Comment: 8 pages, 7 figure

Hydrogeological investigations in the Pampa of Argentina

The author has identified the following significant results. Satellite imagery in combination with ground investigations allows identification and delineation of differences in the conditions of the near surface ground water (depth to ground water, salinity). The degree of precision achieved is greater than that obtainable by conventional ground survey methods alone

Kinematic Structure of Merger Remnants

We use numerical simulations to study the kinematic structure of remnants formed from mergers of equal-mass disk galaxies. In particular, we show that remnants of dissipational mergers, which include the radiative cooling of gas, star formation, feedback from supernovae, and the growth of supermassive black holes, are smaller, rounder, have, on average, a larger central velocity dispersion, and show significant rotation compared to remnants of dissipationless mergers. The increased rotation speed of dissipational remnants owes its origin to star formation that occurs in the central regions during the galaxy merger. We have further quantified the anisotropy, three-dimensional shape, minor axis rotation, and isophotal shape of each merger remnant, finding that dissipational remnants are more isotropic, closer to oblate, have the majority of their rotation along their major axis, and are more disky than dissipationless remnants. Individual remnants display a wide variety of kinematic properties. A large fraction of the dissipational remnants are oblate isotropic rotators. Many dissipational, and all of the dissipationless, are slowly rotating and anisotropic. The remnants of gas-rich major mergers can well-reproduce the observed distribution of projected ellipticities, rotation parameter (V/\sigma)*, kinematic misalignments, Psi, and isophotal shapes. The dissipationless remnants are a poor match to this data. Our results support the merger hypothesis for the origin of low-luminosity elliptical galaxies provided that the progenitor disks are sufficiently gas-rich, however our remnants are a poor match to the bright ellipticals that are slowly rotating and uniformly boxy.Comment: 22 pages, 17 figures, accepted to Ap