Transport coefficients for the shear dynamo problem at small Reynolds numbers

We build on the formulation developed in S. Sridhar and N. K. Singh J. Fluid Mech. 664, 265 (2010)] and present a theory of the shear dynamo problem for small magnetic and fluid Reynolds numbers, but for arbitrary values of the shear parameter. Specializing to the case of a mean magnetic field that is slowly varying in time, explicit expressions for the transport coefficients alpha(il) and eta(iml) are derived. We prove that when the velocity field is nonhelical, the transport coefficient alpha(il) vanishes. We then consider forced, stochastic dynamics for the incompressible velocity field at low Reynolds number. An exact, explicit solution for the velocity field is derived, and the velocity spectrum tensor is calculated in terms of the Galilean-invariant forcing statistics. We consider forcing statistics that are nonhelical, isotropic, and delta correlated in time, and specialize to the case when the mean field is a function only of the spatial coordinate X-3 and time tau; this reduction is necessary for comparison with the numerical experiments of A. Brandenburg, K. H. Radler, M. Rheinhardt, and P. J. Kapyla Astrophys. J. 676, 740 (2008)]. Explicit expressions are derived for all four components of the magnetic diffusivity tensor eta(ij) (tau). These are used to prove that the shear-current effect cannot be responsible for dynamo action at small Re and Rm, but for all values of the shear parameter.

Network approach for capturing ligand-induced subtle global changes in protein structures

Ligand-induced conformational changes in proteins are of immense functional relevance. It is a major challenge to elucidate the network of amino acids that are responsible for the percolation of ligand-induced conformational changes to distal regions in the protein from a global perspective. Functionally important subtle conformational changes (at the level of side-chain noncovalent interactions) upon ligand binding or as a result of environmental variations are also elusive in conventional studies such as those using root-mean-square deviations (r.m.s.d.s). In this article, the network representation of protein structures and their analyses provides an efficient tool to capture these variations (both drastic and subtle) in atomistic detail in a global milieu. A generalized graph theoretical metric, using network parameters such as cliques and/or communities, is used to determine similarities or differences between structures in a rigorous manner. The ligand-induced global rewiring in the protein structures is also quantified in terms of network parameters. Thus, a judicious use of graph theory in the context of protein structures can provide meaningful insights into global structural reorganizations upon perturbation and can also be helpful for rigorous structural comparison. Data sets for the present study include high-resolution crystal structures of serine proteases from the S1A family and are probed to quantify the ligand-induced subtle structural variations.

The role of Compton and Raman scattering in the quasar continuum

There are three ways in which an electromagnetic wave can undergo scattering in a plasma: (i) when the scattering of radiation occurs by a single electron, it is called Compton Scattering (CS); (ii) if it occurs by a longitudinal electron plasma mode, it is called Stimulated Raman Scattering (SRS), and (iii) if it occurs by a highly damped electron plasma mode, it is called Stimulated Compton Scattering (SCS). The non-thermal continuum of quasars is believed to be produced through the combined action of synchrotron and inverse Compton processes, which are essentially single-particle processes. Here, we investigate the role of SRS and SCS in the generation of continuum radiation from these compact objects. It is shown as an example that the complete spectrum of 3C 273 can be reproduced by suitably combining SCS and SRS. The differential contributions of SCS and SRS under different values of the plasma parameters are also calculated.

Potential and economics of forestry options for carbon sequestration in India

There is a need to understand the carbon (C) sequestration potential of the forestry option and its financial implications for each country.In India the C emissions from deforestation are estimated to be nearly offset by C sequestration in forests under succession and tree plantations. India has nearly succeeded in stabilizing the area under forests and has adequate forest conservation strategies. Biomass demands for softwood, hardwood and firewood are estimated to double or treble by the year 2020. A set of forestry options were developed to meet the projected biomass needs, and keeping in mind the features of land categories available, three scenarios were developed: potential; demand-driven; and programme-driven scenarios. Adoption of the demand-driven scenario, targeted at meeting the projected biomass needs, is estimated to sequester 78 Mt of C annually after accounting for all emissions resulting from clearfelling and end use of biomass. The demand-driven scenario is estimated to offset 50% of national C emission at 1990 level. The cost per t of C sequestered for forestry options is lower than the energy options considered. The annual investment required for implementing the demand-driven scenario is estimated to be US$ 2.1 billion for six years and is shown to be feasible. Among forestry options, the ranking based on investment cost per t of C sequestered from least cost to highest cost is; natural regeneration-agro-forestry-enhanced natural regeneration (< US$ 2.5/t C)-timber-community-softwood forestry (US$ 3.3 to 7.3 per t of C).

Stability of spatially developing boundary layers in pressure gradients

A new formulation of the stability of boundary-layer flows in pressure gradients is presented, taking into account the spatial development of the flow and utilizing a special coordinate transformation. The formulation assumes that disturbance wavelength and eigenfunction vary downstream no more rapidly than the boundary-layer thickness, and includes all terms nominally of order R(-1) in the boundary-layer Reynolds number R. In Blasius flow, the present approach is consistent with that of Bertolotti et al. (1992) to O(R(-1)) but simpler (i.e. has fewer terms), and may best be seen as providing a parametric differential equation which can be solved without having to march in space. The computed neutral boundaries depend strongly on distance from the surface, but the one corresponding to the inner maximum of the streamwise velocity perturbation happens to be close to the parallel flow (Orr-Sommerfeld) boundary. For this quantity, solutions for the Falkner-Skan flows show the effects of spatial growth to be striking only in the presence of strong adverse pressure gradients. As a rational analysis to O(R(-1)) demands inclusion of higher-order corrections on the mean flow, an illustrative calculation of one such correction, due to the displacement effect of the boundary layer, is made, and shown to have a significant destabilizing influence on the stability boundary in strong adverse pressure gradients. The effect of non-parallelism on the growth of relatively high frequencies can be significant at low Reynolds numbers, but is marginal in other cases. As an extension of the present approach, a method of dealing with non-similar flows is also presented and illustrated. However, inherent in the transformation underlying the present approach is a lower-order non-parallel theory, which is obtained by dropping all terms of nominal order R(-1) except those required for obtaining the lowest-order solution in the critical and wall layers. It is shown that a reduced Orr-Sommerfeld equation (in transformed coordinates) already contains the major effects of non-parallelism.

Stability of the viscous flow of a fluid through a flexible tube

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (rho VR/eta), the ratio of the viscosities of the wall and fluid eta(r) = (eta(s)/eta), the ratio of radii H and the dimensionless velocity Gamma = (rho V-2/G)(1/2). Here rho is the density of the fluid, G is the coefficient of elasticity of the wall and V is the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter epsilon = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate s((0)), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctuations due to the Reynolds stress. There is an O(epsilon(1/2)) correction to the growth rate, s((1)), due to the presence of a wall layer of thickness epsilon(1/2)R where the viscous stresses are O(epsilon(1/2)) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Gamma and wavenumber k where s((1)) = 0. At these points, the wall layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(epsilon) correction to the growth rate s((2)) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s((2)) increases proportional to (H-1)(-2) for (H-1) much less than 1 (thickness of wall much less than the tube radius), and decreases proportional to H-4 for H much greater than 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube

The solar dynamo with meridional circulation

We show that meridional circulation can have a profound influence on dynamo models for the solar cycle. Motivated by the observed tilt angles of sunspot groups we assume that the generation of the poloidal field takes place near the surface, while a shear layer of radial differential rotation produces the toroidal field at the bottom of the convection zone. Both layers are coupled by a circulation with a poleward directed flow in the upper part and an equatorward flow in the deep layers of the convection zone. The circulation forces the toroidal field belts (which are responsible for the surface activity) to move equatorward. This leads to butterfly diagrams in qualitative agreement with the observations, even if the dynamo wave would propagate poleward in the absence of circulation. This result opens the possibility to construct models for the solar cycle which are based on observational data (tilt angles, differential rotation, and meridional circulation).

Spectroscopic Observation of Oscillations in the Corona During the Total Solar Eclipse of 22 July 2009

We performed high resolution spectroscopy of the solar corona during the total solar eclipse of 22 July 2009 in two emission lines: the green line at 5303 due to Fe xiv and the red line at 6374 due to Fe x, simultaneously from Anji (latitude 30A degrees 28.1' N; longitude 119A degrees 35.4' E; elevation 890 m), China. A two-mirror coelostat with 100 cm focal length lens produced a 9.2 mm image of the Sun. The spectrograph using 140 cm focal length lens in Littrow mode and a grating with 600 lines per millimeter blazed at 2 mu m provided a dispersion of 30 m and 43 m per pixel in the fourth order around the green line and third order around the red line, respectively. Two Peltier cooled 1k x 1k CCD cameras, with a pixel size of 13 mu m square and 14-bit readout at 10 MHz operated in frame transfer mode, were used to obtain the time sequence spectra in two emission lines simultaneously. The duration of totality was 341 s, but we could get spectra for 270 s after a trial exposure at an interval of 5 s. We report here on the detection of intensity, velocity, and line width oscillations with periodicity in the range of 25 -50 s. These oscillations can be interpreted in terms of the presence of fast magnetoacoustic waves or torsional Alfv,n waves. The intensity ratios of green to red emission lines indicate the temperature of the corona to be 1.65 MK in the equatorial region and 1.40 MK in the polar region, relatively higher than the expected temperature during the low activity period. The width variation of the emission lines in different coronal structures suggests different physical conditions in different structures.

Uniform algebras generated by holomorphic and close-to-harmonic functions

The initial motivation for this paper is to discuss a more concrete approach to an approximation theorem of Axler and Shields, which says that the uniform algebra on the closed unit disc (D) over bar generated by z and h, where h is a nowhere-holomorphic harmonic function on D that is continuous up to partial derivative D, equals C((D) over bar). The abstract tools used by Axler and Shields make harmonicity of h an essential condition for their result. We use the concepts of plurisubharmonicity and polynomial convexity to show that, in fact, the same conclusion is reached if h is replaced by h + R, where R is a non-harmonic perturbation whose Laplacian is ``small'' in a certain sense.

Learning the global maximum with parameterized learning automata

A feedforward network composed of units of teams of parameterized learning automata is considered as a model of a reinforcement teaming system. The internal state vector of each learning automaton is updated using an algorithm consisting of a gradient following term and a random perturbation term. It is shown that the algorithm weakly converges to a solution of the Langevin equation implying that the algorithm globally maximizes an appropriate function. The algorithm is decentralized, and the units do not have any information exchange during updating. Simulation results on common payoff games and pattern recognition problems show that reasonable rates of convergence can be obtained.

A two timescale stochastic approximation scheme for simulation-based parametric optimization

A two timescale stochastic approximation scheme which uses coupled iterations is used for simulation-based parametric optimization as an alternative to traditional "infinitesimal perturbation analysis" schemes, It avoids the aggregation of data present in many other schemes. Its convergence is analyzed, and a queueing example is presented.

Carbon recombination lines between 34.5 and 770 MHz toward Cassiopeia A

We present observations of low-frequency recombination lines of carbon toward Cas A near 34.5 MHz (n similar to 575) using the Gauribidanur radio telescope and near 560 MHz (n similar to 225) and 770 MHz (n similar to 205) using the NRAO 140 foot (43 m) telescope in Greenbank. We also present high angular resolution (1') observations of the C270 alpha line near 332 MHz using the Very Large Array in B-configuration. A high signal-to-noise ratio spectrum is obtained at 34.5 MHz, which clearly shows a Voigt profile with distinct Lorentzian wings, resulting from significant pressure and radiation broadening at such high quantum numbers. The emission lines detected near 332, 550, and 770 MHz, on the other hand, are narrow and essentially Doppler-broadened. The measured Lorentzian width at 34.5 MHz constrains the allowed combinations of radiation temperature, electron density, and electron temperature in the line-forming region. Radiation broadening at 34.5 MHz places a lower limit of 115 pc on the separation between Cas A and the line-forming clouds. Modeling the variation in the integrated line-to-continuum ratio with frequency indicates that the region is likely to be associated with the cold atomic hydrogen component of the interstellar medium, and the physical properties of this region are likely to be T-e = 75 K, n(e) = 0.02 cm(-3), T-R100 = 3200 K, and n(H) T-e = 10,000 cm(-3) K. Comparison of the distribution of the C270 alpha recombination line emission across Cas A with that of (CO)-C-12 and H I also supports the above conclusion.

Optical studies of SN 2009jf: a Type Ib supernova with an extremely slow decline and aspherical signature

Optical UBVRI photometry and medium-resolution spectroscopy of the Type Ib supernova SN 2009jf, during the period from similar to -15 to +250 d, with respect to the B maximum are reported. The light curves are broad, with an extremely slow decline. The early post-maximum decline rate in the V band is similar to SN 2008D; however, the late-phase decline rate is slower than other Type Ib supernovae studied. With an absolute magnitude of M-V = -17.96 +/- 0.19 at peak, SN 2009jf is a normally bright supernova. The peak bolometric luminosity and the energy deposition rate via the 56Ni -> 56Co chain indicate that similar to 0.17+0.03(-0.03) M-circle dot of 56Ni was ejected during the explosion. The He i 5876 A line is clearly identified in the first spectrum of day similar to -15, at a velocity of similar to 16 000 km s-1. The O i] 6300-6364 A line seen in the nebular spectrum has a multipeaked and asymmetric emission profile, with the blue peak being stronger. The estimated flux in this line implies that greater than or similar to 1.34 M-circle dot oxygen was ejected. The slow evolution of the light curves of SN 2009jf indicates the presence of a massive ejecta. The high expansion velocity in the early phase and broader emission lines during the nebular phase suggest it to be an explosion with a large kinetic energy. A simple qualitative estimate leads to the ejecta mass of M-ej = 4-9 M-circle dot and kinetic energy E-K = 3-8 x 1051 erg. The ejected mass estimate is indicative of an initial main-sequence mass of greater than or similar to 20-25 M-circle dot.

Lopsidedness in WHISP galaxies II. Morphological lopsidedness

The distribution of stars and gas in many galaxies is asymmetric. This so-called lopsidedness is expected to significantly affect the dynamics and evolution of the disc, including the star formation activity. Here, we measure the degree of lopsidedness for the gas distribution in a selected sample of 70 galaxies from the Westerbork Hi Survey of Spiral and Irregular Galaxies. This complements our earlier work (Paper I) where the kinematic lopsidedness was derived for the same galaxies. The morphological lopsidedness is measured by performing a harmonic decomposition of the surface density maps. The amplitude of lopsidedness A(1), the fractional value of the first Fourier component, is typically quite high (about 0.1) within the optical disc and has a constant phase. Thus, lopsidedness is a common feature in galaxies and indicates a global mode. We measure A(1) out to typically one to four optical radii, sometimes even further. This is, on average, four times larger than the distance to which lopsidedness was measured in the past using near-IR as a tracer of the old stellar component, and therefore provides a new, more stringent constraint on the mechanism for the origin of lopsidedness. Interestingly, the value of A(1) saturates beyond the optical radius. Furthermore, the plot of A(1) versus radius shows fluctuations that we argue are due to local spiral features. We also try to explain the physical origin of this observed disc lopsidedness. No clear trend is found when the degree of lopsidedness is compared to a measure of the isolation or interaction probability of the sample galaxies. However, this does not rule out a tidal origin if the lopsidedness is long-lived. In addition, we find that the early-type galaxies tend to be more morphologically lopsided than the late-type galaxies. Both results together indicate that lopsidedness has a tidal origin.