Optimal low-thrust trajectories to asteroids through an algorithm based on differential dynamic programming

In this paper an optimisation algorithm based on Differential Dynamic Programming is applied to the design of rendezvous and fly-by trajectories to near Earth objects. Differential dynamic programming is a successive approximation technique that computes a feedback control law in correspondence of a fixed number of decision times. In this way the high dimensional problem characteristic of low-thrust optimisation is reduced into a series of small dimensional problems. The proposed method exploits the stage-wise approach to incorporate an adaptive refinement of the discretisation mesh within the optimisation process. A particular interpolation technique was used to preserve the feedback nature of the control law, thus improving robustness against some approximation errors introduced during the adaptation process. The algorithm implements global variations of the control law, which ensure a further increase in robustness. The results presented show how the proposed approach is capable of fully exploiting the multi-body dynamics of the problem; in fact, in one of the study cases, a fly-by of the Earth is scheduled, which was not included in the first guess solution

Second order gauge invariant gravitational perturbations of a Kerr black hole

We investigate higher than the first order gravitational perturbations in the Newman-Penrose formalism. Equations for the Weyl scalar $\psi_4,$ representing outgoing gravitational radiation, can be uncoupled into a single wave equation to any perturbative order. For second order perturbations about a Kerr black hole, we prove the existence of a first and second order gauge (coordinates) and tetrad invariant waveform, $\psi_I$, by explicit construction. This waveform is formed by the second order piece of $\psi_4$ plus a term, quadratic in first order perturbations, chosen to make $\psi_I$ totally invariant and to have the appropriate behavior in an asymptotically flat gauge. $\psi_I$ fulfills a single wave equation of the form ${\cal T}\psi_I=S,$ where ${\cal T}$ is the same wave operator as for first order perturbations and $S$ is a source term build up out of (known to this level) first order perturbations. We discuss the issues of imposition of initial data to this equation, computation of the energy and momentum radiated and wave extraction for direct comparison with full numerical approaches to solve Einstein equations.Comment: 19 pages, REVTEX. Some misprints corrected and changes to improve presentation. Version to appear in PR

Spontaneous Coherence and Collective Modes in Double-Layer Quantum Dot Systems

We study the ground state and the collective excitations of parabolically-confined double-layer quantum dot systems in a strong magnetic field. We identify parameter regimes where electrons form maximum density droplet states, quantum-dot analogs of the incompressible states of the bulk integer quantum Hall effect. In these regimes the Hartree-Fock approximation and the time-dependent Hartree-Fock approximations can be used to describe the ground state and collective excitations respectively. We comment on the relationship between edge excitations of dots and edge magneto-plasmon excitations of bulk double-layer systems.Comment: 20 pages (figures included) and also available at http://fangio.magnet.fsu.edu/~jhu/Paper/qdot_cond.ps, replaced to fix figure

Double-Layer Systems at Zero Magnetic Field

We investigate theoretically the effects of intralayer and interlayer exchange in biased double-layer electron and hole systems, in the absence of a magnetic field. We use a variational Hartree-Fock-like approximation to analyze the effects of layer separation, layer density, tunneling, and applied gate voltages on the layer densities and on interlayer phase coherence. In agreement with earlier work, we find that for very small layer separations and low layer densities, an interlayer-correlated ground state possessing spontaneous interlayer coherence (SILC) is obtained, even in the absence of interlayer tunneling. In contrast to earlier work, we find that as a function of total density, there exist four, rather than three, distinct noncrystalline phases for balanced double-layer systems without interlayer tunneling. The newly identified phase exists for a narrow range of densities and has three components and slightly unequal layer densities, with one layer being spin polarized, and the other unpolarized. An additional two-component phase is also possible in the presence of sufficiently strong bias or tunneling. The lowest-density SILC phase is the fully spin- and pseudospin-polarized ``one-component'' phase discussed by Zheng {\it et al.} [Phys. Rev. B {\bf 55}, 4506 (1997)]. We argue that this phase will produce a finite interlayer Coulomb drag at zero temperature due to the SILC. We calculate the particle densities in each layer as a function of the gate voltage and total particle density, and find that interlayer exchange can reduce or prevent abrupt transfers of charge between the two layers. We also calculate the effect of interlayer exchange on the interlayer capacitance.Comment: 35 pages, 19 figures included. To appear in PR

Gravitational radiation from gamma-ray bursts as observational opportunities for LIGO and VIRGO

Gamma-ray bursts are believed to originate in core-collapse of massive stars. This produces an active nucleus containing a rapidly rotating Kerr black hole surrounded by a uniformly magnetized torus represented by two counter-oriented current rings. We quantify black hole spin-interactions with the torus and charged particles along open magnetic flux-tubes subtended by the event horizon. A major output of Egw=4e53 erg is radiated in gravitational waves of frequency fgw=500 Hz by a quadrupole mass-moment in the torus. Consistent with GRB-SNe, we find (i) Ts=90s (tens of s, Kouveliotou et al. 1993), (ii) aspherical SNe of kinetic energy Esn=2e51 erg (2e51 erg in SN1998bw, Hoeflich et al. 1999) and (iii) GRB-energies Egamma=2e50 erg (3e50erg in Frail et al. 2001). GRB-SNe occur perhaps about once a year within D=100Mpc. Correlating LIGO/Virgo detectors enables searches for nearby events and their spectral closure density 6e-9 around 250Hz in the stochastic background radiation in gravitational waves. At current sensitivity, LIGO-Hanford may place an upper bound around 150MSolar in GRB030329. Detection of Egw thus provides a method for identifying Kerr black holes by calorimetry.Comment: to appear in PRD, 49

Magnetic Field Amplification in Galaxy Clusters and its Simulation

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.Comment: 60 pages, 19 Figure

On the selection of AGN neutrino source candidates for a source stacking analysis with neutrino telescopes

The sensitivity of a search for sources of TeV neutrinos can be improved by grouping potential sources together into generic classes in a procedure that is known as source stacking. In this paper, we define catalogs of Active Galactic Nuclei (AGN) and use them to perform a source stacking analysis. The grouping of AGN into classes is done in two steps: first, AGN classes are defined, then, sources to be stacked are selected assuming that a potential neutrino flux is linearly correlated with the photon luminosity in a certain energy band (radio, IR, optical, keV, GeV, TeV). Lacking any secure detailed knowledge on neutrino production in AGN, this correlation is motivated by hadronic AGN models, as briefly reviewed in this paper. The source stacking search for neutrinos from generic AGN classes is illustrated using the data collected by the AMANDA-II high energy neutrino detector during the year 2000. No significant excess for any of the suggested groups was found.Comment: 43 pages, 12 figures, accepted by Astroparticle Physic

Dynamical properties of liquid Al near melting. An orbital-free molecular dynamics study

The static and dynamic structure of liquid Al is studied using the orbital free ab-initio molecular dynamics method. Two thermodynamic states along the coexistence line are considered, namely T = 943 K and 1323 K for which X-ray and neutron scattering data are available. A new kinetic energy functional, which fulfills a number of physically relevant conditions is employed, along with a local first principles pseudopotential. In addition to a comparison with experiment, we also compare our ab-initio results with those obtained from conventional molecular dynamics simulations using effective interionic pair potentials derived from second order pseudopotential perturbation theory.Comment: 15 pages, 12 figures, 2 tables, submitted to PR

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Measurement of D*+/- meson production in jets from pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

This paper reports a measurement of D*+/- meson production in jets from proton-proton collisions at a center-of-mass energy of sqrt(s) = 7 TeV at the CERN Large Hadron Collider. The measurement is based on a data sample recorded with the ATLAS detector with an integrated luminosity of 0.30 pb^-1 for jets with transverse momentum between 25 and 70 GeV in the pseudorapidity range |eta| < 2.5. D*+/- mesons found in jets are fully reconstructed in the decay chain: D*+ -> D0pi+, D0 -> K-pi+, and its charge conjugate. The production rate is found to be N(D*+/-)/N(jet) = 0.025 +/- 0.001(stat.) +/- 0.004(syst.) for D*+/- mesons that carry a fraction z of the jet momentum in the range 0.3 < z < 1. Monte Carlo predictions fail to describe the data at small values of z, and this is most marked at low jet transverse momentum.Comment: 10 pages plus author list (22 pages total), 5 figures, 1 table, matches published version in Physical Review