Photometric Model Research Articles

Light curve modeling of the eclipsing binary systems V0876 Lyr, V3660 Oph, and V0988 Mon

Abstract Photometric study and light curve modeling for the systems V0876 Lyr, V3660 Oph, and V0988 Mon were carried out by means of their first charge-coupled device observations and the Wilson–Devinney code based on model atmospheres by Kurucz (1993. In: Milone, E. (Ed.), Light curve modeling of eclipsing binary stars. New York: Springer-Verlag. p. 93). The accepted models reveal physical parameters that show that all the primary components of the studied systems are more massive and hotter than the secondary ones. The systems V0876 Lyr and V0988 Mon were classified as over contact, with mass ratios of 0.8940 and 0.3904, respectively, while the system V366Oph was classified as detached, with a mass ratio of 0.1181. The evolutionary state of the system components was investigated, and their spectral types were adopted. Both components of the system V0876 Lyr and secondary components of the systems V3660 Oph and V0988 Mon show lower radii and luminosities than expected for Zero Age Main Sequence (ZAMS). Secondary components of the system, V3660 Oph and V0988 Mon, deviated from ZAMS and the Terminal-age Main Sequence (TAMS) with a higher radius and luminosity than expected for ZAMS and TAMS.

Prior-shape-guided photometric stereo model for 3D damage measurement of worn surfaces

Three-dimensional (3D) damage measurement often encounters limited accuracy due to the non-Lambertian effects on worn surfaces. To address this issue, a prior-shape-guided surface reconstruction model is introduced with photometric stereo vision. Initially, an illumination calibration method is constructed to mitigate the uneven brightness on worn surface images. Subsequently, a prior-shape-guided network is developed that employs the feature aggregation strategy and prior shape uncertainty to guide worn surface reconstruction. Experimental results reveal that the proposed model can effectively recover the typical morphologies of worn surfaces, and the root mean square error is reduced from 7.39 to 1.59 when compared to existing methods. Most importantly, this model has been applied to detect aero-engine bearings and achieve 3D quantitative characterization for worn damages.

A Contact Binary Misclassified as an Ellipsoidal Variable: Complications for Detached Black Hole Searches

Identifying sources exhibiting ellipsoidal variability in large photometric surveys is becoming a promising method to search for candidate detached black holes (BHs) in binaries. This technique aims to exploit the orbital-phase-dependent modulation in optical photometry caused by the BH distorting the shape of the luminous star to constrain the mass ratio of the binary. Without understanding if, or how much, contamination is present in the candidate BH samples produced by this new technique it is hard to leverage them for BH discovery. Here, we follow up one of the best candidates identified from Gaia Data Release 3, Gaia DR3 4042390512917208960, with a radial velocity (RV) campaign. Combined photometric and RV modeling, along with spectral disentangling, suggests that the true mass ratio (the mass of the unseen object divided by the mass of the luminous star) is an order of magnitude smaller than that inferred assuming the modulations arise from ellipsoidal variability. We therefore infer that this system is likely a contact binary, or on the boundary of both stars nearly filling their Roche lobes; however, further observations are required to confidently detect the secondary. We find that the well-known problem of discriminating between ellipsoidal and contact binary light curves results in a larger contamination from contact binaries than previously suggested. Until ellipsoidal variables can be reliably distinguished from contact binaries, samples of BH candidates selected based on ellipsoidal variability are likely to be highly contaminated by contact binaries or similar systems.

The MAGPI Survey: massive slow rotator population in place by z ∼ 0.3

ABSTRACT We use the ‘Middle Ages Galaxy Properties with Integral field spectroscopy’ (MAGPI) survey to investigate whether galaxies have evolved in the distribution of their stellar angular momentum in the past 3–4 Gyr, as probed by the observational proxy for spin, λR. We use 2D stellar kinematics to measure λR along with detailed photometric models to estimate galaxy ellipticity. The combination of these measurements quantifies the kinematic classes of ‘fast rotators’ and the rarer ‘slow rotators’, which show no regular rotation in their line-of-sight velocity fields. We compare 51 MAGPI galaxies with log10 (M⋆/M⊙) > 10 to carefully drawn samples of MaNGA galaxies in the local Universe, selected to represent possible descendants of the MAGPI progenitors. The EAGLE simulations are used to identify possible evolutionary pathways between the two samples, explicitly accounting for progenitor bias in our results and the varied evolutionary pathways a galaxy might take between the two epochs. We find that the occurrence of slow rotating galaxies is unchanged between the MAGPI (z ∼ 0.3) and MaNGA (z ∼ 0) samples, suggesting the massive slow rotator population was already in place ∼4 Gyr ago and has not accumulated since. There is a hint of the MAGPI sample having an excess of high λR galaxies compared to the MaNGA sample, corresponding to more ordered rotation, but statistically the samples are not significantly different. The large-scale stellar kinematics, as quantified through the λR parameter, of galaxies at z ∼ 0.3 have already evolved into the diversity of structures seen today in the local Universe.

Phobos photometric properties from Mars Express HRSC observations

Aims. This study aims to analyze Phobos’ photometric properties using Mars Express mission observations to support the Martian Moons exploration mission (MMX) devoted to the investigation of the Martian system and to the return of Phobos samples. Methods. We analyzed resolved images of Phobos acquired between 2004 and 2022 by the High Resolution Stereo Camera (HRSC) on board the Mars Express spacecraft at a resolution ranging from ~30 m px−1 to 330 m px−1. We used data acquired with the blue, green, red, and IR filters of HRSC and the panchromatic data of the Super Resolution Channel (SRC). The SRC data are unique because they cover small phase angles (0.2–10°), permitting the investigation of the Phobos opposition effect. We simulated illumination and geometric conditions for the different observations using the Marx Express and the camera spice kernels provided by the HRSC team. We performed photometric analysis using the Hapke model for both integrated and disk-resolved data. Results. The Phobos phase function is characterized by a strong opposition effect due to shadow hiding, with an amplitude and a half-width of the opposition surge of 2.28±0.03 and 0.0573±0.0001, respectively. Overall, the surface of Phobos is dark, with a geometric albedo of 6.8% in the green filter and backscattering. Its single-scattering albedo (SSA) value (7.2% in the green filter) is much higher than what has been found for primitive asteroids and cometary nuclei and is close to the values reported in the literature for Ceres. We also found a surface porosity of 87%, indicating the presence of a thick dust mantle or of fractal aggregates on the top surface. The SSA maps revealed high reflectance variability, with the blue unit area in the northeast Stickney rim being up to 65% brighter than average, while the Stickney floor is among the darkest regions, with reflectance 10 to 20% lower than average. Photometric modeling of the regions of interest selected in the red and blue units indicates that red unit terrains have a stronger opposition effect and a smaller SSA value than the blue ones, but they have similar porosity and backscattering properties. Conclusions. The HRSC data provide a unique investigation of the Phobos phase function and opposition surge, which is valuable information for the MMX observational planning. The Phobos opposition surge, surface porosity, phase integral, and spectral slope are very similar to the values observed for the comet 67P and for Jupiter family comets in general. Based on these similarities, we formulate a hypothesis that the Mars satellites might be the results of a binary or bilobated comet captured by Mars.

The temporal variability of Io’s hotspots

We investigate the variability of the power emission of Io’s hotspots by using recent Juno/JIRAM infrared observations. The Jovian Infrared Auroral Mapper (JIRAM) is an imaging spectrometer which began observing Jupiter in August 2016. Although observing Jupiter’s moons is not its primary objective, JIRAM can use the frequent opportunities to observe Io (up to once per orbit) to gather infrared images and spectra of its surface. The present study uses the data acquired by JIRAM during the last 2 years, including the location and morphology of Io’s hotspots, and the temporal variability of the total output. A new photometric model for the hotspots and the dayside surface has been developed, which permits us to disentangle the temporal variability from the changes in the observation geometry. While the latitudinal dependence of the power output is not well constrained, low-latitude hotspots show a significantly more intense temporal variability and greater temperature.

BRDF-Based Photometric Modeling of LEO Constellation Satellite from Massive Observations

BRDF-Based Photometric Modeling of LEO Constellation Satellite from Massive Observations

Multicolour photometry of LEO mega-constellations Starlink and OneWeb

ABSTRACT The development of low earth orbit (LEO) mega-constellation fundamentally threatens ground-based optical astronomical observations. To study the photometric properties of the LEO mega-constellations, we used the Xinglong 50 cm telescope to conduct a large-sample, high-precision, and multicolour target-tracking photometry of two typical LEO mega-constellations: Starlink and OneWeb. Over a three-month observation period starting on 2022 January 1, we collected 1447 light curves of 404 satellites in four typical versions: Starlink v1.0, DarkSat, VisorSat, Starlink v1.5, and OneWeb. According to data statistics, Starlink v1.0 has the smallest median magnitude at clear and Sloan Digital Sky Survey gri band, and OneWeb is the dimmest bus. The brightness of Starlink v1.5 is slightly brighter than VisorSat. We construct a detailed photometric model with solar phase angle variations by calculating the illumination-visibility geometry based on the orbital parameters. Our data analysis shows that the solar phase angle is the significant characteristic which influencing Starlink satellites’ brightness, but it is not sensitive to OneWeb satellites. VisorSat and Starlink v1.5 versions, which are equipped with deployable visors, have significantly reduced scattered light compared to the previous Starlink v1.0 version. The multiband LOWESS and colour index are analysed in characterizing the energy and colour features of LEO mega-constellation satellites. This work found that the proportion of scattered sunlight mitigation achieved with VisorSat and Starlink v1.5 was 55.1 and 40.4 per cent, respectively. The colour index of different buses shows an evident clustering feature. Our observation and analysis could provide valuable quantitative data and photometric models, which can contribute to assessing the impact of LEO mega-constellations on astronomical observations.

Sample data selection based on the feature subspace for photometric correction of the Chang'E-4 in situ observations

The lunar surface exhibits noticeable variations in photometric properties across different regions, influenced by their specific compositional and physical characteristics. These variations are particularly pronounced at the millimeter to centimeter scale. In this study, we propose a feature subspace-based sample selection method to identify Chang'E-4 Visible Near-infrared Imaging Spectrometer (CE-4 VNIS) in situ observations that share similar photometric properties at the centimeter scale. By constructing a feature subspace through principal component transformation, the CE-4 VNIS observation sites with comparable surface photometric properties are able to be selected from a series of observations. The selection results demonstrate that the majority of observation sites are covered by the lunar regolith, with a few exceptions consisting of rocks, soil bulk, shadows, and rough surfaces. These observations should be excluded from the photometric correction process during phase function fitting. Analyzing the photometric performance of the VNIS data with lunar regolith reveals that reflectance decreases as the phase angle increases below 80°, while above 80°, reflectance increases with increasing phase angle. This phenomenon is likely attributed to stronger forward scattering. Thus, a second-order polynomial is employed to fit the phase function for VNIS data. Consequently, a photometric model is developed for lunar regolith observed by CE-4 VNIS, incorporating the fitted phase function. Experimental results demonstrate reduced differences among observations with varying phase angles, and the corrected VNIS spectra exhibit a phase-reddening effect. These findings validate the effectiveness of the proposed model. A comparison between the corrected data with and without sample selection reveals a 33 % decrease in mean standard deviation for VIS/NIR bands and a 19 % decrease for SWIR bands in the corrected spectra after sample selection. These results indicate the potential applicability of our method for future VNIS observations, enabling the acquisition of photometrically consistent data.

The first CCD photometric studies of the member eclipsing binary ZTFJ015003.88+534734.1 in the newly discovered young open cluster UBC 188

We present the first CCD observations of an eclipsing binary, ZTFJ015003.88+534734.1, which is a member in the open star cluster UBC 188. The observations were taken by the 1.88 m telescope at the Kottamia Astronomical Observatory (KAO) in SDSS griz bands. The latest version of the Wilson-Devinney (W-D) code was employed for photometric analysis and light curve modeling of the eclipsing binary. The results indicate that the binary system is in an over-contact configuration. The mass of the primary star (M1) is determined to be 1.293 M⊙, and the mass of the secondary star (M2) is directly derived from the system's estimated mass ratio (q = M2/M1) as 0.340 times the solar mass (M⊙). We investigated the color-magnitude diagram and the membership probability of the open cluster UBC 188 using the Gaia DR3 data. We determined the membership probability of the eclipsing binary ZTFJ015003.88+534734.1 using the pyUPMASK algorithm and found that its membership probability is one.

Spectroscopic r-process Abundance Retrieval for Kilonovae. II. Lanthanides in the Inferred Abundance Patterns of Multicomponent Ejecta from the GW170817 Kilonova

In kilonovae, freshly synthesized r-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the r-process compositions of the merger ejecta is computationally challenging. Vieira et al. introduced Spectroscopic r-process Abundance Retrieval for Kilonovae (SPARK), a software tool to infer elemental abundance patterns of the ejecta and associate spectral features with particular species. Previously, we applied SPARK to the 1.4-day spectrum of AT2017gfo and inferred its abundance pattern for the first time, characterized by electron fraction Y e = 0.31, a substantial abundance of strontium, and a dearth of lanthanides and heavier elements. This ejecta is consistent with wind from a remnant hypermassive neutron star and/or accretion disk. We now extend our inference to spectra at 2.4 and 3.4 days and test the need for multicomponent ejecta, where we stratify the ejecta in composition. The ejecta at 1.4 and 2.4 days is described by the same single blue component. At 3.4 days, a new redder component with lower Y e = 0.16 and a significant abundance of lanthanides emerges. This new redder component is consistent with dynamical ejecta and/or neutron-rich ejecta from a magnetized accretion disk. As expected from photometric modeling, this component emerges as the ejecta expands, the photosphere recedes, and the earlier bluer component dims. At 3.4 days, we find an ensemble of lanthanides, with the presence of cerium most concrete. This presence of lanthanides has important implications for the contribution of kilonovae to the r-process abundances observed in the Universe.

Photometry of the Didymos System across the DART Impact Apparition

On 2022 September 26, the Double Asteroid Redirection Test (DART) spacecraft impacted Dimorphos, the satellite of binary near-Earth asteroid (65803) Didymos. This demonstrated the efficacy of a kinetic impactor for planetary defense by changing the orbital period of Dimorphos by 33 minutes. Measuring the period change relied heavily on a coordinated campaign of lightcurve photometry designed to detect mutual events (occultations and eclipses) as a direct probe of the satellite’s orbital period. A total of 28 telescopes contributed 224 individual lightcurves during the impact apparition from 2022 July to 2023 February. We focus here on decomposable lightcurves, i.e., those from which mutual events could be extracted. We describe our process of lightcurve decomposition and use that to release the full data set for future analysis. We leverage these data to place constraints on the postimpact evolution of ejecta. The measured depths of mutual events relative to models showed that the ejecta became optically thin within the first ∼1 day after impact and then faded with a decay time of about 25 days. The bulk magnitude of the system showed that ejecta no longer contributed measurable brightness enhancement after about 20 days postimpact. This bulk photometric behavior was not well represented by an HG photometric model. An HG 1 G 2 model did fit the data well across a wide range of phase angles. Lastly, we note the presence of an ejecta tail through at least 2023 March. Its persistence implied ongoing escape of ejecta from the system many months after DART impact.

Galaxies decomposition with spiral arms – II. A multiwavelength case study of M 51

ABSTRACT Spiral structure can contribute significantly to a galaxy’s luminosity. However, only rarely are proper photometric models of spiral arms used in decompositions. As we show in the previous work, including the spirals as a separate component in a photometric model of a galaxy would both allow to obtain their structural parameters, and reduce the systematic errors in estimating the parameters of other components. Doing so in different wavebands, one can explore how their properties vary with the wavelength. In this paper, second in this series, we perform decomposition of M 51 in 17 bands, from the far-ultraviolet (UV) to far-infrared, using imaging from the DustPedia project. We use the same 2D photometric model of spiral structure where each arm is modelled independently. The complex and asymmetric spiral structure in M 51 is reproduced relatively well with our model. We analyse the differences between models with and without spiral arms, and investigate how the fit parameters change with wavelength. In particular, we find that the spiral arms demonstrate the largest width in the optical, whereas their contribution to the galaxy luminosity is most significant in the UV. The disc central intensity drops by a factor of 1.25–3 and its exponential scale changes by 5–10 per cent when spiral arms are included, depending on wavelength. Taking into account the full light distribution across the arms, we do not observe the signs of a long-lived density wave in the spiral pattern of M 51 as a whole.

Characterization of the Eclipsing Binary System EPIC 211982753

We present the physical parameters of an eclipsing binary system EPIC 211982753 derived through photometric and radial velocity data modeling. We make use of photometric data from NASA’s K2 mission, ASAS-SN, and 1.3-m Devasthal Fast Optical Telescope (DFOT) while spectroscopic data have been acquired from the HERMES spectrograph at 1.2-m Mercator telescope. The linear ephemeris for the system is updated using the K2 mission data. The synthetic lightcurve and radial velocity curves are generated with the help of eclipsing binary modeling package PHOEBE 1.0. The masses of primary and secondary components are determined as 1.64 ±0.02 and 1.55 ±0.01 M⊙, respectively. The radius for primary and secondary components are estimated as 1.73 ±0.02 and 1.47 ±0.02 R⊙, respectively. The distance of the system is calculated as 238 ± 4 pc. The eclipsing binary is found to be a total eclipsing system with a high mass ratio of q=0.94.

Deep Learning Methods for Calibrated Photometric Stereo and Beyond.

Photometric stereo recovers the surface normals of an object from multiple images with varying shading cues, i.e., modeling the relationship between surface orientation and intensity at each pixel. Photometric stereo prevails in superior per-pixel resolution and fine reconstruction details. However, it is a complicated problem because of the non-linear relationship caused by non-Lambertian surface reflectance. Recently, various deep learning methods have shown a powerful ability in the context of photometric stereo against non-Lambertian surfaces. This paper provides a comprehensive review of existing deep learning-based calibrated photometric stereo methods utilizing orthographic cameras and directional light sources. We first analyze these methods from different perspectives, including input processing, supervision, and network architecture. We summarize the performance of deep learning photometric stereo models on the most widely-used benchmark data set. This demonstrates the advanced performance of deep learning-based photometric stereo methods. Finally, we give suggestions and propose future research trends based on the limitations of existing models.

Discovery of a Dust Sorting Process on Boulders Near the Reiner Gamma Swirl on the Moon

AbstractIn a database of lunar fractured boulders (Rüsch & Bickel, 2023, https://doi.org/10.3847/psj/acd1ef), we found boulders with reflectance features dissimilar to previously known morphologies. We performed a photo‐geologic investigation and determined that the features correspond to a dust mantling on top of boulders with a unique photometric behavior. We next performed a photometric model inversion on the dust mantling using Bayesian inference sampling. Modeling indicates that the dust photometric anomaly is most likely due to a reduced opposition effect, whereas the single scattering albedo is not significantly different from that of the nearby background regolith. This implies a different structure of the dust mantling relative to the normal regolith. We identified and discussed several potential processes to explain the development of such soil. None of these mechanisms can entirely explain the multitude of observational constraints unless evoking anomalous boulder properties. Further study of these boulders can shed light on the workings of a natural dust sorting process potentially involving dust dynamics, a magnetic field, and electrostatic dust transport. The presence of these boulders appears to be limited to the Reiner K crater near the Reiner Gamma magnetic and photometric anomaly. This close spatial relationship further highlights that poorly understood processes occur in this specific region of the Moon.

Spectrophotometric properties of CoPhyLab’s dust mixtures

ABSTRACT Objective: In the framework of the Cometary Physics Laboratory (CoPhyLab) and its sublimation experiments of cometary surface analogues under simulated space conditions, we characterize the properties of intimate mixtures of juniper charcoal and SiO2 chosen as a dust analogue (Lethuillier et al. 2022). We present the details of these investigations for the spectrophotometric properties of the samples. Methods: We measured these properties using a hyperspectral imager and a radio-goniometer. From the samples’ spectra, we evaluated reflectance ratios and spectral slopes. From the measured phase curves, we inverted a photometric model for all samples. Complementary characterizations were obtained using a pycnometer, a scanning electron microscope and an organic elemental analyser. Results: We report the first values for the apparent porosity, elemental composition, and VIS-NIR spectrophotometric properties for juniper charcoal, as well as for intimate mixtures of this charcoal with the SiO2. We find that the juniper charcoal drives the spectrophotometric properties of the intimate mixtures and that its strong absorbance is consistent with its elemental composition. We find that SiO2 particles form large and compact agglomerates in every mixture imaged with the electron microscope, and its spectrophotometric properties are affected by such features and their particle-size distribution. We compare our results to the current literature on comets and other small Solar system bodies and find that most of the characterized properties of the dust analogue are comparable to some extent with the spacecraft-visited cometary nuclei, as well as to Centaurs, Trojans, and the bluest TNOs.

Galaxies decomposition with spiral arms – I: 29 galaxies from S4G

ABSTRACT Spiral structure can occupy a significant part of the galaxy, but properly accounting for it in photometric decomposition is rarely done. This may lead to significant errors in the parameters determined. To estimate how exactly neglecting the presence of spiral arms affects the estimation of galaxy decomposition parameters, we perform fitting of 29 galaxies considering spiral arms as a separate component. In this study, we utilize 3.6 μm-band images from the S4G survey and use a new 2D photometric model where each spiral arm is modelled independently. In our model, the light distribution both along and across the arm and its overall shape can be varied significantly. We analyse the differences between models with and without spiral arms, and show that neglecting spiral arms in decomposition causes errors in estimating the parameters of the disc, the bulge, and the bar. We retrieve different parameters of the spiral arms themselves, including their pitch angles, widths, and spiral-to-total luminosity ratio, and examine various relations between them and other galaxy parameters. In particular, we find that the spiral-to-total ratio is higher for galaxies with more luminous discs and with higher bulge-to-total ratios. We report that the pitch angle of spiral arms decreases with increasing bulge or bar fraction. We measure the width of the spiral arms to be 53 per cent of the disc scale length, on average. We examine the contribution of the spiral arms to the azimuthally averaged brightness profile and find that spiral arms produce a ‘bump’ on this profile with a typical height of 0.3–0.7 mag.

Optical Monitoring of the Didymos–Dimorphos Asteroid System with the Danish Telescope around the DART Mission Impact

The NASA’s Double-Asteroid Redirection Test (DART) was a unique planetary defence and technology test mission, the first of its kind. The main spacecraft of the DART mission impacted the target asteroid Dimorphos, a small moon orbiting the asteroid Didymos (65803), on 2022 September 26. The impact brought up a mass of ejecta which, together with the direct momentum transfer from the collision, caused an orbital period change of 33 ± 1 minutes, as measured by ground-based observations. We report here the outcome of the optical monitoring campaign of the Didymos system from the Danish 1.54 m telescope at La Silla around the time of impact. The observations contributed to the determination of the changes in the orbital parameters of the Didymos–Dimorphos system, as reported by Thomas et al., but in this paper we focus on the ejecta produced by the DART impact. We present photometric measurements from which we remove the contribution from the Didymos–Dimorphos system using an H–G photometric model. Using two photometric apertures we determine the fading rate of the ejecta to be 0.115 ± 0.003 mag day−1 (in a 2″ aperture) and 0.086 ± 0.003 mag day−1 (5″) over the first week postimpact. After about 8 days postimpact we note the fading slows down to 0.057 ± 0.003 mag day−1 (2″ aperture) and 0.068 ± 0.002 mag day−1 (5″). We include deep-stacked images of the system to illustrate the ejecta evolution during the first 18 days, noting the emergence of dust tails formed from ejecta pushed in the antisolar direction, and measuring the extent of the particles ejected Sunward to be at least 4000 km.

The Advantages of Global Photometric Models in Fitting Transit Variations

Estimation of planetary orbital and physical parameters from light-curve data relies heavily on the accurate interpretation of Transit Timing Variation (TTV) measurements. In this letter, we review the process of TTV measurement and compare two fitting paradigms—one that relies on making transit-by-transit timing estimates and then fitting a TTV model to the observed timings, and one that relies on fitting a global flux model to the entire light-curve data set simultaneously. The latter method is achieved either by solving for the underlying planetary motion (often referred to as “photodynamics”), or by using an approximate or empirical shape of the TTV signal. We show that, across a large range of the transit S/N regime, the probability distribution function of the mid-transit time significantly deviates from a Gaussian, even if the flux errors do distribute normally. Treating the timing uncertainties as if they are distributed normally leads, in such a case, to a wrong interpretation of the TTV measurements. We illustrate these points using numerical experiments and conclude that a fitting process that relies on a global flux fitting, rather than the derived TTVs, should be preferred.