We test if Refracted Gravity (RG) can describe the dynamics of disk galaxies without resorting to dark matter. RG is a classical theory of gravity where the standard Poisson equation is modified by the gravitational permittivity, $\epsilon$, a universal monotonic function of the local mass density. We use the rotation curves and the vertical velocity dispersions of 30 galaxies in the DiskMass Survey (DMS) to determine $\epsilon$. RG describes the kinematic profiles with mass-to-light ratios in agreement with SPS models, and disk thicknesses in agreement with observations, once observational biases are considered. Our results rely on setting the three free parameters of RG for each galaxy. However, we show that the differences of these parameters from galaxy to galaxy could be ascribed to statistical fluctuations. We adopt an approximate method to find a single set of parameters that might properly describe the kinematics of the entire sample and suggest the universality of $\epsilon$. We finally show that the RG models of the individual rotation curves can only partly describe the radial acceleration relation (RAR). Evidently, the RG models underestimate the observed accelerations by 0.1-0.3 dex at low Newtonian accelerations. Another problem is the strong correlations, at largely more than 5$\sigma$, between the residuals of the RAR models and three radially-dependent properties of galaxies, whereas the DMS data show considerably less significant correlations, at more than 4$\sigma$, for only two of them. These correlations might originate the non-null intrinsic scatter of the RG models, at odds with the observed intrinsic scatter of galaxy samples, different from DMS, which is consistent with 0. Further studies are required to assess if these discrepancies in the RAR originate from the DMS sample, which might not be ideal for deriving the RAR, or if they are genuine failures of RG.

A distinct population of planetary systems that contain dynamically isolated, Earth-size planets with orbital periods $P_\mathrm{orb}\sim 1\,$day was recently identified in an analysis of data from the Kepler planet candidate catalog. We argue that these objects could represent the remnant rocky cores of giant planets that arrived at the stellar vicinity on high-eccentricity orbits and were rapidly stripped of their gaseous envelopes after crossing their respective Roche limits (RLs) $a_\mathrm{R,p}$. In this picture, objects with $P_\mathrm{orb}\gtrsim1\,$day are mostly "early" cores that originated in planets with an initial periastron distance $a_\mathrm{per,0}\le a_\mathrm{R,p}$; they had high initial eccentricities but their orbits underwent fast tidal circularization after the cores were exposed. Objects with $P_\mathrm{orb}\lesssim 1\,$day are, by contrast, mostly "late" cores that originated in planets with $a_\mathrm{per,0}>a_\mathrm{R,p}$; these planets underwent orbital circularization to a radius $>a_\mathrm{per,0}$ but eventually reached $a_\mathrm{R,p}$ through tidal orbital decay. This picture naturally accounts for the spatial distribution of hot Earths and for the similarity of their inferred occurrence rate to that of hot Jupiters, and it fits well with the interpretation of the so-called sub-Jovian desert in the orbital-period--planetary-mass plane in terms of high-eccentricity planet migration to the vicinity of the RL.

In the high-eccentricity migration (HEM) scenario, close-in planets reach the vicinity of the central star on high-eccentricity orbits that become circularized---with a concomitant decrease in the semimajor axis---through a tidal interaction with the star. Giant planets that arrive with periastron distances that are smaller than the Roche limit $a_\mathrm{R}$ lose their gaseous envelopes, resulting in an inner edge to the surviving planets' distribution. The observational evidence for this effect, while extensive, is nonetheless somewhat ambiguous because of the effect of tidal orbital decay. Here we consider another key prediction of the HEM scenario---the existence of a spatial eccentricity gradient near the location where the circularization time becomes comparable to the planet's age for typical parameters. Previous studies already found evidence for this gradient and demonstrated that its properties are consistent with the circularization process being dominated by tidal dissipation in the planet (encapsulated by the tidal quality factor $Q^\prime_\mathrm{p}$). Our work extends these treatments by constructing explicit model distributions for comparison with the data and by carrying out backward-in-time integrations using observed system parameters. We show that circularization generally occurs outside the distribution's inner edge (which defines the boundary of the so-called sub-Jovian desert) and that typically $Q^\prime_\mathrm{p}\approx10^6$ in the circularization zone (to within a factor of 3). We also find tentative evidence for an eccentricity gradient in lower-mass planets, indicating that formation through HEM may be relevant down to Neptune scales.

Many of the observed spin--orbit alignment properties of exoplanets can be explained in the context of the primordial disk misalignment model, in which an initially aligned protoplanetary disk is torqued by a distant stellar companion on a misaligned orbit, resulting in a precessional motion that can lead to large-amplitude oscillations of the spin--orbit angle. We consider a variant of this model in which the companion is a giant planet with an orbital radius of a few au. Guided by the results of published numerical simulations, we model the dynamical evolution of this system by dividing the disk into inner and outer parts---separated at the location of the planet---that behave as distinct, rigid disks. We show that the planet misaligns the inner disk even as the orientation of the outer disk remains unchanged. In addition to the oscillations induced by the precessional motion, whose amplitude is larger the smaller the initial inner-disk-to-planet mass ratio, the spin--orbit angle also exhibits a secular growth in this case---driven by ongoing mass depletion from the disk---that becomes significant when the inner disk's angular momentum drops below that of the planet. Altogether, these two effects can produce significant misalignment angles for the inner disk, including retrograde configurations. We discuss these results within the framework of the Stranded Hot Jupiter scenario and consider their implications, including to the interpretation of the alignment properties of debris disks.

The recent advent of spatially resolved mm- and cm-wavelength polarimetry in protostellar accretion discs could help clarify the role of magnetic fields in the angular momentum transport in these systems. The best case to date is that of HL~Tau, where the inability to produce a good fit to the 1.25-mm data with a combination of vertical and azimuthal magnetic field components was interpreted as implying that centrifugally driven winds (CDWs) are probably not a significant transport mechanism on the $\sim 10^2\,$au scale probed by the observations. Using synthetic polarization maps of heuristic single-field-component discs and of a post-processed simulation of a wind-driving disc, we demonstrate that a much better fit to the data can be obtained if the radial field component, a hallmark of the CDW mechanism, dominates in the polarized emission region. A similar inference was previously made in modelling the far-infrared polarization map of the pc-scale dust ring in the Galactic centre. To reconcile this interpretation with theoretical models of protostellar discs, which indicate that the wind is launched from a comparatively high elevation above the mid-plane, we propose that most of the polarized emission originates -- with a high ($\ga 10$\%) intrinsic degree of polarization -- in small ($\la 0.1\,$mm) grains that remain suspended above the mid-plane, and that the bulk of the mm-wavelength emission is produced -- with low intrinsic polarization -- by larger grains that have settled to the mid-plane.

We investigate the proof of concept and the implications of \textit{refracted gravity}, a novel modified gravity aimed to solve the discrepancy between the luminous and the dynamical mass of cosmic structures without resorting to dark matter. Inspired by the behavior of electric fields in matter, refracted gravity introduces a gravitational permittivity that depends on the local mass density and modifies the standard Poisson equation. The resulting gravitational field can become more intense than the Newtonian field and can mimic the presence of dark matter. We show that the refracted gravitational field correctly describes (1) the rotation curves and the Tully-Fisher relation of disk galaxies; and (2) the observed temperature profile of the X-ray gas of galaxy clusters. According to these promising results, we conclude that refracted gravity deserves further investigation.

Transit and radial velocity observations indicate a dearth of sub-Jupiter--mass planets on short-period orbits, outlined roughly by two oppositely sloped lines in the period--mass plane. We interpret this feature in terms of high-eccentricity migration of planets that arrive in the vicinity of the Roche limit, where their orbits are tidally circularized, long after the dispersal of their natal disk. We demonstrate that the two distinct segments of the boundary are a direct consequence of the different slopes of the empirical mass--radius relation for small and large planets, and show that this relation also fixes the mass coordinate of the intersection point. The period coordinate of this point, as well as the detailed shape of the lower boundary, can be reproduced with a plausible choice of a key parameter in the underlying migration model. The detailed shape of the upper boundary, on the other hand, is determined by the post-circularization tidal exchange of angular momentum with the star and can be reproduced with a stellar tidal quality factor $Q^\prime_*\sim10^6$.

In this paper we focus on the construction of stellar outflow models emerging from a polar coronal hole-type region surrounded by a magnetosphere in the equatorial regions during phases of quiescent accretion. The models are based on initial analytical solutions. We adopt a meridionally self-similar solution of the time-independent and axisymmetric MHD equations which describes effectively a jet originating from the corona of a star. We modify appropriately this solution in order to incorporate a physically consistent stellar magnetosphere. We find that the closed fieldline region may exhibit different behaviour depending on the associated boundary conditions and the distribution of the heat flux. However, the stellar jet in all final equilibrium states is very similar to the analytical one prescribed in the initial conditions. When the initial net heat flux is maintained, the magnetosphere takes the form of a dynamical helmet streamer with a quasi steady state slow magnetospheric wind. With no heat flux, a small dead zone forms around the equator, surrounded by a weak streamer, source of a low mass flux magnetospheric ejection. An equilibrium with no streamer and a helmet shaped dead zone is found when the coronal hole is larger and extends up to a reduced closed fieldline static zone. The study demonstrates the robustness of self-similar solutions and their ability to provide quasi-steady jet simulations with realistic boundaries and surrounded by equatorial magnetospheres.

Models of planet formation and evolution predict that giant planets form efficiently in protoplanetary disks, that most of these migrate rapidly to the disk's inner edge, and that, if the arriving planet's mass is $\lesssim$ Jupiter's mass, it could remain stranded near that radius. We argue that such planets would be ingested by tidal interaction with the host star on a timescale $\lesssim1\,$Gyr, and that, in the case of a solar-type host, this would cause the stellar spin to approach the direction of the ingested planet's orbital axis even if the two were initially highly misaligned. Primordially misaligned stars whose effective temperatures are $\gtrsim6250\,$K cannot be realigned in this way because, in contrast with solar-type hosts, their angular momenta are typically higher than the orbital angular momentum of the ingested planet as a result of inefficient magnetic braking and of a comparatively large moment of inertia. Hot Jupiters located farther out from the star can contribute to this process, but their effect is weaker because the tidal interaction strength decreases rapidly with increasing semimajor axis. We demonstrate that, if $\sim50\%$ of planetary systems harbored a stranded hot Jupiter, this scenario can in principle account for (1) the good alignment exhibited by planets around cool stars irrespective of the planet's mass or orbital period, (2) the prevalence of misaligned planets around hot stars, (3) the apparent upper bound on the mass of hot Jupiters on retrograde orbits, and (4) the inverse correlation between stellar spin periods and hot-Jupiter masses.

Hot Jupiters are subject to strong irradiation from the host stars and, as a consequence, they do evaporate. They can also interact with the parent stars by means of tides and magnetic fields. Both phenomena have strong implications for the evolution of these systems. Here we present time resolved spectroscopy of HD~189733 observed with the Cosmic Origin Spectrograph (COS) on board to HST. The star has been observed during five consecutive HST orbits, starting at a secondary transit of the planet ($\phi$ ~0.50-0.63). Two main episodes of variability of ion lines of Si, C, N and O are detected, with an increase of line fluxes. Si IV lines show the highest degree of variability. The FUV variability is a signature of enhanced activity in phase with the planet motion, occurring after the planet egress, as already observed three times in X-rays. With the support of MHD simulations, we propose the following interpretation: a stream of gas evaporating from the planet is actively and almost steadily accreting onto the stellar surface, impacting at $70-90\deg$ ahead of the sub-planetary point.

Close-in exoplanets interact with their host stars gravitationally as well as via their magnetized plasma outflows. The rich dynamics that arises may result in distinct observable features. Our objective is to study and classify the morphology of the different types of interaction that can take place between a giant close-in planet (a Hot Jupiter) and its host star, based on the physical parameters that characterize the system. We perform 3D magnetohydrodynamic numerical simulations to model the star--planet interaction, incorporating a star, a Hot Jupiter, and realistic stellar and planetary outflows. We explore a wide range of parameters and analyze the flow structures and magnetic topologies that develop. Our study suggests the classification of star--planet interactions into four general types, based on the relative magnitudes of three characteristic length scales that quantify the effects of the planetary magnetic field, the planetary outflow, and the stellar gravitational field in the interaction region. We describe the dynamics of these interactions and the flow structures that they give rise to, which include bow shocks, cometary-type tails, and inspiraling accretion streams. We point out the distinguishing features of each of the classified cases and discuss some of their observationally relevant properties. The magnetized interactions of star--planet systems can be categorized, and their general morphologies predicted, based on a set of basic stellar, planetary, and orbital parameters.

We investigate the properties of X-ray emission from accretion shocks in classical T Tauri stars (CTTSs), generated where the infalling material impacts the stellar surface. Both observations and models of the accretion process reveal several aspects that are still unclear: the observed X-ray luminosity in accretion shocks is below the predicted value, and the density versus temperature structure of the shocked plasma, with increasing densities at higher temperature, deduced from the observations, is at odds with that proposed in the current picture of accretion shocks. To address these open issues we investigate whether a correct treatment of the local absorption by the surrounding medium is crucial to explain the observations. To this end, we describe the impact of an accretion stream on a CTTS by considering a magnetohydrodynamic model. From the model results we synthesize the X-ray emission from the accretion shock by producing maps and spectra. We perform density and temperature diagnostics on the synthetic spectra, and we directly compare the results with the observations. Our model shows that the X-ray fluxes inferred from the emerging spectra are lower than expected because of the complex local absorption by the optically thick material of the chromosphere and of the unperturbed stream. Moreover, our model including the effects of local absorption explains in a natural way the apparently puzzling pattern of density versus temperature observed in the X-ray emission from accretion shocks.

Young stellar object observations suggest that some jets rotate in the opposite direction with respect to their disk. In a recent study, Sauty et al. (2012) have shown that this does not contradict the magnetocentrifugal mechanism that is believed to launch such outflows. Signatures of motions transverse to the jet axis and in opposite directions have recently been measured in M87 (Meyer et al. 2013). One possible interpretation of this motion is the one of counter rotating knots. Here, we extend our previous analytical derivation of counter-rotation to relativistic jets, demonstrating that counter-rotation can indeed take place under rather general conditions. We show that both the magnetic field and a non-negligible enthalpy are necessary at the origin of counter-rotating outflows, and that the effect is associated with a transfer of energy flux from the matter to the electromagnetic field. This can be realized in three cases : if a decreasing enthalpy causes an increase of the Poynting flux, if the flow decelerates, or, if strong gradients of the magnetic field are present. An illustration of the involved mechanism is given by an example of relativistic MHD jet simulation.

Astronomical observations, analytical solutions and numerical simulations have provided the building blocks to formulate the current theory of young stellar object jets. Although each approach has made great progress independently, it is only during the last decade that significant efforts are being made to bring the separate pieces together. Building on previous work that combined analytical solutions and numerical simulations, we apply a sophisticated cooling function to incorporate optically thin energy losses in the dynamics. On the one hand, this allows a self-consistent treatment of the jet evolution and on the other, it provides the necessary data to generate synthetic emission maps. Firstly, analytical disk and stellar outflow solutions are properly combined to initialize numerical two-component jet models inside the computational box. Secondly, magneto-hydrodynamical simulations are performed in 2.5D, following properly the ionization and recombination of a maximum of $29$ ions. Finally, the outputs are post-processed to produce artificial observational data. The first two-component jet simulations, based on analytical models, that include ionization and optically thin radiation losses demonstrate promising results for modeling specific young stellar object outflows. The generation of synthetic emission maps provides the link to observations, as well as the necessary feedback for the further improvement of the available models.

(abridged) AIMS. We investigate the dynamics and stability of post-shock plasma streaming along nonuniform stellar magnetic fields at the impact region of accretion columns. We study how the magnetic field configuration and strength determine the structure, geometry, and location of the shock-heated plasma. METHODS. We model the impact of an accretion stream onto the chromosphere of a CTTS by 2D axisymmetric magnetohydrodynamic simulations. Our model takes into account the gravity, the radiative cooling, and the magnetic-field-oriented thermal conduction. RESULTS. The structure, stability, and location of the shocked plasma strongly depend on the configuration and strength of the magnetic field. For weak magnetic fields, a large component of B may develop perpendicular to the stream at the base of the accretion column, limiting the sinking of the shocked plasma into the chromosphere. An envelope of dense and cold chromospheric material may also develop around the shocked column. For strong magnetic fields, the field configuration determines the position of the shock and its stand-off height. If the field is strongly tapered close to the chromosphere, an oblique shock may form well above the stellar surface. In general, a nonuniform magnetic field makes the distribution of emission measure vs. temperature of the shocked plasma lower than in the case of uniform magnetic field. CONCLUSIONS. The initial strength and configuration of the magnetic field in the impact region of the stream are expected to influence the chromospheric absorption and, therefore, the observability of the shock-heated plasma in the X-ray band. The field strength and configuration influence also the energy balance of the shocked plasma, its emission measure at T > 1 MK being lower than expected for a uniform field. The above effects contribute in underestimating the mass accretion rates derived in the X-ray band.

Context. Theoretical arguments and numerical simulations of radiative shocks produced by the impact of the accreting gas onto young stars predict quasi-periodic oscillations in the emitted radiation. However, observational data do not show evidence of such periodicity. Aims. We investigate whether physically plausible perturbations in the accretion column or in the chromosphere could disrupt the shock structure influencing the observability of the oscillatory behavior. Methods. We performed local 2D magneto-hydrodynamical simulations of an accretion shock impacting a chromosphere, taking optically thin radiation losses and thermal conduction into account. We investigated the effects of several perturbation types, such as clumps in the accretion stream or chromospheric fluctuations, and also explored a wide range of plasma-\beta values. Results. In the case of a weak magnetic field, the post-shock region shows chaotic motion and mixing, smoothing out the perturbations and retaining a global periodic signature. On the other hand, a strong magnetic field confines the plasma in flux tubes, which leads to the formation of fibrils that oscillate independently. Realistic values for the amplitude, length, and time scales of the perturbation are capable of bringing the fibril oscillations out of phase, suppressing the periodicity of the emission. Conclusions. The strength of a locally uniform magnetic field in YSO accretion shocks determines the structure of the post-shock region, namely, whether it will be somewhat homogeneous or if it will split up to form a collection of fibrils. In the second case, the size and shape of the fibrils is found to depend strongly on the plasma-\beta value but not on the perturbation type. Therefore, the actual value of the protostellar magnetic field is expected to play a critical role in the time dependence of the observable emission.

Rotation measurement in jets from T Tauri stars is a rather difficult task. Some jets seem to be rotating in a direction opposite to that of the underlying disk, although it is not yet clear if this affects the totality or part of the outflows. On the other hand, Ulysses data also suggest that the solar wind may rotate in two opposite ways between the northern and southern hemispheres. We show that this result is not as surprising as it may seem and that it emerges naturally from the ideal MHD equations. Specifically, counterrotating jets neither contradict the magnetocentrifugal driving of the flow nor prevent extraction of angular momentum from the disk. The demonstration of this result is shown by combining the ideal MHD equations for steady axisymmetric flows. Provided that the jet is decelerated below some given threshold beyond the Alfven surface, the flow will change its direction of rotation locally or globally. Counterrotation is also possible for only some layers of the outflow at specific altitudes along the jet axis. We conclude that the counterrotation of winds or jets with respect to the source, star or disk, is not in contradiction with the magnetocentrifugal driving paradigm. This phenomenon may affect part of the outflow, either in one hemisphere, or only in some of the outflow layers. From a time-dependent simulation, we illustrate this effect and show that it may not be permanent.

It is a well established fact that some YSO jets (e.g. RW Aur) display different propagation speeds between their blue and red shifted parts, a feature possibly associated with the central engine or the environment in which the jet propagates. In order to understand the origin of asymmetric YSO jet velocities, we investigate the efficiency of two candidate mechanisms, one based on the intrinsic properties of the system and one based on the role of the external medium. In particular, a parallel or anti-parallel configuration between the protostellar magnetosphere and the disk magnetic field is considered and the resulting dynamics are examined both in an ideal and a resistive magneto-hydrodynamical (MHD) regime. Moreover, we explore the effects of a potential difference in the pressure of the environment, as a consequence of the non-uniform density distribution of molecular clouds. Ideal and resistive axisymmetric numerical simulations are carried out for a variety of models, all of which are based on a combination of two analytical solutions, a disk wind and a stellar outflow. We find that jet velocity asymmetries can indeed occur both when multipolar magnetic moments are present in the star-disk system as well as when non-uniform environments are considered. The latter case is an external mechanism that can easily explain the large time scale of the phenomenon, whereas the former one naturally relates it to the YSO intrinsic properties. [abridged]

Recent observations as well as theoretical studies of YSO jets suggest the presence of two steady components: a disk wind type outflow needed to explain the observed high mass loss rates and a stellar wind type outflow probably accounting for the observed stellar spin down. In this framework, we construct numerical two-component jet models by properly mixing an analytical disk wind solution with a complementary analytically derived stellar outflow. Their combination is controlled by both spatial and temporal parameters, in order to address different physical conditions and time variable features. We study the temporal evolution and the interaction of the two jet components on both small and large scales. The simulations reach steady state configurations close to the initial solutions. Although time variability is not found to considerably affect the dynamics, flow fluctuations generate condensations, whose large scale structures have a strong resemblance to observed YSO jet knots.

Theoretical arguments along with observational data of YSO jets suggest the presence of two steady components: a disk wind type outflow needed to explain the observed high mass loss rates and a stellar wind type outflow probably accounting for the observed stellar spin down. Each component's contribution depends on the intrinsic physical properties of the YSO-disk system and its evolutionary stage. The main goal of this paper is to understand some of the basic features of the evolution, interaction and co-existence of the two jet components over a parameter space and when time variability is enforced. Having studied separately the numerical evolution of each type of the complementary disk and stellar analytical wind solutions in Paper I of this series, we proceed here to mix together the two models inside the computational box. The evolution in time is performed with the PLUTO code, investigating the dynamics of the two-component jets, the modifications each solution undergoes and the potential steady state reached.

(Abridged) Finite radius accretion disks are a strong candidate for launching astrophysical jets from their inner parts and disk-winds are considered as the basic component of such magnetically collimated outflows. The only available analytical MHD solutions for describing disk-driven jets are those characterized by the symmetry of radial self-similarity. Radially self-similar MHD models, in general, have two geometrical shortcomings, a singularity at the jet axis and the non-existence of an intrinsic radial scale, i.e. the jets formally extend to radial infinity. Hence, numerical simulations are necessary to extend the analytical solutions towards the axis and impose a physical boundary at finite radial distance. We focus here on studying the effects of imposing an outer radius of the underlying accreting disk (and thus also of the outflow) on the topology, structure and variability of a radially self-similar analytical MHD solution. The initial condition consists of a hybrid of an unchanged and a scaled-down analytical solution, one for the jet and the other for its environment. In all studied cases, we find at the end steady two-component solutions.

Observations of collimated outflows in young stellar objects indicate that several features of the jets can be understood by adopting the picture of a two-component outflow, wherein a central stellar component around the jet axis is surrounded by an extended disk-wind. The precise contribution of each component may depend on the intrinsic physical properties of the YSO-disk system as well as its evolutionary stage. In this context, the present article starts a systematic investigation of two-component jet models via time-dependent simulations of two prototypical and complementary analytical solutions, each closely related to the properties of stellar-outflows and disk-winds. These models describe a meridionally and a radially self-similar exact solution of the steady-state, ideal hydromagnetic equations, respectively. By using the PLUTO code to carry out the simulations, the study focuses on the topological stability of each of the two analytical solutions, which are successfully extended to all space by removing their singularities. In addition, their behavior and robustness over several physical and numerical modifications is extensively examined. It is found that radially self-similar solutions (disk-winds) always reach a final steady-state while maintaining all their well-defined properties. The different ways to replace the singular part of the solution around the symmetry axis, being a first approximation towards a two-component outflow, lead to the appearance of a shock at the super-fast domain corresponding to the fast magnetosonic separatrix surface. Conversely, the asymptotic configuration and the stability of meridionally self-similar models (stellar-winds) is related to the heating processes at the base of the wind.

We present a new numerical code, PLUTO, for the solution of hypersonic flows in 1, 2 and 3 spatial dimensions and different systems of coordinates. The code provides a multi-physics, multi-algorithm modular environment particularly oriented towards the treatment of astrophysical flows in presence of discontinuities. Different hydrodynamic modules and algorithms may be independently selected to properly describe Newtonian, relativistic, MHD or relativistic MHD fluids. The modular structure exploits a general framework for integrating a system of conservation laws, built on modern Godunov-type shock-capturing schemes. Although a plethora of numerical methods has been successfully developed over the past two decades, the vast majority shares a common discretization recipe, involving three general steps: a piecewise polynomial reconstruction followed by the solution of Riemann problems at zone interfaces and a final evolution stage. We have checked and validated the code against several benchmarks available in literature. Test problems in 1, 2 and 3 dimensions are discussed.