An Unsplit, Cell-Centered Godunov Method for Ideal MHD

We present a second-order Godunov algorithm for multidimensional, ideal MHD. Our algorithm is based on the unsplit formulation of Colella (J. Comput. Phys. vol. 87, 1990), with all of the primary dependent variables centered at the same location. To properly represent the divergence-free condition of the magnetic fields, we apply a discrete projection to the intermediate values of the field at cell faces, and apply a filter to the primary dependent variables at the end of each time step. We test the method against a suite of linear and nonlinear tests to ascertain accuracy and stability of the scheme under a variety of conditions. The test suite includes rotated planar linear waves, MHD shock tube problems, low-beta flux tubes, and a magnetized rotor problem. For all of these cases, we observe that the algorithm is second-order accurate for smooth solutions, converges to the correct weak solution for problems involving shocks, and exhibits no evidence of instability or loss of accuracy due to the possible presence of non-solenoidal fields.Comment: 37 Pages, 9 Figures, submitted to Journal of Computational Physic

Long-Lived Non-Equilibrium Interstitial-Solid-Solutions in Binary Mixtures

We perform particle resolved experimental studies on the heterogeneous crystallisation process of two compo- nent mixtures of hard spheres. The components have a size ratio of 0.39. We compared these with molecular dynamics simulations of homogenous nucleation. We find for both experiments and simulations that the final assemblies are interstitial solid solutions, where the large particles form crystalline close-packed lattices, whereas the small particles occupy random interstitial sites. This interstitial solution resembles that found at equilibrium when the size ratios are 0.3 [Filion et al., Phys. Rev. Lett. 107, 168302 (2011)] and 0.4 [Filion, PhD Thesis, Utrecht University (2011)]. However, unlike these previous studies, for our system sim- ulations showed that the small particles are trapped in the octahedral holes of the ordered structure formed by the large particles, leading to long-lived non-equilibrium structures in the time scales studied and not the equilibrium interstitial solutions found earlier. Interestingly, the percentage of small particles in the crystal formed by the large ones rapidly reaches a maximum of around 14% for most of the packing fractions tested, unlike previous predictions where the occupancy of the interstitial sites increases with the system concentration. Finally, no further hopping of the small particles was observed

Advection-dominated Inflow/Outflows from Evaporating Accretion Disks

In this Letter we investigate the properties of advection-dominated accretion flows (ADAFs) fed by the evaporation of a Shakura-Sunyaev accretion disk (SSD). In our picture the ADAF fills the central cavity evacuated by the SSD and extends beyond the transition radius into a coronal region. We find that, because of global angular momentum conservation, a significant fraction of the hot gas flows away from the black hole forming a transsonic wind, unless the injection rate depends only weakly on radius (if $r^2\dot\sigma\propto r^{-\xi}$, $\xi< 1/2$). The Bernoulli number of the inflowing gas is negative if the transition radius is $\lesssim 100$ Schwarzschild radii, so matter falling into the hole is gravitationally bound. The ratio of inflowing to outflowing mass is $\approx 1/2$, so in these solutions the accretion rate is of the same order as in standard ADAFs and much larger than in advection-dominated inflow/outflow models (ADIOS). The possible relevance of evaporation-fed solutions to accretion flows in black hole X-ray binaries is briefly discussed.Comment: 5 pages Latex with 2 ps figures. Accepted for publication in ApJ Letter

Temporal characterization of individual harmonics of an attosecond pulse train by THz streaking

We report on the global temporal pulse characteristics of individual harmonics in an attosecond pulse train by means of photo-electron streaking in a strong low-frequency transient. The scheme allows direct retrieval of pulse durations and first order chirp of individual harmonics without the need of temporal scanning. The measurements were performed using an intense THz field generated by tilted phase front technique in LiNbO_3 . Pulse properties for harmonics of order 23, 25 and 27 show that the individual pulse durations and linear chirp are decreasing by the harmonic order

Supersymmetrization of the Radiation Damping

We construct a supersymmetrized version of the model to the radiation damping \cite{03} introduced by the present authors \cite{ACWF}. We dicuss its symmetries and the corresponding conserved Noether charges. It is shown this supersymmetric version provides a supersymmetric generalization of the Galilei algebra obtained in \cite{ACWF}. We have shown that the supersymmetric action can be splited into dynamically independent external and internal sectors.Comment: 9 page

Extensive chaos in Rayleigh-Bénard convection

Using large-scale numerical calculations we explore spatiotemporal chaos in Rayleigh-Bénard convection for experimentally relevant conditions. We calculate the spectrum of Lyapunov exponents and the Lyapunov dimension describing the chaotic dynamics of the convective fluid layer at constant thermal driving over a range of finite system sizes. Our results reveal that the dynamics of fluid convection is truly chaotic for experimental conditions as illustrated by a positive leading-order Lyapunov exponent. We also find the chaos to be extensive over the range of finite-sized systems investigated as indicated by a linear scaling between the Lyapunov dimension of the chaotic attractor and the system size

Reversing a granular flow on a vibratory conveyor

Experimental results are presented on the transport properties of granular materials on a vibratory conveyor. For circular oscillations of the shaking trough a non-monotonous dependence of the transport velocity on the normalized acceleration is observed. Two maxima are separated by a regime, where the granular flow is much slower and, in a certain driving range, even reverses its direction. A similar behavior is found for a single solid body with a low coefficient of restitution, whereas an individual glass bead of 1 mm diameter is propagated in the same direction for all accelerations.Comment: 4 pages, 5 figures, submitted to Applied Physics Letter

Decoherence in a Two Slit Diffraction Experiment with Massive Particles

Matter-wave interferometry has been largely studied in the last few years. Usually, the main problem in the analysis of the diffraction experiments is to establish the causes for the loss of coherence observed in the interference pattern. In this work, we use different type of environmental couplings to model a two slit diffraction experiment with massive particles. For each model, we study the effects of decoherence on the interference pattern and define a visibility function that measures the loss of contrast of the interference fringes on a distant screen. Finally, we apply our results to the experimental reported data on massive particles $C_{70}$.Comment: 6 pages, 3 figure

Uncertainties in water chemistry in disks: An application to TW Hya

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.Comment: accepted for publication in A&