Nuclear three-body problem in the complex energy plane: Complex-Scaling-Slater method

The physics of open quantum systems is an interdisciplinary area of research. The nuclear "openness" manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape towards the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a Complex Scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow Shell Model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6He viewed as an $\alpha$+ n + n cluster system. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly-bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method, if it can be applied, can provide an additional information about partial energy widths associated with individual thresholds.Comment: 15 pages, 16 figure

Depolarization volume and correlation length in the homogenization of anisotropic dielectric composites

In conventional approaches to the homogenization of random particulate composites, both the distribution and size of the component phase particles are often inadequately taken into account. Commonly, the spatial distributions are characterized by volume fraction alone, while the electromagnetic response of each component particle is represented as a vanishingly small depolarization volume. The strong-permittivity-fluctuation theory (SPFT) provides an alternative approach to homogenization wherein a comprehensive description of distributional statistics of the component phases is accommodated. The bilocally-approximated SPFT is presented here for the anisotropic homogenized composite which arises from component phases comprising ellipsoidal particles. The distribution of the component phases is characterized by a two-point correlation function and its associated correlation length. Each component phase particle is represented as an ellipsoidal depolarization region of nonzero volume. The effects of depolarization volume and correlation length are investigated through considering representative numerical examples. It is demonstrated that both the spatial extent of the component phase particles and their spatial distributions are important factors in estimating coherent scattering losses of the macroscopic field.Comment: Typographical error in eqn. 16 in WRM version is corrected in arxiv versio

Electron-phonon coupling in the C60 fullerene within the many-body GW approach

We study the electron-phonon coupling in the C60 fullerene within the first-principles GW approach, focusing on the lowest unoccupied t1u three-fold electronic state which is relevant for the superconducting transition in electron doped fullerides. It is shown that the strength of the coupling is significantly enhanced as compared to standard density functional theory calculations with (semi)local functionals, with a 48% increase of the electron-phonon potential Vep. The calculated GW value for the contribution from the Hg modes of 93 meV comes within 4% of the most recent experimental values. The present results call for a reinvestigation of previous density functional based calculations of electron-phonon coupling in covalent systems in general.Comment: 4 pages, 0 figur

Random walk on surfaces with hyperbolic cusps

We consider the operator associated to a random walk on finite volume surfaces with hyperbolic cusps. We study the spectral gap (upper and lower bound) associated to this operator and deduce some rate of convergence of the iterated kernel towards its stationary distribution.Comment: 28 page

Speeding up simulations of relativistic systems using an optimal boosted frame

It can be computationally advantageous to perform computer simulations in a Lorentz boosted frame for a certain class of systems. However, even if the computer model relies on a covariant set of equations, it has been pointed out that algorithmic difficulties related to discretization errors may have to be overcome in order to take full advantage of the potential speedup. We summarize the findings, the difficulties and their solutions, and show that the technique enables simulations important to several areas of accelerator physics that are otherwise problematic, including self-consistent modeling in three-dimensions of laser wakefield accelerator stages at energies of 10 GeV and above.Comment: To be published in the proceedings of DPF-2009, Detroit, MI, July 2009, eConf C09072

Ekman layers and the damping of inertial r-modes in a spherical shell: application to neutron stars

Recently, eigenmodes of rotating fluids, namely inertial modes, have received much attention in relation to their destabilization when coupled to gravitational radiation within neutron stars. However, these modes have been known for a long time by fluid dynamicists. We give a short account of their history and review our present understanding of their properties. Considering the case of a spherical container, we then give the exact solution of the boundary (Ekman) layer flow associated with inertial r-modes and show that previous estimations all underestimated the dissipation by these layers. We also show that the presence of an inner core has little influence on this dissipation. As a conclusion, we compute the window of instability in the Temperature/rotation plane for a crusted neutron star when it is modeled by an incompressible fluid.Comment: 7 pages, 2 figures, revised version to appear in ApJ, March 1

Depolarization regions of nonzero volume in bianisotropic homogenized composites

In conventional approaches to the homogenization of random particulate composites, the component phase particles are often treated mathematically as vanishingly small, point-like entities. The electromagnetic responses of these component phase particles are provided by depolarization dyadics which derive from the singularity of the corresponding dyadic Green functions. Through neglecting the spatial extent of the depolarization region, important information may be lost, particularly relating to coherent scattering losses. We present an extension to the strong-property-fluctuation theory in which depolarization regions of nonzero volume and ellipsoidal geometry are accommodated. Therein, both the size and spatial distribution of the component phase particles are taken into account. The analysis is developed within the most general linear setting of bianisotropic homogenized composite mediums (HCMs). Numerical studies of the constitutive parameters are presented for representative examples of HCM; both Lorentz-reciprocal and Lorentz-nonreciprocal HCMs are considered. These studies reveal that estimates of the HCM constitutive parameters in relation to volume fraction, particle eccentricity, particle orientation and correlation length are all significantly influenced by the size of the component phase particles

Effects of Hyperbolic Rotation in Minkowski Space on the Modeling of Plasma Accelerators in a Lorentz Boosted Frame

Laser driven plasma accelerators promise much shorter particle accelerators but their development requires detailed simulations that challenge or exceed current capabilities. We report the first direct simulations of stages up to 1 TeV from simulations using a Lorentz boosted calculation frame resulting in a million times speedup, thanks to a frame boost as high as gamma=1300. Effects of the hyperbolic rotation in Minkowski space resulting from the frame boost on the laser propagation in the plasma is shown to be key in the mitigation of a numerical instability that was limiting previous attempts