An important tool for the development of the next generation of extremely large telescopes (ELTs) is a robust Systems Engineering (SE) methodology. GMACS is a first-generation multi-object spectrograph that will work at visible wavelengths on the Giant Magellan Telescope (GMT). In this paper, we discuss the application of SE to the design of next-generation instruments for ground-based astronomy and present the ongoing development of SE products for the GMACS spectrograph, currently in its Conceptual Design phase. SE provides the means to assist in the management of complex projects, and in the case of GMACS, to ensure its operational success, maximizing the scientific potential of GMT.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

We propose a new technique to study fast transitions during inflation, by studying the dynamics of quantum quenches in an $O(N)$ scalar field theory in de Sitter spacetime. We compute the time evolution of the system using a non-perturbative large-$N$ limit approach. We derive the self-consistent mass equation for several physically relevant transitions of the parameters of the theory, in a slow motion approximation. Our computations reveal that the effective mass after the quench evolves in the direction of recovering its value before the quench, but stopping at a different asymptotic value, in which the mass squared is strictly positive. Furthermore, we tentatively find situations in which the effective mass squared can be temporarily negative, thus breaking the $O(N)$ symmetry of the system for a certain time, only to then come back to a positive value, restoring the symmetry. We argue the relevance of our new method in a cosmological scenario.

Attempts to modify gravity in the infrared typically require a screening mechanism to ensure consistency with local tests of gravity. These screening mechanisms fit into three broad classes; we investigate theories which are capable of exhibiting more than one type of screening. Specifically, we focus on a simple model which exhibits both Vainshtein and kinetic screening. We point out that due to the two characteristic length scales in the problem, the type of screening that dominates depends on the mass of the sourcing object, allowing for different phenomenology at different scales. We consider embedding this double screening phenomenology in a broader cosmological scenario and show that the simplest examples that exhibit double screening are radiatively stable.

In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bigravity, a non-minimal matter coupling involving both metrics generically re-introduces the Boulware--Deser (BD) ghost. A non-minimal matter coupling via a simple, yet specific composite metric has been proposed, which eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formulation and for arbitrary spacetime dimensions, we show that this composite metric is the unique consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD ghost from matter quantum loop corrections.

If the amplitude of primordial gravitational waves is measured in the near-future, what could it tell us about bigravity? To address this question, we study massive bigravity theories by focusing on a region in parameter space which is safe from known instabilities. Similarly to investigations on late time constraints, we implicitly assume there is a successful implementation of the Vainshtein mechanism which guarantees that standard cosmological evolution is largely unaffected. We find that viable bigravity models are subject to far less stringent constraints than massive gravity, where there is only one set of (massive) tensor modes. In principle sensitive to the effective graviton mass at the time of recombination, we find that in our setup the primordial tensor spectrum is more responsive to the dynamics of the massless tensor sector rather than its massive counterpart. We further show there are intriguing windows in the parameter space of the theory which could potentially induce distinctive signatures in the B-modes spectrum.

We investigate the stability of a class of derivative theories known as $P(X)$ and Galileons against corrections generated by quantum effects. We use an exact renormalisation group approach to argue that these theories are stable under quantum corrections at all loops in regions where the kinetic term is large compared to the strong coupling scale. This is the regime of interest for screening or Vainshtein mechanisms, and in inflationary models that rely on large kinetic terms. Next, we clarify the role played by the symmetries. While symmetries protect the form of the quantum corrections, theories equipped with more symmetries do not necessarily have a broader range of scales for which they are valid. We show this by deriving explicitly the regime of validity of the classical solutions for $P(X)$ theories including Dirac-Born-Infeld (DBI) models, both in generic and for specific background field configurations. Indeed, we find that despite the existence of an additional symmetry, the DBI effective field theory has a regime of validity similar to an arbitrary $P(X)$ theory. We explore the implications of our results for both early and late universe contexts. Conversely, when applied to static and spherical screening mechanisms, we deduce that the regime of validity of typical power-law $P(X)$ theories is much larger than that of DBI.

We investigate the coupling to matter in ghost-free massive (bi-)gravity. When species in the matter sector couple covariantly to only one metric, we show that at one--loop these couplings do not spoil the special structure of the graviton potential. When the same species couples directly to both metrics we show that a ghost is present at the classical level and that loops destroy the special structure of the potential at an unacceptably low scale. We then propose a new `composite' effective metric built out of both metrics. When matter fields couple covariantly to this effective metric, the would be Boulware--Deser ghost is absent in different representative limits. At one--loop such couplings do not detune the special structure of the potential. We conjecture that matter can couple covariantly to that effective metric in all generality without introducing any Boulware--Deser ghost below a cut-off scale parametrically larger than the strong coupling scale. We also discuss alternative couplings to matter where the kinetic and potential terms of the matter field couple to different metrics. In both cases we discuss preliminary implications for cosmology.

Recently, several works have investigated the coupling to matter in ghost-free massive (bi-&multi-)gravity and a new effective coupling to matter has been proposed. In this note we clarify some confusion on the existence and the implications of a ghost above the strong coupling scale. We confirm that the standard constraint which is otherwise typically present in this type of theories disappears on generic backgrounds as soon as this new coupling is considered. This implies the re-emergence of the Boulware-Deser ghost. Nevertheless the absence of ghost in the decoupling limit implies that the cut-off scale (if identified with the scale at which the ghost enters) is higher than the strong coupling scale. Therefore there is a valid interesting region of applicability for these couplings at scales below the cut-off.

The latest results from PLANCK impose strong constraints on features in the spectrum of the curvature perturbations from inflation. We analyse the possibility of particle production induced by sharp turns of the trajectory in field space in inflation models with multiple fields. Although the evolution of the background fields can be altered by particle production, we find rather modest changes in the power spectrum even for the most extreme case in which the entire kinetic energy of the scalar fields is converted into particles.

The effective field theory (EFT) of inflation has become the preferred method for computing cosmological correlation functions of the curvature fluctuation, $\zeta$. It makes explicit use of the soft breaking of time diffeomorphisms by the inflationary background to organize the operators expansion in the action of the Goldstone mode $\pi$ associated with this breaking. Despite its ascendancy, there is another method for calculating $\zeta$ correlators, involving the direct calculation of the so-called Horndeski action order by order in powers of $\zeta$ and its derivatives. The question we address in this work is whether or not the $\zeta$ correlators calculated in these seemingly different ways are in fact the same. The answer is that the actions to cubic order in either set of variables do indeed give rise to the same $\zeta$ bispectra, but that to make this equivalence manifest requires a careful understanding of the non-linear transformations relating $\pi$ to $\zeta$ and how boundary terms in the actions are affected by imposing this relation. As a by product of our study we find that the calculations in the $\pi$ language can be simplified considerably in a way that allows us to use only the linear part of the $\pi-\zeta$ relation simply by changing the coefficients of some of the operators in the EFT. We also note that a proper accounting of the boundary terms will be of the greatest importance when computing the bispectrum for more general initial states than the Bunch-Davies one.

We derive the 'separate universe' method for the inflationary bispectrum, beginning directly from a field-theory calculation. We work to tree-level in quantum effects but to all orders in the slow-roll expansion, with masses accommodated perturbatively. Our method provides a systematic basis to account for novel sources of time-dependence in inflationary correlation functions, and has immediate applications. First, we use our result to obtain the correct matching prescription between the 'quantum' and 'classical' parts of the separate universe computation. Second, we elaborate on the application of this method in situations where its validity is not clear. As a by-product of our calculation we give the leading slow-roll corrections to the three-point function of field fluctuations on spatially flat hypersurfaces in a canonical, multiple-field model.

Until recently our knowledge of the primordial curvature perturbation was relatively modest. Ever since COBE delivered its map of data we know the scalar spectrum of primordial perturbations is approximately flat, with the power being only slightly stronger at larger scales. Most inflationary models predict an approximately scale-invariant spectrum, which therefore cannot be used as a distinctive signature. To distinguish between different inflationary microphysics we need to study higher point statistics of the primordial perturbation, which can encode non-gaussian data. In the first part of this thesis we study the bispectrum in all single-field models with a well-defined quantum field theory during a quasi-de Sitter inflationary phase. Any single-field models without ghost-like instabilities fall into this description: from canonical, to Dirac-Born-Infeld inflation and galileon inflation theories. We investigate the scale and shape- dependences of the bispectrum to next-order in the slow-roll approximation. We illustrate our results by applying them to different models and argue these corrections must be taken into account to keep the theoretical error below the observational precision set by the Planck satellite. We then explore the ability of using bispectrum shapes to distinguish between inflationary models more efficiently. We further extend the study of the bispectrum of single-field models beyond the slow-roll approximation, demanding the spectral index to be close to, but not exactly, unity. In the second part of this thesis we explore the process by which the universe is repopulated with matter particles at the end of a Dirac-Born-Infeld inflation phase. We place some mild bounds on the reheating temperature of these models. We argue that the constraints arising from the preheating analysis are complementary to those derived from the primordial perturbation.

We compute the one-loop quantum corrections to the potential of ghost-free massive gravity. We show how the mass of external matter fields contribute to the running of the cosmological constant, but do not change the ghost-free structure of the massive gravity potential at one-loop. When considering gravitons running in the loops, we show how the structure of the potential gets destabilized at the quantum level, but in a way which would never involve a ghost with a mass smaller than the Planck scale. This is done by explicitly computing the one-loop effective action and supplementing it with the Vainshtein mechanism. We conclude that to one-loop order the special mass structure of ghost-free massive gravity is technically natural.

We present the first 3+1-dimensional numerical simulations of scalar fields with non-minimal kinetic terms. As an example, we examine the existence and stability of preheating in the presence of a Dirac-Born-Infeld (DBI) inflaton coupled to a canonical matter field. The simulations represent the full non-linear theory in the presence of an expanding Universe. We show that parametric resonance in the matter field, along with self-resonance in the inflaton, repopulate the Universe with matter particles as efficiently as in traditional preheating.

We provide a formula for the scaling behaviour of the inflationary bispectrum in the 'squeezed' limit where one momentum becomes much smaller than the other two. This determines the scaling of the halo bias at low wavenumber and will be an important observable for the next generation of galaxy surveys. Our formula allows it to be predicted for the first time for a generic inflationary model with multiple light, canonically-normalized scalar fields.

A pressing problem in comparing inflationary models with observation is the accurate calculation of correlation functions. One approach is to evolve them using ordinary differential equations ("transport equations"), analogous to the Schwinger-Dyson hierarchy of in-out quantum field theory. We extend this approach to the complete set of momentum space correlation functions. A formal solution can be obtained using raytracing techniques adapted from geometrical optics. We reformulate inflationary perturbation theory in this language, and show that raytracing reproduces the familiar "delta N" Taylor expansion. Our method produces ordinary differential equations which allow the Taylor coefficients to be computed efficiently. We use raytracing methods to express the gauge transformation between field fluctuations and the curvature perturbation, zeta, in geometrical terms. Using these results we give a compact expression for the nonlinear gauge-transform part of fNL in terms of the principal curvatures of uniform energy-density hypersurfaces in field space.

If the expansion of the early Universe was not close to de Sitter, the statistical imprints of the primordial density perturbation on the cosmic microwave background can be quite different from those derived in slow-roll inflation. In this paper we study the inflationary signatures of all single-field models which are free of ghost-like instabilities. We allow for a rapid change of the Hubble parameter and the speed of sound of scalar fluctuations, in a way that is compatible with a nearly scale-invariant spectrum of perturbations, as supported by current cosmological observations. Our results rely on the scale-invariant approximation, which is different from the standard slow-roll approximation. We obtain the propagator of scalar fluctuations and compute the bispectrum, keeping next-order corrections proportional to the deviation of the spectral index from unity. These theories offer an explicit example where the shape and scale-dependences of the bispectrum are highly non-trivial whenever slow-roll is not a good approximation.

Artificial neural networks (ANN) have different applications in Astronomy, including data reduction and data mining. In this work we propose the use ANNs in the identification of stellar model solutions. We illustrate this method, by applying an ANN to the 0.8M$_\odot$ star CG Cyg B. Our ANN was trained using 60,000 different 0.8M$_\odot$ stellar models. With this approach we identify the models which reproduce CG Cyg B's position in the HR diagram. We observe a correlation between the model's initial metal and helium abundance which, in most cases, does not agree with a helium to metal enrichment ratio $\Delta$Y/$\Delta$Z=2. Moreover, we identify a correlation between the model's initial helium/metal abundance and both its age and mixing-length parameter. Additionally, every model found has a mixing-length parameter below 1.3. This means that CG Cyg B's mixing-length parameter is clearly smaller than the solar one. From this study we conclude that ANNs are well suited to deal with the degeneracy of model solutions of solar type stars.

We investigate the validity of effective field theory methods and the decoupling of heavy fields during inflation. Considering models of inflation in which the inflaton is coupled to a heavy (super-Hubble) degree of freedom initially in its vacuum state, we find that violations of decoupling are absent unless there is a breakdown of the slow-roll conditions. Next we allow for a temporary departure from inflation resulting in a period of non-adiabaticity during which effective field theory methods are known to fail. We find that the locality of the event and energy conservation lead to a tight bound on the size of the effects of the heavy field. We discuss the implications for the power spectrum and non-gaussianity, and comment on the connection with recent studies of the dynamics of multi-field inflation models. Our results further motivate the use of effective field theory methods to characterize cosmic inflation, and focus the question of observability of additional degrees of freedom during inflation to near the Hubble scale or below - as anticipated from the Wilsonian notions of decoupling and naturalness.

We re-examine the importance of slow-roll corrections during the evolution of cosmological perturbations in models of multi-field inflation. We find that in many instances the presence of light degrees of freedom leads to situations in which next to leading order slow-roll corrections become significant. Examples where we expect such corrections to be crucial include models in which modes exit the Hubble radius while the inflationary trajectory undergoes an abrupt turn in field space, or during a phase transition. We illustrate this with two examples -- hybrid inflation and double quadratic inflation. Utilizing both analytic estimates and full numerical results, we find that corrections can be as large as 20%. Our results have implications for many existing models in the literature, as these corrections must be included to obtain accurate observational predictions -- particularly given the level of accuracy expected from CMB experiments such as Planck

Galileon fields arise naturally from the decoupling limit of massive gravities, and possess special self-interactions which are protected by a spacetime generalization of Galilean symmetry. We briefly revisit the inflationary phenomenology of Galileon theories. Working from recent computations of the fluctuation Lagrangian to cubic order in the most general model with second-order equations of motion, we show that a distinct shape is present but with suppressed amplitude. A similar shape has been found in other higher-derivative models. It may be visible in a theory tuned to suppress the leading-order shapes, or if the overall bispectrum has large amplitude. Using a partial-wave expansion of the bispectrum, we suggest a possible origin for the frequent appearance of this shape. It follows that models with very disparate microphysics can produce very similar bispectra. We argue that it may be more profitable to distinguish these models by searching for relations between the amplitudes of these common shapes. We illustrate this method using the example of DBI and k-inflation.

Nongaussian statistics are a powerful discriminant between inflationary models, particularly those with noncanonical kinetic terms. Focusing on theories where the Lagrangian is an arbitrary Lorentz-invariant function of a scalar field and its first derivatives, we review and extend the calculation of the observable three-point function. We compute the "next-order" slow-roll corrections to the bispectrum in closed form, and obtain quantitative estimates of their magnitude in DBI and power-law k-inflation. In the DBI case our results enable us to estimate corrections from the shape of the potential and the warp factor: these can be of order several tens of percent. We track the possible sources of large logarithms which can spoil ordinary perturbation theory, and use them to obtain a general formula for the scale dependence of the bispectrum. Our result satisfies the next-order version of Maldacena's consistency condition and an equivalent consistency condition for the scale dependence. We identify a new bispectrum shape available at next-order, which is similar to a shape encountered in Galileon models. If fNL is sufficiently large this shape may be independently detectable.

We study how the universe reheats following an era of chaotic Dirac-Born-Infeld inflation, and compare the rate of particle production with that in models based on canonical kinetic terms. Particle production occurs through non-perturbative resonances whose structure is modified by the nonlinearities of the Dirac-Born-Infeld action. We investigate these modifications and show that the reheating process may be efficient. We estimate the initial temperature of the subsequent hot, radiation-dominated phase.

The ubiquitous role of the cyber-infrastructures, such as the WWW, provides myriad opportunities for machine learning and its broad spectrum of application domains taking advantage of digital communication. Pattern classification and feature extraction are among the first applications of machine learning that have received extensive attention. The most remarkable achievements have addressed data sets of moderate-to-large size. The 'data deluge' in the last decade or two has posed new challenges for AI researchers to design new, effective and accurate algorithms for similar tasks using ultra-massive data sets and complex (natural or synthetic) dynamical systems. We propose a novel principled approach to feature extraction in hybrid architectures comprised of humans and machines in networked communication, who collaborate to solve a pre-assigned pattern recognition (feature extraction) task. There are two practical considerations addressed below: (1) Human experts, such as plant biologists or astronomers, often use their visual perception and other implicit prior knowledge or expertise without any obvious constraints to search for the significant features, whereas machines are limited to a pre-programmed set of criteria to work with; (2) in a team collaboration of collective problem solving, the human experts have diverse abilities that are complementary, and they learn from each other to succeed in cognitively complex tasks in ways that are still impossible imitate by machines.

Understanding the structure of relationships between objects in a given database is one of the most important problems in the field of data mining. The structure can be defined for a set of single objects (clustering) or a set of groups of objects (network mapping). We propose a method for discovering relationships between individuals (single or groups) that is based on what we call the empirical topology, a system-theoretic measure of functional proximity. To illustrate the suitability and efficiency of the method, we apply it to an astronomical data base.

We study preheating in DBI hybrid inflation. Preheating occurs since the matter fields are non-minimally coupled to the inflaton. Despite the coupling being small, preheating happens as a narrow parametric resonance in which matter fields and inflaton modes are exponentially created. We show that, within the approximations used, the process is very efficient at late times. The resonance stage gives way to the decay of the inflaton into bosonic or fermionic degrees of freedom and the corresponding reheating temperature is determined.

Cosmological solutions for a pressureless fluid in the Brans-Dicke theory exhibit asymptotical accelerated phase for some range of values of the parameter $\omega$, interpolating a matter dominated phase and an inflationary phase. The effective gravitational coupling is negative. We test this model against the supernovae type Ia data. The fitting of the observational data is slightly better than that obtained from the $\Lambda CDM$ model. We make some speculations on how to reconcile the negative gravitational coupling in large scale with the local tests. Some considerations on the structure formation problem in this context are also made.

The possibility that the present acceleration of the universe is driven by a kind of viscous fluid is exploited. At background level this model is similar to the generalized Chaplygin gas model (GCGM). But, at perturbative level, the viscous fluid exhibits interesting properties. In particular the oscillations in the power spectrum that plagues the GCGM are not present. Possible fundamental descriptions for this viscous dark energy are discussed.

The generalized Chaplygin gas model is characterized by the equation of state $p = - \frac{A}{\rho^\alpha}$. It is generally stated that the case $\alpha = 0$ is equivalent to a model with cosmological constant and dust ($\Lambda CDM$). In this work we show that, if this is true for the background equations, this is not true for the perturbation equations. Hence, the mass spectrum predicted for both models may differ.