Measure the Sloan g' magnitudes of the Starlink's STARLINK-1130 (Darksat) and 1113 LEO communication satellites and determine the effectiveness of the Darksat darkening treatment at 475.4\,nm. Two observations of the Starlink's Darksat LEO communication satellite were conducted on 2020/02/08 and 2020/03/06 using a Sloan r' and g' filter respectively. While a second satellite, STARLINK-1113 was observed on 2020/03/06 using a Sloan g' filter. The initial observation on 2020/02/08 was a test observation when Darksat was still manoeuvring to its nominal orbit and orientation. Based on the successful test observation, the first main observation was conducted on 2020/03/06 along with an observation of the second Starlink satellite. The calibration, image processing and analysis of the Darksat Sloan g' image gives an estimated Sloan g' magnitude of $7.46\pm0.04$ at a range of 976.50\,km. For STARLINK-1113 an estimated Sloan g' magnitude of $6.59\pm0.05$ at a range of 941.62\,km was found. When scaled to a range of 550\,km and corrected for the solar and observer phase angles, a reduction by a factor of two is seen in the reflected solar flux between Darksat and STARLINK-1113. The data and results presented in this work, show that the special darkening coating used by Starlink for Darksat has darkened the Sloan g' magnitude by $0.77\pm0.05$\,mag, when the range is equal to a nominal orbital height (550\,km). This result will serve members of the astronomical community modelling the satellite mega-constellations, to ascertain their true impact on both the amateur and professional astronomical communities. Concurrent and further observations are planned to cover the full optical and NIR spectrum, from an ensemble of instruments, telescopes and observatories.

We present radiation-magnetohydrodynamic simulations aimed at studying evolutionary properties of H\,{\normalsize II} regions in turbulent, magnetised, and collapsing molecular clouds formed by converging flows in the warm neutral medium. We focus on the structure, dynamics and expansion laws of these regions. Once a massive star forms in our highly structured clouds, its ionising radiation eventually stops the accretion (through filaments) toward the massive star-forming regions. The new over-pressured H\,{\normalsize II} regions push away the dense gas, thus disrupting the more massive collapse centres. Also, because of the complex density structure in the cloud, the H\,{\normalsize II} regions expand in a hybrid manner: they virtually do not expand toward the densest regions (cores), while they expand according to the classical analytical result towards the rest of the cloud, and in an accelerated way, as a blister region, towards the diffuse medium. Thus, the ionised regions grow anisotropically, and the ionising stars generally appear off-centre of the regions. Finally, we find that the hypotheses assumed in standard H\,{\normalsize II}-region expansion models (fully embedded region, blister-type, or expansion in a density gradient) apply simultaneously in different parts of our simulated H\,{\normalsize II} regions, producing a net expansion law ($R \propto t^\alpha$, with $\alpha$ in the range of 0.93-1.47 and a mean value of $1.2 \pm 0.17$) that differs from any of those of the standard models.

We present a method for automatic detection and classification of galaxies which includes a novel data-augmentation procedure to make trained models more robust against the data taken from different instruments and contrast-stretching functions. This method is shown as part of AstroCV, a growing open source computer vision repository for processing and analyzing big astronomical datasets, including high performance Python and C++ algorithms used in the areas of image processing and computer vision. The underlying models were trained using convolutional neural networks and deep learning techniques, which provide better results than methods based on manual feature engineering and SVMs in most of the cases where training datasets are large. The detection and classification methods were trained end-to-end using public datasets such as the Sloan Digital Sky Survey (SDSS), the Galaxy Zoo, and private datasets such as the Next Generation Virgo (NGVS) and Fornax (NGFS) surveys. Training results are strongly bound to the conversion method from raw FITS data for each band into a 3-channel color image. Therefore, we propose data augmentation for the training using 5 conversion methods. This greatly improves the overall galaxy detection and classification for images produced from different instruments, bands and data reduction procedures. The detection and classification methods were trained using the deep learning framework DARKNET and the real-time object detection system YOLO. These methods are implemented in C language and CUDA platform, and makes intensive use of graphical processing units (GPU). Using a single high-end Nvidia GPU card, it can process a SDSS image in 50 milliseconds and a DECam image in less than 3 seconds.

Nanoparticles are ubiquitous in nature and are increasingly important for technology. They are subject to bombardment by ionizing radiation in a diverse range of environments. In particular, nanodiamonds represent a variety of nanoparticles of significant fundamental and applied interest. Here we present a combined experimental and computational study of the behaviour of nanodiamonds under irradiation by xenon ions. Unexpectedly, we observed a pronounced size effect on the radiation resistance of the nanodiamonds: particles larger than 8 nm behave similarly to macroscopic diamond (i.e. characterized by high radiation resistance) whereas smaller particles can be completely destroyed by a single impact from an ion in a defined energy range. This latter observation is explained by extreme heating of the nanodiamonds by the penetrating ion. The obtained results are not limited to nanodiamonds, making them of interest for several fields, putting constraints on processes for the controlled modification of nanodiamonds, on the survival of dust in astrophysical environments, and on the behaviour of actinides released from nuclear waste into the environment.

We present hydrodynamical simulations of interacting Coronal Mass Ejections in the Interplanetary medium (ICMEs). In these events, two consecutive CMEs are launched from the Sun in similar directions within an interval of time of a few hours. In our numerical model, we assume that the ambient solar wind is characterized by its velocity and mass-loss rate. Then, the CMEs are generated when the flow velocity and mass-loss rate suddenly change, with respect to the ambient solar wind conditions during two intervals of time, which correspond to the durations of the CMEs. After their interaction, a merged region is formed and evolve as a single structure into the interplanetary medium. In this work, we are interested in the general morphology of this merged region, which depends on the initial parameters of the ambient solar wind and each of the CMEs involved. In order to understand this morphology, we have performed a parametric study in which we characterize the effects of the initial parameters variations on the density and velocity profiles at 1 AU, using as reference the well-documented event of July 25th, 2004. Based on this parametrization we were able to reproduce with a high accuracy the observed profiles. Then, we apply the parametrization results to the interaction events of May 23, 2010; August 1, 2010; and November 9, 2012. With this approach and using values for the input parameters within the CME observational errors, our simulated profiles reproduce the main features observed at 1 AU. Even though we do not take into account the magnetic field, our models give a physical insight into the propagation and interaction of ICMEs.

We study the correlation between the specific star formation rate of central galaxies and neighbour galaxies, also known as 'galactic conformity', out to 20 Mpc/h using three semi-analytic models (SAMs, one from L-GALAXIES and other two from GALFORM). The aim is to establish whether SAMs are able to show galactic conformity using different models and selection criteria. In all the models, when the selection of primary galaxies is based on an isolation criterion in real space, the mean fraction of quenched galaxies around quenched primary galaxies is higher than that around star-forming primary galaxies of the same stellar mass. The overall signal of conformity decreases when we remove satellites selected as primary galaxies, but the effect is much stronger in GALFORM models compared with the L-GALAXIES model. We find this difference is partially explained by the fact that in GALFORM once a galaxy becomes a satellite remains as such, whereas satellites can become centrals at a later time in L-GALAXIES. The signal of conformity decreases down to 60% in the L-GALAXIES model after removing central galaxies that were ejected from their host halo in the past. Galactic conformity is also influenced by primary galaxies at fixed stellar mass that reside in dark matter haloes of different masses. Finally, we explore a proxy of conformity between distinct haloes. In this case the conformity is weak beyond ~ 3 Mpc/h (<3% in L-GALAXIES, <1-2% in GALFORM models). Therefore, it seems difficult that conformity is directly related with a long-range effect.

We study the stellar populations of bulges of Milky Way-like (MW-like) galaxies with the aim of identifying the physical processes involved in the formation of the bulge of our Galaxy. We use the semi-analytic model of galaxy formation and evolution SAG adapted to this aim; this kind of models can trace the properties of galaxies and their components like stellar discs, bulges and halos, but resolution limits prevent them from reaching the scale of stellar populations (SPs). Properties of groups of stars formed during single star formation events are stored and tracked in the model and results are compared with observations of stars in the galactic bulge. MW-like galaxies are selected using two different criteria. One of them considers intrinsic photo-metric properties and the second is focused on the cosmological context of the local group of galaxies (LG). We compare our model results with spectroscopic and photometric stellar metallicity distributions. We find that 87% of stars in bulges of MWtype galaxies in our model are accreted and formed in starbursts during disc instability events. Mergers contribute to 13% of the mass budget of the bulge and are responsible for the low metallicity tail of the distribution. Abundance ratios of {\alpha} elements with respect to iron, [{\alpha}/Fe], are measured in SPs of model galaxies. The patterns found in the model for SPs with different origins help to explain the lack of a gradient of [{\alpha}/Fe] ratios in observed stars along the minor axis of the bulge.

We observed with the JVLA at 3.6 and 1.3 cm a sample of 11 proto-brown dwarf candidates in Taurus in a search for thermal radio jets driven by the most embedded brown dwarfs. We detected for the first time four thermal radio jets in proto-brown dwarf candidates. We compiled data from UKIDSS, 2MASS, Spitzer, WISE and Herschel to build the Spectral Energy Distribution (SED) of the objects in our sample, which are similar to typical Class~I SEDs of Young Stellar Objects (YSOs). The four proto-brown dwarf candidates driving thermal radio jets also roughly follow the well-known trend of centimeter luminosity against bolometric luminosity determined for YSOs, assuming they belong to Taurus, although they present some excess of radio emission compared to the known relation for YSOs. Nonetheless, we are able to reproduce the flux densities of the radio jets modeling the centimeter emission of the thermal radio jets using the same type of models applied to YSOs, but with corresponding smaller stellar wind velocities and mass-loss rates, and exploring different possible geometries of the wind or outflow from the star. Moreover, we also find that the modeled mass outflow rates for the bolometric luminosities of our objects agree reasonably well with the trends found between the mass outflow rates and bolometric luminosities of YSOs, which indicates that, despite the "excess" centimeter emission, the intrinsic properties of proto-brown dwarfs are consistent with a continuation of those of very low mass stars to a lower mass range. Overall, our study favors the formation of brown dwarfs as a scaled-down version of low-mass stars.

In order to elucidate the origin of spin in both dark matter and baryons in galaxies, we have performed hydrodynamical simulations from cosmological initial conditions. We study atomic cooling haloes in the redshift range $100 > z > 9$ with masses of order $10^9{\rm M_{\odot}}$ at redshift $z=10$. We assume that the gas has primordial composition and that ${\rm H_2}$-cooling and prior star-formation in the haloes have been suppressed. We present a comprehensive analysis of the gas and dark matter properties of four halos with very low ($\lambda \approx 0.01$), low ($\lambda \approx 0.04$), high ($\lambda \approx 0.06$) and very high ($\lambda \approx 0.1$) spin parameter. Our main conclusion is that the spin orientation and magnitude is initially well described by tidal torque linear theory, but later on is determined by the merging and accretion history of each halo. We provide evidence that the topology of the merging region, i.e. the number of colliding filaments, gives an accurate prediction for the spin of dark matter and gas: halos at the center of knots will have low spin while those in the center of filaments will have high spin. The spin of a halo is given by $\lambda \approx 0.05 \times \left(\frac{7.6}{\rm number\,\,\, of \,\,\, filaments}\right)^{5.1}$.

(Abridged). We present numerical simulations of isothermal, MHD, supersonic turbulence, designed to test various hypotheses frequently assumed in star formation(SF) theories. We consider three simulations, each with a different combination of physical size, rms sonic Mach number, and Jeans parameter, but chosen as to give the same value of the virial parameter and to conform with Larson's scaling relations. As in the non-magnetic case: we find no simultaneously subsonic and super-Jeans structures in our MHD simulations. We find that the fraction of small-scale super-Jeans structures increases when self gravity is turned on, and that the production of gravitationally unstable dense cores by turbulence alone is very low. This implies that self-gravity is in general necessary not only to induce the collapse of Jeans-unstable cores, but also to form them. We find that denser regions tend to have more negative values of the velocity divergence, implying a net inwards flow towards the regions' centers. We compare the results from our simulations with the predictions from the recent SF theories by Krumholz & McKee, Padoan & Nordlund, and Hennebelle & Chabrier, using the expressions recently provided by Federrath & Klessen. We find that none of these theories reproduces the dependence of the SFEff with Ms observed in our simulations in the MHD case. The SFEff predicted by the theories ranges between half and one order of magnitude larger than what we observe in the simulations in both the HD and the MHD cases. We conclude that the type of flow used in simulations like the ones presented here and assumed in recent SF theories, may not correctly represent the flow within actual clouds, and that theories that assume it does may be missing a fundamental aspect of the flow. We suggest that a more realistic regime may be that of hierarchical gravitational collapse.

We present a formalism of the dynamics of internal shocks in relativistic jets where the source has a time-dependent injection velocity and mass-loss rate. The variation of the injection velocity produces a two-shock wave structure, the working surface, that moves along the jet. This new formalism takes into account the fact that momentum conservation is not valid for relativistic flows where the relativistic mass lost by radiation must be taken into account, in contrast to the classic regime. We find analytic solutions for the working surface velocity and radiated energy for the particular case of a step function variability of the injection parameters. We model two cases: a pulse of fast material and a pulse of slow material (with respect to the mean flow). Applying these models to gamma ray burst light curves, one can determine the ratio of the Lorentz factors gamma_2 / gamma_1 and the ratio of the mass-loss rates dot{m_2} / dot{m_1} of the upstream and downstream flows. As an example, we apply this model to the sources GRB 080413B and GRB 070318 and find the values of these ratios. Assuming a Lorentz factor gamma_1=100, we further estimate jet mass-loss rates between dot{m_1} ~ 10^{-5}-1 Msun.yr^{-1}. We also calculate the fraction of the injected mass lost by radiation. For GRB 070318 this fraction is ~7%. In contrast, for GRB 080413B this fraction is larger than 50%; in this case radiation losses clearly affect the dynamics of the internal shocks.

We measure the evolution of galaxy clustering out to a redshift of z~1.5 using data from two MUSYC fields, the Extended Hubble Deep Field South (EHDF-S) and the Extended Chandra Deep Field South (ECDF-S). We use photometric redshift information to calculate the projected-angular correlation function, omega(sigma), from which we infer the projected correlation function Xi(sigma). We demonstrate that this technique delivers accurate measurements of clustering even when large redshift measurement errors affect the data. To this aim we use two mock MUSYC fields extracted from a LambdaCDM simulation populated with GALFORM semi-analytic galaxies which allow us to assess the degree of accuracy of our estimates of Xi(sigma) and to identify and correct for systematic effects in our measurements. We study the evolution of clustering for volume limited subsamples of galaxies selected using their photometric redshifts and rest-frame r-band absolute magnitudes. We find that the real-space correlation length r_0 of bright galaxies, M_r<-21 (rest-frame) can be accurately recovered out to z~1.5, particularly for ECDF-S given its near-infrared photometric coverage. There is mild evidence for a luminosity dependent clustering in both fields at the low redshift samples (up to <z>=0.57), where the correlation length is higher for brighter galaxies by up to 1Mpc/h between median rest-frame r-band absolute magnitudes of -18 to -21.5. As a result of the photometric redshift measurement, each galaxy is assigned a best-fit template; we restrict to E and E+20%Sbc types to construct subsamples of early type galaxies (ETGs). Our ETG samples show a strong increase in r_0 as the redshift increases, making it unlikely (95% level) that ETGs at median redshift z_med=1.15 are the direct progenitors of ETGs at z_med=0.37 with equivalent passively evolved luminosities. (ABRIDGED)

We present here new results of two-dimensional hydrodynamical simulations of the eruptive events of the 1840s (the great) and the 1890s (the minor) eruptions suffered by the massive star $\eta$ Car. The two bipolar nebulae commonly known as the Homunculus and the little Homunculus were formed from the interaction of these eruptive events with the underlying stellar wind. As in previous work (Gonzalez et al. 2004a, 2004b), we assume here an interacting, nonspherical multiple-phase wind scenario to explain the shape and the kinematics of both Homunculi, but adopt a more realistic parametrization of the phases of the wind. During the 1890s eruptive event, the outflow speed {\it decreased} for a short period of time. This fact suggests that the little Homunculus is formed when the eruption ends, from the impact of the post-outburst $\eta$ Car wind (that follows the 1890s event) with the eruptive flow (rather than by the collision of the eruptive flow with the pre-outburst wind, as claimed in previous models; Gonzalez et al. 2004a, 2004b). Our simulations reproduce quite well the shape and the observed expansion speed of the large Homunculus. The little Homunculus (which is embedded within the large Homunculus) becomes Rayleigh-Taylor unstable and develop filamentary structures that resembles the spatial features observed in the polar caps. In addition, we find that the interior cavity between the two Homunculi is partially filled by material that is expelled during the decades following the great eruption. This result may be connected with the observed double-shell structure in the polar lobes of the $\eta$ Car nebula. Finally, as in previous work, we find the formation of tenuous, equatorial, high-speed features that seem to be related to the observed equatorial skirt of $\eta$ Car.

The opacity of a spiral disk due to dust absorption influences every measurement we make of it in the UV and optical. Two separate techniques directly measure the total absorption by dust in the disk: calibrated distant galaxy counts and overlapping galaxy pairs. The main results from both so far are a semi-transparent disk with more opaque arms, and a relation between surface brightness and disk opacity. In the Spitzer era, SED models of spiral disks add a new perspective on the role of dust in spiral disks. Combined with the overall opacity from galaxy counts, they yield a typical optical depth of the dusty ISM clouds: 0.4 that implies a size of ~60 pc. Work on galaxy counts is currently ongoing on the ACS fields of M51, M101 and M81. Occulting galaxies offer the possibility of probing the history of disk opacity from higher redshift pairs. Evolution in disk opacity could influence distance measurements (SN1a, Tully-Fisher relation). Here, we present first results from spectroscopically selected occulting pairs in the SDSS. The redshift range for this sample is limited, but does offer a first insight into disk opacity evolution as well as a reference for higher redshift measurements.

In this paper we explore the relation between dust extinction and stellar light distribution in disks of spiral galaxies. Extinction influences our dynamical and photometric perception of disks, since it can distort our measurement of the contribution of the stellar component. To characterize the total extinction by a foreground disk, Gonzalez et al. (1998) proposed the ``Synthetic Field Method'' (SFM), which uses the calibrated number of distant galaxies seen through the foreground disk as a direct indication of extinction. The method is described in Gonzalez et al. (1998) and Holwerda et al. (2005a). To obtain good statistics, the method was applied to a set of HST/WFPC2 fields Holwerda et al. (2005b) and radial extinction profiles were derived, based on these counts. In the present paper, we explore the relation of opacity with surface brightness or color from 2MASS images, as well as the relation between the scalelengths for extinction and light in the I band. We find that there is indeed a relation between the opacity (A_I) and the surface brightness, particularly at the higher surface brightnesses. No strong relation between near infrared (H-J, H-K) color and opacity is found. The scalelengths of the extinction are uncertain for individual galaxies but seem to indicate that the dust distribution is much more extended than the stellar light. The results from the distant galaxy counts are also compared to the reddening derived from the Cepheids light-curves Freedman et al. (2001). The extinction values are consistent, provided the selection effect against Cepheids with higher values of A_I is taken into account. The implications from these relations for disk photometry, M/L conversion and galaxy dynamical modeling are briefly discussed.

The opacity of spiral galaxy disks, from counts of distant galaxies, is compared to HI column densities. The opacity measurements are calibrated using the ``Synthetic Field Method'' from Gonzalez et al (1998) and Holwerda et al. (2005a). When compared for individual disks, the HI column density and dust opacity do not seem to be correlated as HI and opacity follow different radial profiles. To improve statistics, an average radial opacity profile is compared to an average HI profile. Compared to dust-to-HI estimates from the literature, more extinction is found in this profile. This difference may be accounted for by an underestimate of the dust in earlier measurements due to their dependence on dust temperature. Since the SFM is insensitive to the dust temperature, the ratio between the SFM opacity and HI could very well be indicative of the true ratio. Earlier claims for a radially extended cold dust disk were based on sub-mm observations. A comparison between sub-mm observations and counts of distant galaxies is therefore desirable. We present the best current example of such a comparison, M51, for which the measurements seem to agree. However, this remains an area where improved counts of distant galaxies, sub-mm observations and our understanding of dust emissivity are needed.

We here present two-dimensional, time-dependent radiatively cooling hydrodynamical simulations of the large and little Homunculus nebulae around $\eta$ Carinae. We employ an alternative scenario to previous interacting stellar wind models which is supported by both theoretical and observational evidence, where a non-spherical outburst wind (with a latitudinal velocity dependence that matches the observations of the large Homunculus), which is expelled for 20 years, interacts with a pre-eruptive slow wind also with a toroidal density distribution, but with a much smaller equator-to-polar density contrast than that assumed in previous models. A second eruptive wind with spherical shape is ejected about 50 years after the first outburst, and causes the development of the little internal nebula. We find that, as a result of an appropriate combination of the parameters that control the degree of asymmetry of the interacting winds, we are able to produce not only the structure and kinematics of both Homunculus, but also the high-velocity equatorial ejecta. These arise from the impact between the non-spherical outburst and the pre-outburst winds in the equatorial plane.