Multidimensional modelling of X-ray spectra for AGN accretion disc outflows - II

Highly ionized fast accretion disc winds have been suggested as an explanation for a variety of observed absorption and emission features in the X-ray spectra of active galactic nuclei. Simple estimates have suggested that these flows may be massive enough to carry away a significant fraction of the accretion energy and could be involved in creating the link between supermassive black holes and their host galaxies. However, testing these hypotheses, and quantifying the outflow signatures, requires high-quality theoretical spectra for comparison with observations. Here, we describe extensions of our Monte Carlo radiative transfer code that allow us to generate realistic theoretical spectra for a much wider variety of disc wind models than that was possible in our previous work. In particular, we have expanded the range of atomic physics simulated by the code so that L- and M-shell ions can now be included. We have also substantially improved our treatment of both ionization and radiative heating such that we are now able to compute spectra for outflows containing far more diverse plasma conditions. We present example calculations that illustrate the variety of spectral features predicted by parametrized outflow models and demonstrate their applicability to the interpretation of data by comparison with observations of the bright quasar PG1211+143. We find that the major features in the observed 2-10 keV spectrum of this object can be well reproduced by our spectra, confirming that it likely hosts a massive outflow. © 2010 The Authors. Journal compilation © 2010 RAS.

Harmonic generation by noble-gas atoms in the near-IR regime using ab initio time-dependent R-matrix theory

We demonstrate the capability of ab initio time-dependent R-matrix theory to obtain accurate harmonic generation spectra of noble-gas atoms at near-IR wavelengths between 1200 and 1800 nm and peak intensities up to 1.8 × 10^(14) W/cm^(2). To accommodate the excursion length of the ejected electron, we use an angular-momentum expansion up to Lmax=279. The harmonic spectra show evidence of atomic structure through the presence of a Cooper minimum in harmonic generation for Kr, and of multielectron interaction through the giant resonance for Xe. The theoretical spectra agree well with those obtained experimentally.

Analysis of the Second Harmonic Effect on Power Amplifier Intermodulation Products

A theoretical analysis is reported in this paper to investigate the effect that a second harmonic signal which might be present at an amplifier’s input has on generating additional intermodulation products, particularly the third-order intermodulation (IM3) products. The analysis shows that the amplitude of an extra generated IM3 component is equal to the product of the fundamental amplitude, the second harmonic amplitude, and the second order Taylor series coefficient. The effect of the second order harmonic on the IM3 is examined through a simulated example of a 2.22-GHz 10-W Class-EF amplifier whereby the IM3 levels have been reduced by 2-3 dB after employing a second harmonic termination stub at the input.

Atmospheric parameters and rotational velocities for a sample of Galactic B-type supergiants

High-resolution optical spectra of 57 Galactic B-type supergiant stars have been analysed to determine their rotational and macroturbulent velocities. In addition, their atmospheric parameters (effective temperature, surface gravity and microturbulent velocity) and surface nitrogen abundances have been estimated using a non-local thermodynamic equilibrium grid of model atmospheres. Comparisons of the projected rotational velocities have been made with the predictions of stellar evolutionary models and in general good agreement was found. However, for a small number of targets, their observed rotational velocities were significantly larger than predicted, although their nitrogen abundances were consistent with the rest of the sample. We conclude that binarity may have played a role in generating their large rotational velocities. No correlation was found between nitrogen abundances and the current projected rotational velocities. However, a correlation was found with the inferred projected rotational velocities of the main-sequence precursors of our supergiant sample. This correlation is again in agreement with the predictions of single star evolutionary models that incorporate rotational mixing. The origin of the macroturbulence and microturbulent velocity fields is discussed and our results support previous theoretical studies that link the former to subphotospheric convection and the latter to non-radial gravity mode oscillations. In addition, we have attempted to identify differential rotation in our most rapidly rotating targets.

Combined R-matrix eigenstate basis set and finite-difference propagation method for the time-dependent Schrodinger equation: The one-electron case

In this work we present the theoretical framework for the solution of the time-dependent Schrödinger equation (TDSE) of atomic and molecular systems under strong electromagnetic fields with the configuration space of the electron’s coordinates separated over two regions; that is, regions I and II. In region I the solution of the TDSE is obtained by an R-matrix basis set representation of the time-dependent wave function. In region II a grid representation of the wave function is considered and propagation in space and time is obtained through the finite-difference method. With this, a combination of basis set and grid methods is put forward for tackling multiregion time-dependent problems. In both regions, a high-order explicit scheme is employed for the time propagation. While, in a purely hydrogenic system no approximation is involved due to this separation, in multielectron systems the validity and the usefulness of the present method relies on the basic assumption of R-matrix theory, namely, that beyond a certain distance (encompassing region I) a single ejected electron is distinguishable from the other electrons of the multielectron system and evolves there (region II) effectively as a one-electron system. The method is developed in detail for single active electron systems and applied to the exemplar case of the hydrogen atom in an intense laser field.

Breit interaction in dielectronic recombination of hydrogenlike uranium

Absolute rate coefficients for dielectronic recombination (DR) of H-like U91+ ions have been measured. The electron-ion merged-beam technique at a heavy-ion storage ring was employed using a stochastically cooled ion beam. Thereby, the previously accessible electron-ion collision energies could be greatly extended to the range 63-90 keV. High-resolution DR spectra were measured covering all KLL and KLM resonances. For the resonance strengths, excellent agreement between relativistic theory and experiment is found only if the Breit contribution to the electron-electron interaction is included in the calculations. For the KL1/2L1/2 and KL1/2M1/2 groups the Breit contribution amounts to 44% of their total resonance strengths.

Analysis of luminescent samples using subtracted shifted Raman spectroscopy

A novel method of obtaining high-quality Raman spectra of luminescent samples was tested using cyclohexane solutions which had been treated with a fluorescent dye. The method involves removing the fixed pattern irregularity found in the spectra taken with CCD detectors by subtracting spectra taken at several different, closely spaced spectrometer positions. It is conceptually similar to SERDS (shifted excitation Raman difference spectroscopy) but has the distinct experimental advantage that it does not require a tunable laser source. The subtracted spectra obtained as the raw data are converted into a more recognisable and conventional form by iterative fitting of appropriate double Lorentzian functions whose peak parameters are then used to 'reconstruct' a conventional representation of the spectrum. Importantly, it is shown that the degree of uncertainty in the resultant 'reconstructed' spectra can be gauged reliably by comparing reconstructed spectra obtained at two different spectrometer shifts (delta and 2 delta), The method was illustrated and validated using a solvent (cyclohexane) the spectrum of which is well known and which contains both regions with complex overlapping bands and regions with isolated bands, Possible sources of error are discussed and it is shown that, provided the degree of uncertainty in the data is correctly characterised, it is completely valid to draw conclusions about the spectra of the sample on the basis of the reconstructed data. The acronym SSRS (subtracted shifted Raman spectroscopy; pronounced scissors) is proposed for this method, to distinguish it from the SERDS technique.

Metallicities of 4 blue supergiants near the Galactic centre

High resolution optical spectra of four blue supergiants (HD148422 B0.5 Ib; HD178487 B0.5 Ib; HD179407 B1 Ib, HD163522 B1 Ib) which lie within 4.5 kpc of the Galactic centre are presented. Careful differential LTE model atmosphere analyses are used to quantify the differences in photospheric metal abundances between these stars and MK spectral standards in the solar neighborhood. A detailed non-LTE model atmosphere analysis of one star (HD163522) confirms that the LTE differential abundances should be reliable, provided we use a comparison star with similar atmospheric parameters.

Label-free detection of nanomolar unmodified single- and double-stranded dna by using surface-enhanced raman spectroscopy on Ag and Au colloids

Unlabelled single- and double-stranded DNA (ssDNA and dsDNA, respectively) has been detected at concentrations =10-9?M by surface-enhanced Raman spectroscopy. Under appropriate conditions the sequences spontaneously adsorbed to the surface of both Ag and Au colloids through their nucleobases; this allowed highly reproducible spectra with good signal-to-noise ratios to be recorded on completely unmodified samples. This eliminated the need to promote absorption by introducing external linkers, such as thiols. The spectra of model ssDNA sequences contained bands of all the bases present and showed systematic changes when the overall base composition was altered. Initial tests also showed that small but reproducible changes could be detected between oligonucleotides with the same bases arranged in a different order. The spectra of five ssDNA sequences that correspond to different strains of the Escherichia coli bacterium were found to be sufficiently composition-dependent so that they could be differentiated without the need for any advanced multivariate data analysis techniques.

Proton ordering energetics in ice phases

Results from first-principles calculations on the subtle energetics of proton ordering in ice phases are shown only to depend on the electrostatic components of the total energy. Proton ordered ice phases can therefore be predicted using electronic structure methods or a tailored potential model. However, analysis of the electron density reveals that high order multipole components, up to hexadecapole, are needed to adequately capture total energy differences between proton ordered and disordered phases. This suggests that current potential models may be unable to reproduce the position of proton ordered ice phases in the phase diagram without extensions to describe high order electrostatics. (c) 2006 Elsevier B.V. All rights reserved.

Assessing the Nonequilibrium Thermodynamics in a Quenched Quantum Many-Body System via Single Projective Measurements

We analyze the nature of the statistics of the work done on or by a quantum many-body system brought out of equilibrium. We show that, for the sudden quench and for an initial state that commutes with the initial Hamiltonian, it is possible to retrieve the whole nonequilibrium thermodynamics via single projective measurements of observables. We highlight, in a physically clear way, the qualitative implications for the statistics of work coming from considering processes described by operators that either commute or do not commute with the unperturbed Hamiltonian of a given system. We consider a quantum many-body system and derive an expression that allows us to give a physical interpretation, for a thermal initial state, to all of the cumulants of the work in the case of quenched operators commuting with the unperturbed Hamiltonian. In the commuting case, the observables that we need to measure have an intuitive physical meaning. Conversely, in the noncommuting case, we show that, although it is possible to operate fully within the single-measurement framework irrespectively of the size of the quench, some difficulties are faced in providing a clear-cut physical interpretation to the cumulants. This circumstance makes the study of the physics of the system nontrivial and highlights the nonintuitive phenomenology of the emergence of thermodynamics from the fully quantum microscopic description. We illustrate our ideas with the example of the Ising model in a transverse field showing the interesting behavior of the high-order statistical moments of the work distribution for a generic thermal state and linking them to the critical nature of the model itself.

Reducing Parametric Uncertainty in Limit Cycle Oscillations

The assimilation of discrete higher fidelity data points with model predictions can be used to achieve a reduction in the uncertainty of the model input parameters which generate accurate predictions. The problem investigated here involves the prediction of limit-cycle oscillations using a High-Dimensional Harmonic Balance method (HDHB). The efficiency of the HDHB method is exploited to enable calibration of structural input parameters using a Bayesian inference technique. Markov-chain Monte Carlo is employed to sample the posterior distributions. Parameter estimation is carried out on both a pitch/plunge aerofoil and Goland wing configuration. In both cases significant refinement was achieved in the distribution of possible structural parameters allowing better predictions of their
true deterministic values.

A Continuous Time Markov Chain Based Sequential Analytical Approach for Composite Power System Reliability Assessment

This paper proposes a continuous time Markov chain (CTMC) based sequential analytical approach for composite generation and transmission systems reliability assessment. The basic idea is to construct a CTMC model for the composite system. Based on this model, sequential analyses are performed. Various kinds of reliability indices can be obtained, including expectation, variance, frequency, duration and probability distribution. In order to reduce the dimension of the state space, traditional CTMC modeling approach is modified by merging all high order contingencies into a single state, which can be calculated by Monte Carlo simulation (MCS). Then a state mergence technique is developed to integrate all normal states to further reduce the dimension of the CTMC model. Moreover, a time discretization method is presented for the CTMC model calculation. Case studies are performed on the RBTS and a modified IEEE 300-bus test system. The results indicate that sequential reliability assessment can be performed by the proposed approach. Comparing with the traditional sequential Monte Carlo simulation method, the proposed method is more efficient, especially in small scale or very reliable power systems.

Testing quasar unification: radiative transfer in clumpy winds

Various unification schemes interpret the complex phenomenology of quasars and luminous active galactic nuclei (AGN) in terms of a simple picture involving a central black hole, an accretion disc and an associated outflow. Here, we continue our tests of this paradigm by comparing quasar spectra to synthetic spectra of biconical disc wind models, produced with our state-of-the-art Monte Carlo radiative transfer code. Previously, we have shown that we could produce synthetic spectra resembling those of observed broad absorption line (BAL) quasars, but only if the X-ray luminosity was limited to 1043 erg s-1. Here, we introduce a simple treatment of clumping, and find that a filling factor of ˜0.01 moderates the ionization state sufficiently for BAL features to form in the rest-frame UV at more realistic X-ray luminosities. Our fiducial model shows good agreement with AGN X-ray properties and the wind produces strong line emission in, e.g., Lyα and C IV 1550 Å at low inclinations. At high inclinations, the spectra possess prominent LoBAL features. Despite these successes, we cannot reproduce all emission lines seen in quasar spectra with the correct equivalent-width ratios, and we find an angular dependence of emission line equivalent width despite the similarities in the observed emission line properties of BAL and non-BAL quasars. Overall, our work suggests that biconical winds can reproduce much of the qualitative behaviour expected from a unified model, but we cannot yet provide quantitative matches with quasar properties at all viewing angles. Whether disc winds can successfully unify quasars is therefore still an open question.

How curvature-generating proteins build scaffolds on membrane nanotubes

Bin/Amphiphysin/Rvs (BAR) domain proteins control the curvature of lipid membranes in endocytosis, trafficking, cell motility, the formation of complex sub-cellular structures, and many other cellular phenomena. They form three-dimensional assemblies, which act as molecular scaffolds to reshape the membrane and alter its mechanical properties. It is unknown, however, how a protein scaffold forms and how BAR domains interact in these assemblies at protein densities relevant for a cell. In this work, we employ various experimental, theoretical and simulation approaches to explore how BAR proteins organize to form a scaffold on a membrane nanotube. By combining quantitative microscopy with analytical modeling, we demonstrate that a highly curving BAR protein endophilin nucleates its scaffolds at the ends of a membrane tube, contrary to a weaker curving protein centaurin, which binds evenly along the tube’s length. Our work implies that the nature of local protein-membrane interactions can affect the specific localization of proteins on membrane-remodeling sites. Furthermore, we show that amphipathic helices are dispensable in forming protein scaffolds. Finally, we explore a possible molecular structure of a BAR-domain scaffold using coarse-grained molecular dynamics simulations. Together with fluorescence microscopy, the simulations show that proteins need only to cover 30–40% of a tube’s surface to form a rigid assembly. Our work provides mechanical and structural insights into the way BAR proteins may sculpt the membrane as a high-order cooperative assembly in important biological processes.