Thermal instability and the feedback regulation of hot haloes in clusters, groups and galaxies

We present global multidimensional numerical simulations of the plasma that pervades the dark matter haloes of clusters, groups and massive galaxies (the intracluster medium; ICM). Observations of clusters and groups imply that such haloes are roughly in global thermal equilibrium, with heating balancing cooling when averaged over sufficiently long time- and length-scales; the ICM is, however, very likely to be locally thermally unstable. Using simple observationally motivated heating prescriptions, we show that local thermal instability (TI) can produce a multiphase medium with similar to 104 K cold filaments condensing out of the hot ICM only when the ratio of the TI time-scale in the hot plasma (tTI) to the free-fall time-scale (tff) satisfies tTI/tff? 10. This criterion quantitatively explains why cold gas and star formation are preferentially observed in low-entropy clusters and groups. In addition, the interplay among heating, cooling and TI reduces the net cooling rate and the mass accretion rate at small radii by factors of similar to 100 relative to cooling-flow models. This dramatic reduction is in line with observations. The feedback efficiency required to prevent a cooling flow is similar to 10-3 for clusters and decreases for lower mass haloes; supernova heating may be energetically sufficient to balance cooling in galactic haloes. We further argue that the ICM self-adjusts so that tTI/tff? 10 at all radii. When this criterion is not satisfied, cold filaments condense out of the hot phase and reduce the density of the ICM. These cold filaments can power the black hole and/or stellar feedback required for global thermal balance, which drives tTI/tff? 10. In comparison to clusters, groups have central cores with lower densities and larger radii. This can account for the deviations from self-similarity in the X-ray luminositytemperature () relation. The high-velocity clouds observed in the Galactic halo can be due to local TI producing multiphase gas close to the virial radius if the density of the hot plasma in the Galactic halo is >rsim 10-5 cm-3 at large radii.

Energy-dependent orbital modulation of X-rays and constraints on emission of the jet in Cyg X-3

We study the orbital modulation of X-rays from Cyg X-3, using data from Swift, INTEGRAL and RXTE. Using the wealth of data presently available and an improved averaging method, we obtain energy-dependent folded and averaged light curves with unprecedented accuracy. We find that above similar to 5?keV the modulation depth decreases with increasing energy, which is consistent with the modulation being caused by both boundfree absorption and Compton scattering in the stellar wind of the donor, with minima corresponding to the highest optical depth, which occurs around the superior conjunction. We find a decrease of the depth below similar to 3?keV, which appears to be due to re-emission of the absorbed continuum by the wind in soft X-ray lines. Based on the shape of the folded light curves, any X-ray contribution from the jet in Cyg X-3, which emits ?-rays detected at energies >0.1?GeV in the soft spectral states, is found to be minor up to similar to 100?keV. This implies the presence of a rather sharp low-energy break in the jet MeV-range spectrum. We also calculate phase-resolved RXTE X-ray spectra and show that the difference between the spectra corresponding to phases around superior and inferior conjunctions can indeed be accounted for by the combined effect of boundfree absorption in an ionized medium and Compton scattering.

Adsorption on Edge-Functionalized Bilayer Graphene Nanoribbons: Assessing the Role of Functional Groups in Methane Uptake

A combination of ab initio and classical Monte Carlo simulations is used to investigate the effects of functional groups on methane binding. Using Moller-Plesset (MP2) calculations, we obtain the binding energies for benzene functionalized with NH2, OH, CH3, COOH, and H2PO3 and identify the methane binding sites. In all cases, the preferred binding sites are located above the benzene plane in the vicinity of the benzene carbon atom attached to the functional group. Functional groups enhance methane binding relative to benzene (-6.39 kJ/mol), with the largest enhancement observed for H2PO3 (-8.37 kJ/mol) followed by COOH and CH3 (-7.77 kJ/mol). Adsorption isotherms are obtained for edge-functionalized bilayer graphene nanoribbons using grand canonical Monte Carlo simulations with a five-site methane model. Adsorbed excess and heats of adsorption for pressures up to 40 bar and 298 K are obtained with functional group concentrations ranging from 3.125 to 6.25 mol 96 for graphene edges functionalized with OH, NH2, and COOH. The functional groups are found to act as preferred adsorption sites, and in the case of COOH the local methane density in the vicinity of the functional group is found to exceed that of bare graphene. The largest enhancement of 44.5% in the methane excess adsorbed is observed for COOH-functionalized nanoribbons when compared to H terminated ribbons. The corresponding enhancements for OH- and NH2-functionalized ribbons are 10.5% and 3.7%, respectively. The excess adsorption across functional groups reflects the trends observed in the binding energies from MP2 calculations. Our study reveals that specific site functionalization can have a significant effect on the local adsorption characteristics and can be used as a design strategy to tailor materials with enhanced methane storage capacity.

Cavitation in materials with distributed weak zones: Implications on the origin of brittle fracture in metallic glasses

Several experimental studies have shown that fracture surfaces in brittle metallic glasses (MGs) generally exhibit nanoscale corrugations which may be attributed to the nucleation and coalescence of nanovoids during crack propagation. Recent atomistic simulations suggest that this phenomenon is due to large spatial fluctuations in material properties in a brittle MG, which leads to void nucleation in regions of low atomic density and then catastrophic fracture through void coalescence. To explain this behavior, we propose a model of a heterogeneous solid containing a distribution of weak zones to represent a brittle MG. Plane strain continuum finite element analysis of cavitation in such an elastic-plastic solid is performed with the weak zones idealized as periodically distributed regions having lower yield strength than the background material. It is found that the presence of weak zones can significantly reduce the critical hydrostatic stress for the onset of cavitation which is controlled uniquely by the local yield properties of these zones. Also, the presence of weak zones diminishes the sensitivity of the cavitation stress to the volume fraction of a preexisting void. These results provide plausible explanations for the observations reported in recent atomistic simulations of brittle MGs. An analytical solution for a composite, incompressible elastic-plastic solid with a weak inner core is used to investigate the effect of volume fraction and yield strength of the core on the nature of cavitation bifurcation. It is shown that snap-cavitation may occur, giving rise to sudden formation of voids with finite size, which does not happen in a homogeneous plastic solid. (c) 2012 Elsevier Ltd. All rights reserved.

Monte Carlo simulation of light transport in tissue for optimizing light delivery in photoacoustic imaging of the sentinel lymph node

Noninvasive or minimally invasive identification of sentinel lymph node (SLN) is essential to reduce the surgical effects of SLN biopsy. Photoacoustic (PA) imaging of SLN in animal models has shown its promise for clinical use in the future. Here, we present a Monte Carlo simulation for light transport in the SLN for various light delivery configurations with a clinical ultrasound probe. Our simulation assumes a realistic tissue layer model and also can handle the transmission/reflectance at SLN-tissue boundary due to the mismatch of refractive index. Various light incidence angles show that for deeply situated SLNs the maximum absorption of light in the SLN is for normal incidence. We also show that if a part of the diffused reflected photons is reflected back into the skin using a reflector, the absorption of light in the SLN can be increased significantly to enhance the PA signal. (C) 2013 Society of Photo-Optical Instrumentation Engineers (SPIE)

Uplift Resistance of Long Pipelines in the Presence of Seismic Forces

The vertical uplift resistance of long pipes buried in sands and subjected to pseudostatic seismic forces has been computed by using the lower-bound theorem of the limit analysis in conjunction with finite elements and nonlinear optimization. The soil mass is assumed to follow the Mohr-Coulomb failure criterion and an associated flow rule. The failure load is expressed in the form of a nondimensional uplift factor F-gamma. The variation of F-gamma is plotted as a function of the embedment ratio of the pipe, horizontal seismic acceleration coefficient (k(h)), and soil friction angle (phi). The magnitude of F-gamma is found to decrease continuously with an increase in the horizontal seismic acceleration coefficient. The reduction in the uplift resistance becomes quite significant, especially for greater values of embedment ratios and lower values of friction angle. The predicted uplift resistance was found to compare well with the existing results reported from the literature. (C) 2014 American Society of Civil Engineers.

Quantitative estimation of density variation in high-speed flows through inversion of the measured wavefront distortion

A simple method employing an optical probe is presented to measure density variations in a hypersonic flow obstructed by a test model in a typical shock tunnel. The probe has a plane light wave trans-illuminating the flow and casting a shadow of a random dot pattern. Local slopes of the distorted wavefront are obtained from shifts of the dots in the pattern. Local shifts in the dots are accurately measured by cross-correlating local shifted shadows with the corresponding unshifted originals. The measured slopes are suitably unwrapped by using a discrete cosine transform based phase unwrapping procedure and also through iterative procedures. The unwrapped phase information is used in an iterative scheme for a full quantitative recovery of density distribution in the shock around the model through refraction tomographic inversion. Hypersonic flow field parameters around a missile shaped body at a free-stream Mach number of 5.8 measured using this technique are compared with the numerically estimated values. (C) 2014 Society of Photo-Optical Instrumentation Engineers (SPIE)

The cold mode: a phenomenological model for the evolution of density perturbations in the intracluster medium

Cool cluster cores are in global thermal equilibrium but are locally thermally unstable. We study a non-linear phenomenological model for the evolution of density perturbations in the intracluster medium (ICM) due to local thermal instability and gravity. We have analysed and extended a model for the evolution of an overdense blob in the ICM. We find two regimes in which the overdense blobs can cool to thermally stable low temperatures. One for large t(cool)/t(ff) (t(cool) is the cooling time and t(ff) is the free-fall time), where a large initial overdensity is required for thermal runaway to occur; this is the regime which was previously analysed in detail. We discover a second regime for t(cool)/t(ff) less than or similar to 1 (in agreement with Cartesian simulations of local thermal instability in an external gravitational field), where runaway cooling happens for arbitrarily small amplitudes. Numerical simulations have shown that cold gas condenses out more easily in a spherical geometry. We extend the analysis to include geometrical compression in weakly stratified atmospheres such as the ICM. With a single parameter, analogous to the mixing length, we are able to reproduce the results from numerical simulations; namely, small density perturbations lead to the condensation of extended cold filaments only if t(cool)/t(ff) less than or similar to 10.

Stability of Filaments in Star-Forming Clouds and the Formation of Prestellar Cores in Them

The exact process(es) that generate(s) dense filaments which then form prestellar cores within them is unclear. Here we study the formation of a dense filament using a relatively simple set-up of a pressure-confined, uniform-density cylinder. We examine if its propensity to form a dense filament and further, to the formation of prestellar cores along this filament, bears on the gravitational state of the initial volume of gas. We report a radial collapse leading to the formation of a dense filamentary cloud is likely when the initial volume of gas is at least critically stable (characterised by the approximate equality between the mass line-density for this volume and its maximum value). Though self-gravitating, this volume of gas, however, is not seen to be in free-fall. This post-collapse filament then fragments along its length due to the growth of a Jeans-like instability to form prestellar cores. We suggest dense filaments in typical star-forming clouds classified as gravitationally super-critical under the assumption of: (i) isothermality when in fact, they are not, and (ii) extended radial profiles as against pressure-truncated, that significantly over-estimates their mass line-density, are unlikely to experience gravitational free-fall. The radial density and temperature profile derived for this post-collapse filament is consistent with that deduced for typical filamentary clouds mapped in recent surveys of nearby star-forming regions.

Evolution of social behaviour in the primitively eusocial wasp Ropalidia marginata: do we need to look beyond kin selection?

Ropalidia marginata is a primitively eusocial wasp widely distributed in peninsular India. Although solitary females found a small proportion of nests, the vast majority of new nests are founded by small groups of females. In suchmultiple foundress nests, a single dominant female functions as the queen and lays eggs, while the rest function as sterile workers and care for the queen's brood. Previous attempts to understand the evolution of social behaviour and altruism in this species have employed inclusive fitness theory (kin selection) as a guiding framework. Although inclusive fitness theory is quite successful in explaining the high propensity of the wasps to found nests in groups, several features of their social organization suggest that forces other than kin selection may also have played a significant role in the evolution of this species. These features include lowering of genetic relatedness owing to polyandry and serial polygyny, nest foundation by unrelated individuals, acceptance of young non-nest-mates, a combination of well-developed nest-mate recognition and lack of intra-colony kin recognition, a combination of meek and docile queens and a decentralized self-organized work force, long reproductive queues with cryptic heir designates and conflict-free queen succession, all resulting in extreme intra-colony cooperation and inter-colony conflict.