Identification of nonlinear systems using generalized kernel models

Nonlinear system identification is considered using a generalized kernel regression model. Unlike the standard kernel model, which employs a fixed common variance for all the kernel regressors, each kernel regressor in the generalized kernel model has an individually tuned diagonal covariance matrix that is determined by maximizing the correlation between the training data and the regressor using a repeated guided random search based on boosting optimization. An efficient construction algorithm based on orthogonal forward regression with leave-one-out (LOO) test statistic and local regularization (LR) is then used to select a parsimonious generalized kernel regression model from the resulting full regression matrix. The proposed modeling algorithm is fully automatic and the user is not required to specify any criterion to terminate the construction procedure. Experimental results involving two real data sets demonstrate the effectiveness of the proposed nonlinear system identification approach.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Simulated and observed fluxes of sensible and latent heat and CO2 at the WLEF-TV tower using SiB2.5

Three years of meteorological data collected at the WLEF-TV tower were used to drive a revised version of the Simple Biosphere (SiB 2.5) Model. Physiological properties and vegetation phenology were specified from satellite imagery. Simulated fluxes of heat, moisture, and carbon were compared to eddy covariance measurements taken onsite as a means of evaluating model performance on diurnal, synoptic, seasonal, and interannual time scales. The model was very successful in simulating variations of latent heat flux when compared to observations, slightly less so in the simulation of sensible heat flux. The model overestimated peak values of sensible heat flux on both monthly and diurnal scales. There was evidence that the differences between observed and simulated fluxes might be linked to wetlands near the WLEF tower, which were not present in the SiB simulation. The model overestimated the magnitude of the net ecosystem exchange of CO2 in both summer and winter. Mid-day maximum assimilation was well represented by the model, but late afternoon simulations showed excessive carbon uptake due to misrepresentation of within-canopy shading in the model. Interannual variability was not well simulated because only a single year of satellite imagery was used to parameterize the model.

Impact of city changes and weather on anthropogenic heat flux in Europe 1995–2015

How people live, work, move from place to place, consume and the technologies they use all affect heat emissions in a city which influences urban weather and climate. Here we document changes to a global anthropogenic heat flux (QF) model to enhance its spatial (30′′ × 30′′ to 0.5° × 0.5°) resolution and temporal coverage (historical, current and future). QF is estimated across Europe (1995–2015), considering changes in temperature, population and energy use. While on average QF is small (of the order 1.9–4.6 W m−2 across all the urban areas of Europe), significant spatial variability is documented (maximum 185 W m−2). Changes in energy consumption due to changes in climate are predicted to cause a 13% (11%) increase in QF on summer (winter) weekdays. The largest impact results from changes in temperature conditions which influences building energy use; for winter, with the coldest February on record, the mean flux for urban areas of Europe is 4.56 W m−2 and for summer (warmest July on record) is 2.23 W m−2. Detailed results from London highlight the spatial resolution used to model the QF is critical and must be appropriate for the application at hand, whether scientific understanding or decision making.

Curl free pressure gradients over orography in a solution of the fully compressible Euler equations with implicit treatment of acoustic and gravity waves

Steep orography can cause noisy solutions and instability in models of the atmosphere. A new technique for modelling flow over orography is introduced which guarantees curl free gradients on arbitrary grids, implying that the pressure gradient term is not a spurious source of vorticity. This mimetic property leads to better hydrostatic balance and better energy conservation on test cases using terrain following grids. Curl-free gradients are achieved by using the co-variant components of velocity over orography rather than the usual horizontal and vertical components. In addition, gravity and acoustic waves are treated implicitly without the need for mean and perturbation variables or a hydrostatic reference profile. This enables a straightforward description of the implicit treatment of gravity waves. Results are presented of a resting atmosphere over orography and the curl-free pressure gradient formulation is advantageous. Results of gravity waves over orography are insensitive to the placement of terrain-following layers. The model with implicit gravity waves is stable in strongly stratified conditions, with N∆t up to at least 10 (where N is the Brunt-V ̈ais ̈al ̈a frequency). A warm bubble rising over orography is simulated and the curl free pressure gradient formulation gives much more accurate results for this test case than a model without this mimetic property.

Could the cosmic acceleration be transient? A cosmographic evaluation

A possible slowing down of the cosmic expansion is investigated through a cosmographic approach. By expanding the luminosity distance to fourth order and fitting the SN Ia data from the most recent compilations (Union, Constitution and Union 2), the marginal likelihood distributions for the deceleration parameter today suggest a recent reduction of the cosmic acceleration and indicate that there is a considerable probability for q(0) > 0. Also in contrast to the prediction of the Lambda CDM model, the cosmographic q(z) reconstruction permits a cosmic expansion history where the cosmic acceleration could already have peaked and be presently slowing down, which would imply that the recent accelerated expansion of the universe is a transient phenomenon. It is also shown that to describe a transient acceleration the luminosity distance needs to be expanded at least to fourth order. The present cosmographic results depend neither on the validity of general relativity nor on the matter-energy contents of the universe.

Instabilities in thermal gravity with a cosmological constant

It is shown that in quantum gravity at finite temperature, the effective potential evaluated in the tadpole approximation can have a local minimum below a certain critical temperature. However, when the leading higher order thermal loop corrections are included, one finds that no static solution exists at high temperature. (C) 2008 Elsevier B.V. All rights reserved.

An alternative approach for general covariant Horava-Lifshitz gravity and matter coupling

Recently, in [3] Horava and Melby-Thompson proposed a nonrelativistic gravity theory with extended gauge symmetry that is free of the spin-0 graviton. We propose a minimal substitution recipe to implement this extended gauge symmetry which reproduces the results obtained by them. Our prescription has the advantage of being manifestly gauge invariant and immediately generalizable to other fields, like matter. We briefly discuss the coupling of gravity with scalar and vector fields found by our method. We show also that the extended gauge invariance in gravity does not force the value of. to be lambda = 1 as claimed in [3]. However, the spin-0 graviton is eliminated even for general lambda.

Action of the gravitational field on the dynamical Casimir effect

In this paper, we analyze the action of the gravitational field on the dynamical Casimir effect. We consider a massless scalar field confined in a cuboid cavity placed in a gravitational field described by a static and diagonal metric. With one of the plane mirrors of the cavity allowed to move, we compute the average number of particles created inside the cavity by means of the Bogoliubov coefficients computed through perturbative expansions. We apply our result to the case of an oscillatory motion of the mirror, assuming a weak gravitational field described by the Schwarzschild metric. The regime of parametric amplification is analyzed in detail, demonstrating that our computed result for the mean number of particles created agrees with specific associated cases in the literature. Our results, obtained in the framework of the perturbation theory, are restricted, under resonant conditions, to a short-time limit.

Hamiltonian constraint analysis of vector field theories with spontaneous Lorentz symmetry breaking

Recent investigations of various quantum-gravity theories have revealed a variety of possible mechanisms that lead to Lorentz violation. One of the more elegant of these mechanisms is known as Spontaneous Lorentz Symmetry Breaking (SLSB), where a vector or tensor field acquires a nonzero vacuum expectation value. As a consequence of this symmetry breaking, massless Nambu-Goldstone modes appear with properties similar to the photon in Electromagnetism. This thesis considers the most general class of vector field theories that exhibit spontaneous Lorentz violation-known as bumblebee models-and examines their candidacy as potential alternative explanations of E&M, offering the possibility that Einstein-Maxwell theory could emerge as a result of SLSB rather than of local U(1) gauge invariance. With this aim we employ Dirac's Hamiltonian Constraint Analysis procedure to examine the constraint structures and degrees of freedom inherent in three candidate bumblebee models, each with a different potential function, and compare these results to those of Electromagnetism. We find that none of these models share similar constraint structures to that of E&M, and that the number of degrees of freedom for each model exceeds that of Electromagnetism by at least two, pointing to the potential existence of massive modes or propagating ghost modes in the bumblebee theories.

Gravitational constraints of dS branes in an AdS Einstein-Brans-Dicke bulk

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Consistent modified gravity: dark energy, acceleration and the absence of cosmic doomsday

We discuss modified gravity which includes negative and positive powers of curvature and provides gravitational dark energy. It is shown that in GR plus a term containing a negative power of curvature, cosmic speed-up may be achieved while the effective phantom phase (with w less than -1) follows when such a term contains a fractional positive power of curvature. Minimal coupling with matter makes the situation more interesting: even 1/R theory coupled with the usual ideal fluid may describe the (effective phantom) dark energy. The account of the R(2) term (consistent modified gravity) may help to escape cosmic doomsday.

Quantization of the superstring in Ramond-Ramond backgrounds

Sigma model actions are constructed for the type II superstring compactified to four-and six-dimensional curved backgrounds which can contain non-vanishing Ramond-Ramond fields. These actions are N = 2 worldsheet superconformally invariant and can be covariantly quantized preserving manifest spacetime supersymmetry. They are constructed using a hybrid Version of superstring variables which combines features of the Ramond-Neveu-Schwarz and Green-Schwarz formalisms. For the AdS(2) x S-2 and AdS(3) x S-3 backgrounds, these actions differ from the classical Green-Schwarz actions by a crucial kinetic term for the fermions. Parts of this work have been done in collaborations with M Bershadsky, T Hauer, W Siegel, C Vafa, E Witten, S Zhukov and B Zwiebach.

Renormalized new solutions for the massless Thirring model

We present a nonperturbative study of the (1 + 1)-dimensional massless Thirring model by using path integral methods. The regularization ambiguities - coming from the computation of the fermionic determinant - allow to find new solution types for the model. At quantum level the Ward identity for the 1PI 2-point function for the fermionic current separates such solutions in two phases or sectors, the first one has a local gauge symmetry that is implemented at quantum level and the other one without this symmetry. The symmetric phase is a new solution which is unrelated to the previous studies of the model and, in the nonsymmetric phase there are solutions that for some values of the ambiguity parameter are related to well-known solutions of the model. We construct the Schwinger-Dyson equations and the Ward identities. We make a detailed analysis of their UV divergence structure and, after, we perform a nonperturbative regularization and renormalization of the model.

Regge calculus in teleparallel gravity

In the context of the teleparallel equivalent of general relativity, the Weitzenbock manifold is considered as the limit of a suitable sequence of discrete lattices composed of an increasing number of smaller and smaller simplices, where the interior of each simplex (Delaunay lattice) is assumed to be flat. The link lengths l between any pair of vertices serve as independent variables, so that torsion turns out to be localized in the two-dimensional hypersurfaces (dislocation triangle, or hinge) of the lattice. Assuming that a vector undergoes a dislocation in relation to its initial position as it is parallel transported along the perimeter of the dual lattice (Voronoi polygon), we obtain the discrete analogue of the teleparallel action, as well as the corresponding simplicial vacuum field equations.