Models of environment and T_1 relaxation in Josephson Charge Qubits

A theoretical interpretation of the recent experiments of Astafiev et. al. on the T_1-relaxation rate in Josephson Charge Qubits is proposed. The experimentally observed reproducible nonmonotonic dependence of T_1 on the splitting E_J of the qubit levels suggests further specification of the previously proposed models of the background charge noise. From our point of view the most promising is the ``Andreev fluctuator'' model of the noise. In this model the fluctuator is a Cooper pair that tunnels from a superconductor and occupies a pair of localized electronic states. Within this model one can naturally explain both the average linear T_1(E_J) dependence and the irregular fluctuations. The role of fluctuators in the formation of strong resonant peaks in this dependence is also discussed.Comment: 4 pages, 3 figure

Non-anticommutative deformation of N=(1,1) hypermultiplets

We study the SO(4)xSU(2) invariant and N=(1,0) supersymmetry-preserving nilpotent (non-anticommutative) Moyal deformation of hypermultiplets interacting with an abelian gauge multiplet, starting from their off-shell formulation in Euclidean N=(1,1) harmonic superspace. The deformed version of a neutral or a charged hypermultiplet corresponds to the `adjoint' or the `fundamental' representation of the deformed U(1) gauge group on the superfields involved. The neutral hypermultiplet action is invariant under N=(2,0) supersymmetry and describes a deformed N=(2,2) gauge theory. For both the neutral and the charged hypermultiplet we present the corresponding component actions and explicitly give the Seiberg-Witten-type transformations to the undeformed component fields. Mass terms for the hypermultiplets can be generated via the Scherk-Schwarz mechanism and Fayet-Iliopoulos term in analogy to the undeformed case.Comment: 1+17 pages; v2: introduction and conclusions expanded, published versio

Nilpotent deformations of N=2 superspace

We investigate deformations of four-dimensional N=(1,1) euclidean superspace induced by nonanticommuting fermionic coordinates. We essentially use the harmonic superspace approach and consider nilpotent bi-differential Poisson operators only. One variant of such deformations (termed chiral nilpotent) directly generalizes the recently studied chiral deformation of N=(1/2,1/2) superspace. It preserves chirality and harmonic analyticity but generically breaks N=(1,1) to N=(1,0) supersymmetry. Yet, for degenerate choices of the constant deformation matrix N=(1,1/2) supersymmetry can be retained, i.e. a fraction of 3/4. An alternative version (termed analytic nilpotent) imposes minimal nonanticommutativity on the analytic coordinates of harmonic superspace. It does not affect the analytic subspace and respects all supersymmetries, at the expense of chirality however. For a chiral nilpotent deformation, we present non(anti)commutative euclidean analogs of N=2 Maxwell and hypermultiplet off-shell actions.Comment: 1+16 pages; v2: discussion of (pseudo)conjugations extended, version to appear in JHE

Algorithms for computing parsimonious evolutionary scenarios for genome evolution, the last universal common ancestor and dominance of horizontal gene transfer in the evolution of prokaryotes

BACKGROUND: Comparative analysis of sequenced genomes reveals numerous instances of apparent horizontal gene transfer (HGT), at least in prokaryotes, and indicates that lineage-specific gene loss might have been even more common in evolution. This complicates the notion of a species tree, which needs to be re-interpreted as a prevailing evolutionary trend, rather than the full depiction of evolution, and makes reconstruction of ancestral genomes a non-trivial task. RESULTS: We addressed the problem of constructing parsimonious scenarios for individual sets of orthologous genes given a species tree. The orthologous sets were taken from the database of Clusters of Orthologous Groups of proteins (COGs). We show that the phyletic patterns (patterns of presence-absence in completely sequenced genomes) of almost 90% of the COGs are inconsistent with the hypothetical species tree. Algorithms were developed to reconcile the phyletic patterns with the species tree by postulating gene loss, COG emergence and HGT (the latter two classes of events were collectively treated as gene gains). We prove that each of these algorithms produces a parsimonious evolutionary scenario, which can be represented as mapping of loss and gain events on the species tree. The distribution of the evolutionary events among the tree nodes substantially depends on the underlying assumptions of the reconciliation algorithm, e.g. whether or not independent gene gains (gain after loss after gain) are permitted. Biological considerations suggest that, on average, gene loss might be a more likely event than gene gain. Therefore different gain penalties were used and the resulting series of reconstructed gene sets for the last universal common ancestor (LUCA) of the extant life forms were analysed. The number of genes in the reconstructed LUCA gene sets grows as the gain penalty increases. However, qualitative examination of the LUCA versions reconstructed with different gain penalties indicates that, even with a gain penalty of 1 (equal weights assigned to a gain and a loss), the set of 572 genes assigned to LUCA might be nearly sufficient to sustain a functioning organism. Under this gain penalty value, the numbers of horizontal gene transfer and gene loss events are nearly identical. This result holds true for two alternative topologies of the species tree and even under random shuffling of the tree. Therefore, the results seem to be compatible with approximately equal likelihoods of HGT and gene loss in the evolution of prokaryotes. CONCLUSIONS: The notion that gene loss and HGT are major aspects of prokaryotic evolution was supported by quantitative analysis of the mapping of the phyletic patterns of COGs onto a hypothetical species tree. Algorithms were developed for constructing parsimonious evolutionary scenarios, which include gene loss and gain events, for orthologous gene sets, given a species tree. This analysis shows, contrary to expectations, that the number of predicted HGT events that occurred during the evolution of prokaryotes might be approximately the same as the number of gene losses. The approach to the reconstruction of evolutionary scenarios employed here is conservative with regard to the detection of HGT because only patterns of gene presence-absence in sequenced genomes are taken into account. In reality, horizontal transfer might have contributed to the evolution of many other genes also, which makes it a dominant force in prokaryotic evolution

Revealing the intensity of turbulent energy transfer in planetary atmospheres

Images of the giant planets Jupiter and Saturn show highly turbulent storms and swirling Q23 clouds that reflect the intensity of turbulence in their atmospheres. Quantifying planetary turbulence is inaccessible to conventional tools, however, since they require large quantities of spatially and temporally resolved data. Here we show, using experiments, observations, and simulations, that potential vorticity (PV) is a straightforward and universal diagnostic that can be used to estimate turbulent energy transfer in a stably stratified atmosphere. We use the conservation of PV to define a length scale, LM, representing a typical distance over which PV is mixed by planetary turbulence. LM increases as the turbulent intensity increases and can be estimated from any latitudinal PV profile. Using this principle, we estimate LM within Jupiter's and Saturn's tropospheres, showing for the first time that turbulent energy transfer in Saturn's atmosphere is four times less intense than Jupiter'

Cassini observations reveal a regime of zonostrophic macroturbulence on Jupiter

In December 2000, the Cassini fly-by near Jupiter delivered high-resolution images of Jupiter’s clouds over the entire planet in a band between 50°N and 50°S. Three daily-averaged two-dimensional velocity snapshots extracted from these images are used to perform spectral analysis of jovian atmospheric macroturbulence. A similar analysis is also performed on alternative data documented by Choi and Showman (Choi, D., Showman, A. [2011]. Icarus 216, 597–609), based on a different method of image processing. The inter-comparison of the products of both analyses ensures a better constraint of the spectral estimates. Both analyses reveal strong anisotropy of the kinetic energy spectrum. The zonal spectrum is very steep and most of the kinetic energy resides in slowly evolving, alternating zonal (west–east) jets, while the non-zonal, or residual spectrum obeys the Kolmogorov–Kraichnan law specific to two-dimensional turbulence in the range of the inverse energy cascade. The spectral data is used to estimate the inverse cascade rate ∊ and the zonostrophy index Rβ for the first time. Although both datasets yield somewhat different values of ∊, it is estimated to be in the range 0.5–1.0 × 10−5 m2 s−3. The ensuing values of Rβ ≳ 5 belong well in the range of zonostrophic turbulence whose threshold corresponds to Rβ ≃ 2.5. We infer that the large-scale circulation is maintained by an anisotropic inverse energy cascade. The removal of the Great Red Spot from both datasets has no significant effect upon either the spectra or the inverse cascade rate. The spectral data are used to compute the rate of the energy exchange, W, between the non-zonal structures and the large-scale zonal flow. It is found that instantaneous values of W may exceed ∊ by an order of magnitude. Previous numerical simulations with a barotropic model suggest that W and ∊ attain comparable values only after averaging of W over a sufficiently long time. Near-instantaneous values of W that have been routinely used to infer the rate of the kinetic energy supply to Jupiter’s zonal flow may therefore significantly overestimate ∊. This disparity between W and ∊ may resolve the long-standing conundrum of an unrealistically high rate of energy transfer to the zonal flow. The meridional diffusivity Kϕ in the regime of zonostrophic turbulence is given by an expression that depends on ∊. The value of Kϕ estimated from the spectra is compared against data from the dispersion of stratospheric gases and debris resulting from the Shoemaker-Levy 9 comet and Wesley asteroid impacts in 1994 and 2009 respectively. Not only is Kϕ found to be consistent with estimates for both impacts, but the eddy diffusivity found from observations appears to be scale-independent. This behaviour could be a consequence of the interaction between anisotropic turbulence and Rossby waves specific to the regime of zonostrophic macroturbulence

Large eddy simulation of two-dimensional isotropic turbulence

Large eddy simulation (LES) of forced, homogeneous, isotropic, two-dimensional (2D) turbulence in the energy transfer subrange is the subject of this paper. A difficulty specific to this LES and its subgrid scale (SGS) representation is in that the energy source resides in high wave number modes excluded in simulations. Therefore, the SGS scheme in this case should assume the function of the energy source. In addition, the controversial requirements to ensure direct enstrophy transfer and inverse energy transfer make the conventional scheme of positive and dissipative eddy viscosity inapplicable to 2D turbulence. It is shown that these requirements can be reconciled by utilizing a two-parametric viscosity introduced by Kraichnan (1976) that accounts for the energy and enstrophy exchange between the resolved and subgrid scale modes in a way consistent with the dynamics of 2D turbulence; it is negative on large scales, positive on small scales and complies with the basic conservation laws for energy and enstrophy. Different implementations of the two-parametric viscosity for LES of 2D turbulence were considered. It was found that if kept constant, this viscosity results in unstable numerical scheme. Therefore, another scheme was advanced in which the two-parametric viscosity depends on the flow field. In addition, to extend simulations beyond the limits imposed by the finiteness of computational domain, a large scale drag was introduced. The resulting LES exhibited remarkable and fast convergence to the solution obtained in the preceding direct numerical simulations (DNS) by Chekhlov et al. (1994) while the flow parameters were in good agreement with their DNS counterparts. Also, good agreement with the Kolmogorov theory was found. This LES could be continued virtually indefinitely. Then, a simplifiedComment: 34 pages plain tex + 18 postscript figures separately, uses auxilary djnlx.tex fil

Gas Giants

The gas giants (Jupiter and Saturn) and icy giants (Uranus and Neptune) are fluid planets with atmospheres primarily made of hydrogen and helium. The part of their atmospheres accessible to remote sensing occupies only a small fraction of their radii (0.05%). Clouds and hazes form around the 1 bar altitude pressure level and extend vertically, according to the thermochemical models, in a layer with a thickness of 200_500 km where temperature increases with depth (usually known as the "weather layer"). Clouds made of NH3, NH4SH, H2O (in Jupiter and Saturn), with the addition of CH4 (in Uranus and Neptune), cover the planet in stratified layers that are mixed with unknown hromophore agents. Dynamical phenomena in the weather layer shape different cloud patterns that define the visible appearance of these planets. In the thermal part of the spectrum clouds act as opacity sources providing brightness contrasts. The ensemble of cloud morphologies in terms of shapes, sizes and albedos allows their use as tracers of the atmospheric motions in the weather layer (Fig. 4.1). This is the main tool employed so far to study the winds on these fourplanets

Branch-and-lift algorithm for deterministic global optimization in nonlinear optimal control

This paper presents a branch-and-lift algorithm for solving optimal control problems with smooth nonlinear dynamics and potentially nonconvex objective and constraint functionals to guaranteed global optimality. This algorithm features a direct sequential method and builds upon a generic, spatial branch-and-bound algorithm. A new operation, called lifting, is introduced, which refines the control parameterization via a Gram-Schmidt orthogonalization process, while simultaneously eliminating control subregions that are either infeasible or that provably cannot contain any global optima. Conditions are given under which the image of the control parameterization error in the state space contracts exponentially as the parameterization order is increased, thereby making the lifting operation efficient. A computational technique based on ellipsoidal calculus is also developed that satisfies these conditions. The practical applicability of branch-and-lift is illustrated in a numerical example. © 2013 Springer Science+Business Media New York