TURBULENCE AND THE FORMATION OF FILAMENTS, LOOPS, AND SHOCK FRONTS IN NGC 1275

NGC 1275, the central galaxy in the Perseus cluster, is the host of gigantic hot bipolar bubbles inflated by active galactic nucleus (AGN) jets observed in the radio as Perseus A. It presents a spectacular H alpha-emitting nebulosity surrounding NGC 1275, with loops and filaments of gas extending to over 50 kpc. The origin of the filaments is still unknown, but probably correlates with the mechanism responsible for the giant buoyant bubbles. We present 2.5 and three-dimensional magnetohydrodynamical (MHD) simulations of the central region of the cluster in which turbulent energy, possibly triggered by star formation and supernovae (SNe) explosions, is introduced. The simulations reveal that the turbulence injected by massive stars could be responsible for the nearly isotropic distribution of filaments and loops that drag magnetic fields upward as indicated by recent observations. Weak shell-like shock fronts propagating into the intracluster medium (ICM) with velocities of 100-500 km s(-1) are found, also resembling the observations. The isotropic outflow momentum of the turbulence slows the infall of the ICM, thus limiting further starburst activity in NGC 1275. As the turbulence is subsonic over most of the simulated volume, the turbulent kinetic energy is not efficiently converted into heat and additional heating is required to suppress the cooling flow at the core of the cluster. Simulations combining the MHD turbulence with the AGN outflow can reproduce the temperature radial profile observed around NGC 1275. While the AGN mechanism is the main heating source, the SNe are crucial to isotropize the energy distribution.

Recovery of the star formation history of the LMC from the VISTA survey of the Magellanic system

The VISTA near infrared survey of the Magellanic System (VMC) will provide deep YJK(s) photometry reaching stars in the oldest turn-off point throughout the Magellanic Clouds (MCs). As part of the preparation for the survey, we aim to access the accuracy in the star formation history (SFH) that can be expected from VMC data, in particular for the Large Magellanic Cloud (LMC). To this aim, we first simulate VMC images containing not only the LMC stellar populations but also the foreground Milky Way (MW) stars and background galaxies. The simulations cover the whole range of density of LMC field stars. We then perform aperture photometry over these simulated images, access the expected levels of photometric errors and incompleteness, and apply the classical technique of SFH-recovery based on the reconstruction of colour-magnitude diagrams (CMD) via the minimisation of a chi-squared-like statistics. We verify that the foreground MW stars are accurately recovered by the minimisation algorithms, whereas the background galaxies can be largely eliminated from the CMD analysis due to their particular colours and morphologies. We then evaluate the expected errors in the recovered star formation rate as a function of stellar age, SFR(t), starting from models with a known age-metallicity relation (AMR). It turns out that, for a given sky area, the random errors for ages older than similar to 0.4 Gyr seem to be independent of the crowding. This can be explained by a counterbalancing effect between the loss of stars from a decrease in the completeness and the gain of stars from an increase in the stellar density. For a spatial resolution of similar to 0.1 deg(2), the random errors in SFR(t) will be below 20% for this wide range of ages. On the other hand, due to the lower stellar statistics for stars younger than similar to 0.4 Gyr, the outer LMC regions will require larger areas to achieve the same level of accuracy in the SFR( t). If we consider the AMR as unknown, the SFH-recovery algorithm is able to accurately recover the input AMR, at the price of an increase of random errors in the SFR(t) by a factor of about 2.5. Experiments of SFH-recovery performed for varying distance modulus and reddening indicate that these parameters can be determined with (relative) accuracies of Delta(m-M)(0) similar to 0.02 mag and Delta E(B-V) similar to 0.01 mag, for each individual field over the LMC. The propagation of these errors in the SFR(t) implies systematic errors below 30%. This level of accuracy in the SFR(t) can reveal significant imprints in the dynamical evolution of this unique and nearby stellar system, as well as possible signatures of the past interaction between the MCs and the MW.

Modeling Two-Oscillator Circadian Systems Entrained by Two Environmental Cycles

Several experimental studies have altered the phase relationship between photic and non-photic environmental, 24 h cycles (zeitgebers) in order to assess their role in the synchronization of circadian rhythms. To assist in the interpretation of the complex activity patterns that emerge from these ""conflicting zeitgeber'' protocols, we present computer simulations of coupled circadian oscillators forced by two independent zeitgebers. This circadian system configuration was first employed by Pittendrigh and Bruce (1959), to model their studies of the light and temperature entrainment of the eclosion oscillator in Drosophila. Whereas most of the recent experiments have restricted conflicting zeitgeber experiments to two experimental conditions, by comparing circadian oscillator phases under two distinct phase relationships between zeitgebers (usually 0 and 12 h), Pittendrigh and Bruce compared eclosion phase under 12 distinct phase relationships, spanning the 24 h interval. Our simulations using non-linear differential equations replicated complex non-linear phenomena, such as ""phase jumps'' and sudden switches in zeitgeber preferences, which had previously been difficult to interpret. Our simulations reveal that these phenomena generally arise when inter-oscillator coupling is high in relation to the zeitgeber strength. Manipulations in the structural symmetry of the model indicated that these results can be expected to apply to a wide range of system configurations. Finally, our studies recommend the use of the complete protocol employed by Pittendrigh and Bruce, because different system configurations can generate similar results when a ""conflicting zeitgeber experiment'' incorporates only two phase relationships between zeitgebers.

Particle competition for complex network community detection

In many real situations, randomness is considered to be uncertainty or even confusion which impedes human beings from making a correct decision. Here we study the combined role of randomness and determinism in particle dynamics for complex network community detection. In the proposed model, particles walk in the network and compete with each other in such a way that each of them tries to possess as many nodes as possible. Moreover, we introduce a rule to adjust the level of randomness of particle walking in the network, and we have found that a portion of randomness can largely improve the community detection rate. Computer simulations show that the model has good community detection performance and at the same time presents low computational complexity. (C) 2008 American Institute of Physics.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

X-ray phase measurements as a probe of small structural changes in doped nonlinear optical crystals

X-ray multiple diffraction experiments with synchrotron radiation were carried out on pure and doped nonlinear optical crystals: NH(4)H(2)PO(4) and KH(2)PO(4) doped with Ni and Mn, respectively. Variations in the intensity profiles were observed from pure to doped samples, and these variations correlated with shifts in the structure factor phases, also known as triplet phases. This result demonstrates the potential of X-ray phase measurements to study doping in this type of single crystal. Different methodologies for probing structural changes were developed. Dynamical diffraction simulations and curve fitting procedures were also necessary for accurate phase determination. Structural changes causing the observed phase shifts are discussed.

Rapid formation of isoprene photo-oxidation products observed in Amazonia

Isoprene represents the single most important reactive hydrocarbon for atmospheric chemistry in the tropical atmosphere. It plays a central role in global and regional atmospheric chemistry and possible climate feedbacks. Photo-oxidation of primary hydrocarbons (e. g. isoprene) leads to the formation of oxygenated VOCs (OVOCs). The evolution of these intermediates affects the oxidative capacity of the atmosphere (by reacting with OH) and can contribute to secondary aerosol formation, a poorly understood process. An accurate and quantitative understanding of VOC oxidation processes is needed for model simulations of regional air quality and global climate. Based on field measurements conducted during the Amazonian Aerosol Characterization Experiment (AMAZE-08) we show that the production of certain OVOCs (e. g. hydroxyacetone) from isoprene photo-oxidation in the lower atmosphere is significantly underpredicted by standard chemistry schemes. Recently reported fast secondary production could explain 50% of the observed discrepancy with the remaining part possibly produced via a novel primary production channel, which has been proposed theoretically. The observations of OVOCs are also used to test a recently proposed HO(x) recycling mechanism via degradation of isoprene peroxy radicals. If generalized our observations suggest that prompt photochemical formation of OVOCs and other uncertainties in VOC oxidation schemes could result in uncertainties of modelled OH reactivity, potentially explaining a fraction of the missing OH sink over forests which has previously been largely attributed to a missing source of primary biogenic VOCs.

Multistability and Self-Similarity in the Parameter-Space of a Vibro-Impact System

The dynamics of a dissipative vibro-impact system called impact-pair is investigated. This system is similar to Fermi-Ulam accelerator model and consists of an oscillating one-dimensional box containing a point mass moving freely between successive inelastic collisions with the rigid walls of the box. In our numerical simulations, we observed multistable regimes, for which the corresponding basins of attraction present a quite complicated structure with smooth boundary. In addition, we characterize the system in a two-dimensional parameter space by using the largest Lyapunov exponents, identifying self-similar periodic sets. Copyright (C) 2009 Silvio L.T. de Souza et al.

Hidden vorticity in binary Bose-Einstein condensates

We consider a binary Bose-Einstein condensate (BEC) described by a system of two-dimensional (2D) Gross-Pitaevskii equations with the harmonic-oscillator trapping potential. The intraspecies interactions are attractive, while the interaction between the species may have either sign. The same model applies to the copropagation of bimodal beams in photonic-crystal fibers. We consider a family of trapped hidden-vorticity (HV) modes in the form of bound states of two components with opposite vorticities S(1,2) = +/- 1, the total angular momentum being zero. A challenging problem is the stability of the HV modes. By means of a linear-stability analysis and direct simulations, stability domains are identified in a relevant parameter plane. In direct simulations, stable HV modes feature robustness against large perturbations, while unstable ones split into fragments whose number is identical to the azimuthal index of the fastest growing perturbation eigenmode. Conditions allowing for the creation of the HV modes in the experiment are discussed too. For comparison, a similar but simpler problem is studied in an analytical form, viz., the modulational instability of an HV state in a one-dimensional (1D) system with periodic boundary conditions (this system models a counterflow in a binary BEC mixture loaded into a toroidal trap or a bimodal optical beam coupled into a cylindrical shell). We demonstrate that the stabilization of the 1D HV modes is impossible, which stresses the significance of the stabilization of the HV modes in the 2D setting.

Two-dimensional supersonic nonlinear Schrodinger flow past an extended obstacle

Supersonic flow of a superfluid past a slender impenetrable macroscopic obstacle is studied in the framework of the two-dimensional (2D) defocusing nonlinear Schroumldinger (NLS) equation. This problem is of fundamental importance as a dispersive analog of the corresponding classical gas-dynamics problem. Assuming the oncoming flow speed is sufficiently high, we asymptotically reduce the original boundary-value problem for a steady flow past a slender body to the one-dimensional dispersive piston problem described by the nonstationary NLS equation, in which the role of time is played by the stretched x coordinate and the piston motion curve is defined by the spatial body profile. Two steady oblique spatial dispersive shock waves (DSWs) spreading from the pointed ends of the body are generated in both half planes. These are described analytically by constructing appropriate exact solutions of the Whitham modulation equations for the front DSW and by using a generalized Bohr-Sommerfeld quantization rule for the oblique dark soliton fan in the rear DSW. We propose an extension of the traditional modulation description of DSWs to include the linear ""ship-wave"" pattern forming outside the nonlinear modulation region of the front DSW. Our analytic results are supported by direct 2D unsteady numerical simulations and are relevant to recent experiments on Bose-Einstein condensates freely expanding past obstacles.

Observation of double electron-positron pair production by gamma rays reexamined

An experiment was conducted to observe triple- and quadruple-escape peaks, at a photon energy equal to 6.128 MeV, in the spectra recorded with a high-purity Ge detector working in coincidence with six bismuth germanate detectors. The peak intensities may be explained having recourse to only the bremsstrahlung cascade process of consecutive electron-positron pair creation; i.e., the contribution of simultaneous double pair formation (and other cascade effects) is much smaller. The experimental peak areas are in reasonably good agreement with those predicted by Monte Carlo simulations done with the general-purpose radiation-tran sport code PENELOPE.

Critical behavior of the susceptible-infected-recovered model on a square lattice

By means of numerical simulations and epidemic analysis, the transition point of the stochastic asynchronous susceptible-infected-recovered model on a square lattice is found to be c(0)=0.176 500 5(10), where c is the probability a chosen infected site spontaneously recovers rather than tries to infect one neighbor. This point corresponds to an infection/recovery rate of lambda(c)=(1-c(0))/c(0)=4.665 71(3) and a net transmissibility of (1-c(0))/(1+3c(0))=0.538 410(2), which falls between the rigorous bounds of the site and bond thresholds. The critical behavior of the model is consistent with the two-dimensional percolation universality class, but local growth probabilities differ from those of dynamic percolation cluster growth, as is demonstrated explicitly.

Diffusion anomaly and dynamic transitions in the Bell-Lavis water model

In this paper we investigate the dynamic properties of the minimal Bell-Lavis (BL) water model and their relation to the thermodynamic anomalies. The BL model is defined on a triangular lattice in which water molecules are represented by particles with three symmetric bonding arms interacting through van der Waals and hydrogen bonds. We have studied the model diffusivity in different regions of the phase diagram through Monte Carlo simulations. Our results show that the model displays a region of anomalous diffusion which lies inside the region of anomalous density, englobed by the line of temperatures of maximum density. Further, we have found that the diffusivity undergoes a dynamic transition which may be classified as fragile-to-strong transition at the critical line only at low pressures. At higher densities, no dynamic transition is seen on crossing the critical line. Thus evidence from this study is that relation of dynamic transitions to criticality may be discarded. (C) 2010 American Institute of Physics. [doi:10.1063/1.3479001]

Hydrogen bond interactions between acetone and supercritical water

Hydrogen bond interactions between acetone and supercritical water are investigated using a combined and sequential Monte Carlo/quantum mechanics (S-MC/QM) approach. Simulation results show a dominant presence of con. gurations with one hydrogen bond for different supercritical states, indicating that this specific interaction plays an important role on the solvation properties of acetone in supercritical water. Using QM MP2/aug-cc-pVDZ the calculated average interaction energy reveals that the hydrogen-bonded acetone-water complex is energetically more stable under supercritical conditions than ambient conditions and its stability is little affected by variations of temperature and/or pressure. All average results reported here are statistically converged.

Electron collisions with the CH(2)O-H(2)O complex

We report cross sections for elastic collisions of low-energy electrons with the CH(2)O-H(2)O complex. We employed the Schwinger multichannel method with pseudopotentials in the static-exchange and in the static-exchange-polarization approximations for energies from 0.1 to 20 eV. We considered four different hydrogen-bonded structures for the complex that were generated by classical Monte Carlo simulations. Our aim is to investigate the effect of the water molecule on the pi* shape resonance of formaldehyde. Previous studies reported a pi* shape resonance for CH(2)O at around 1 eV. The resonance positions of the complexes appear at lower energies in all cases due to the mutual polarization between the two molecules. This indicates that the presence of water may favor dissociation by electron impact and may lead to an important effect on strand breaking in wet DNA by electron impact.