Some orbital characteristics of lunar artificial satellites

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Sub-Optimal Orbital Maneuvers for Artificial Satellites

The determination of a specific orbit and the procedure to calculate orbital maneuvers of artificial satellites are problems of extreme importance in the study of orbital mechanics. Therefore, the transferring problem of a spaceship from one orbit to another, and the attention due to this subject has in increased during the last years. Many applications can be found in several space activities, for example, to put a satellite in a geostationary orbit, to change the position of a spaceship, to maintain a specific satellite's orbit, in the design of an interplanetary mission, and others. The Brazilian Satellite SCD-1 (Data Collecting Satellite) will be used as example in this paper. It is the first satellite developed entirely in Brazil, and it remains in operation to this date. SCD-1 was designed, developed, built, and tested by Brazilian scientists, engineers, and technicians working at INPE (National Institute for Space Research, and in Brazilian Industries. During the lifetime, it might be necessary do some complementary maneuvers, being this one either an orbital transferring, or just to make periodical corrections. The purpose of transferring problem is to change the position, velocity and the satellite's mass to a new pre determined state. This transfer can be totally linked (in the case of "Rendezvous") or partially free (free time, free final velocity, etc). In the global case, the direction, the orientation and the magnitude of the thrust to be applied must be chosen, respecting the equipment's limit. In order to make this transferring, either sub-optimal or optimal maneuvers may be used. In the present study, only the sub-optimal will be shown. Hence, this method will simplify the direction of thrust application, to allow a fast calculation that may be used in real time, with a very fast processing. The thrust application direction to be applied will be assumed small and constant, and the purpose of this paper is to find the time interval that the thrust is applied. This paper is basically divided into three parts: during the first one the sub-optimal maneuver is explained and detailed, the second presents the Satellite SCD-1, and finally the last part shows the results using the sub-optimal maneuver applied to the Brazilian Satellite.

Nebular gas drag and co-orbital system dynamics

Aims. We study trajectories of planetesimals whose orbits decay due to gas drag in a primordial solar nebula and are perturbed by the gravity of the secondary body on an eccentric orbit whose mass ratio takes values from mu(2) = 10(-7) to mu(2) = 10(-3) increasing ten times at each step. Each planetesimal ultimately suffers one of the three possible fates: (1) trapping in a mean motion resonance with the secondary body; (2) collision with the secondary body and consequent increase of its mass; or (3) diffusion after crossing the orbit of the secondary body.Methods. We take the Burlirsh-Stoer numerical algorithm in order to integrate the Newtonian equations of the planar, elliptical restricted three-body problem with the secondary body and the planetesimal orbiting the primary. It is assumed that there is no interaction among planetesimals, and also that the gas does not affect the orbit of the secondary body.Results. The results show that the optimal value of the gas drag constant k for the 1: 1 resonance is between 0.9 and 1.25, representing a meter size planetesimal for each AU of orbital radius. In this study, the conditions of the gas drag are such that in theory, L4 no longer exists in the circular case for a critical value of k that defines a limit size of the planetesimal, but for a secondary body with an eccentricity larger than 0.05 when mu(2) = 10(-6), it reappears. The decrease of the cutoff collision radius increase the difusions but does not affect the distribution of trapping. The contribution to the mass accretion of the secondary body is over 40% with a collision radius 0.05R(Hill) and less than 15% with 0.005R(Hill) for mu(2) = 10(-7). The trappings no longer occur when the drag constant k reachs 30. That means that the size limit of planetesimal trapping is 0.2 m per AU of orbital radius. In most cases, this accretion occurs for a weak gas drag and small secondary eccentricity. The diffusions represent most of the simulations showing that gas drag is an efficient process in scattering planetesimals and that the trapping of planetesimals in the 1: 1 resonance is a less probable fate. These results depend on the specific drag force chosen.

Ação de forças gravitacionais e não gravitacionais sobre o movimento orbital de satélites artificiais

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Pertubação orbital devida a um terceiro corpo com distribuição não uniforme de massa e em órbita elíptica

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Application of effective theory and Grassmann variables to model compact binaries with spin

Pós-graduação em Física - IFT

Capacitance spectroscopy and density functional theory

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Geocenter coordinates estimated from GNSS data as viewed by perturbation theory

Time series of geocenter coordinates were determined with data of two global navigation satellite systems (GNSSs), namely the U.S. GPS (Global Positioning System) and the Russian GLONASS (Global’naya Nawigatsionnaya Sputnikowaya Sistema). The data was recorded in the years 2008–2011 by a global network of 92 permanently observing GPS/GLONASS receivers. Two types of daily solutions were generated independently for each GNSS, one including the estimation of geocenter coordinates and one without these parameters. A fair agreement for GPS and GLONASS was found in the geocenter x- and y-coordinate series. Our tests, however, clearly reveal artifacts in the z-component determined with the GLONASS data. Large periodic excursions in the GLONASS geocenter z-coordinates of about 40 cm peak-to-peak are related to the maximum elevation angles of the Sun above/below the orbital planes of the satellite system and thus have a period of about 4 months (third of a year). A detailed analysis revealed that the artifacts are almost uniquely governed by the differences of the estimates of direct solar radiation pressure (SRP) in the two solution series (with and without geocenter estimation). A simple formula is derived, describing the relation between the geocenter z-coordinate and the corresponding parameter of the SRP. The effect can be explained by first-order perturbation theory of celestial mechanics. The theory also predicts a heavy impact on the GNSS-derived geocenter if once-per-revolution SRP parameters are estimated in the direction of the satellite’s solar panel axis. Specific experiments using GPS observations revealed that this is indeed the case. Although the main focus of this article is on GNSS, the theory developed is applicable to all satellite observing techniques. We applied the theory to satellite laser ranging (SLR) solutions using LAGEOS. It turns out that the correlation between geocenter and SRP parameters is not a critical issue for the SLR solutions. The reasons are threefold: The direct SRP is about a factor of 30–40 smaller for typical geodetic SLR satellites than for GNSS satellites, allowing it in most cases to not solve for SRP parameters (ruling out the correlation between these parameters and the geocenter coordinates); the orbital arc length of 7 days (which is typically used in SLR analysis) contains more than 50 revolutions of the LAGEOS satellites as compared to about two revolutions of GNSS satellites for the daily arcs used in GNSS analysis; the orbit geometry is not as critical for LAGEOS as for GNSS satellites, because the elevation angle of the Sun w.r.t. the orbital plane is usually significantly changing over 7 days.

Bare-tether sheath and current: comparison of asymptotic theory and kinetic simulations in stationary plasma

Analytical expressions for current to a cylindrical Langmuir probe at rest in unmagnetized plasma are compared with results from both steady-state Vlasov and particle-in-cell simulations. Probe bias potentials that are much greater than plasma temperature (assumed equal for ions and electrons), as of interest for bare conductive tethers, are considered. At a very high bias, both the electric potential and the attracted-species density exhibit complex radial profiles; in particular, the density exhibits a minimum well within the plasma sheath and a maximum closer to the probe. Excellent agreement is found between analytical and numerical results for values of the probe radiusR close to the maximum radius Rmax for orbital-motion-limited (OML) collection at a particular bias in the following number of profile features: the values and positions of density minimum and maximum, position of sheath boundary, and value of a radius characterizing the no-space-charge behavior of a potential near the high-bias probe. Good agreement between the theory and simulations is also found for parametric laws jointly covering the following three characteristic R ranges: sheath radius versus probe radius and bias for Rmax; density minimum versus probe bias for Rmax; and (weakly bias-dependent) current drop below the OML value versus the probe radius for R > Rmax.

Spectrum of 100-kyr glacial cycle: Orbital inclination, not eccentricity

Spectral analysis of climate data shows a strong narrow peak with period ≈100 kyr, attributed by the Milankovitch theory to changes in the eccentricity of the earth’s orbit. The narrowness of the peak does suggest an astronomical origin; however the shape of the peak is incompatible with both linear and nonlinear models that attribute the cycle to eccentricity or (equivalently) to the envelope of the precession. In contrast, the orbital inclination parameter gives a good match to both the spectrum and bispectrum of the climate data. Extraterrestrial accretion from meteoroids or interplanetary dust is proposed as a mechanism that could link inclination to climate, and experimental tests are described that could prove or disprove this hypothesis.

Transition-state-theory calculations for reactions of O(3P) with halogenated olefins

The potential energy surfaces for the reactions of atomic oxygen in its ground electronic state, O(P-3), with the olefins: CF2=CCl2 and CF2=CF - CF3, have been characterized using ab initio molecular orbital calculations. Geometry optimization and vibrational frequency calculations were performed for reactants, transition states and products at the MP2 and QCISD levels of theory using the 6-31G(d) basis set. This database was then used to calculate the rate constants by means of Transition-State-Theory. To obtain a better reference and to test the reliability of the activation barriers we have also carried out computations using the CCSD(T)(fc)/6-311Gdagger, MP4(SDQ)(fc)/CBSB4 and MP2(fc)/CBSB3 single point energy calculations at both of the above levels of theory, as well as with the composite CBS-RAD procedure ( P. M. Mayer, C. J. Parkinson, D. M. Smith and L. Radom, J. Chem. Phys., 1998, 108, 604) and a modi. cation of this approach, called: CBS-RAD( MP2, MP2). It was found that the kinetic parameters obtained in this work particularly with the CBS-RAD ( MP2, MP2) procedure are in reasonable agreement with the experimental values. For both reactions it is found that the channels leading to the olefin double-bond addition predominates with respect to any other reaction pathway. However, on account of the different substituents in the alkenes we have located, at all levels of theory, two transition states for each reaction. Moreover, we have found that, for the reactions studied, a correlation exists between the activation energies and the electronic structure of the transition states which can explain the influence of the substituent effect on the reactivity of the halo-olefins.

Nickel(II) complexes of whole set of the simple ethylenediaminetetracarboxylate ligands:DFT study of complexes invoking molecular orbital and configurational analysis

The whole set of the nickel(II) complexes with no derivatized edta-type hexadentate ligands has been investigated from their structural and electronic properties. Two more complexes have been prepared in order to complete the whole set: trans(O5)-[Ni(ED3AP)]2- and trans(O5O6)-[Ni(EDA3P)]2- complexes. trans(O5) geometry has been verified crystallographically and trans(O5O6) geometry of the second complex has been predicted by the DFT theory and spectral analysis. Mutual dependance has been established between: the number of the five-membered carboxylate rings, octahedral/tetrahedral deviation of metal-ligand/nitrogen-neighbour-atom angles and charge-transfer energies (CTE) calculated by the Morokuma’s energetic decomposition analysis; energy of the absorption bands and HOMO–LUMO gap.

Complete spin and orbital evolution of close-in bodies using a Maxwell viscoelastic rheology

In this paper, we present a formalism designed to model tidal interaction with a viscoelastic body made of Maxwell material. Our approach remains regular for any spin rate and orientation, and for any orbital configuration including high eccentricities and close encounters. The method is to integrate simultaneously the rotation and the position of the planet as well as its deformation. We provide the equations of motion both in the body frame and in the inertial frame. With this study, we generalize preexisting models to the spatial case and to arbitrary multipole orders using a formalism taken from quantum theory. We also provide the vectorial expression of the secular tidal torque expanded in Fourier series. Applying this model to close-in exoplanets, we observe that if the relaxation time is longer than the revolution period, the phase space of the system is characterized by the presence of several spin-orbit resonances, even in the circular case. As the system evolves, the planet spin can visit different spin-orbit configurations. The obliquity is decreasing along most of these resonances, but we observe a case where the planet tilt is instead growing. These conclusions derived from the secular torque are successfully tested with numerical integrations of the instantaneous equations of motion on HD 80606 b. Our formalism is also well adapted to close-in super-Earths in multiplanet systems which are known to have non-zero mutual inclinations.

Single layer lead iodide : computational exploration of structural, electronic and optical properties, strain induced band modulation and the role of spin-orbital-coupling

Graphitic like layered materials exhibit intriguing electronic structures and thus the search for new types of two-dimensional (2D) monolayer materials is of great interest for developing novel nano-devices. By using density functional theory (DFT) method, here we for the first time investigate the structure, stability, electronic and optical properties of monolayer lead iodide (PbI2). The stability of PbI2 monolayer is first confirmed by phonon dispersion calculation. Compared to the calculation using generalized gradient approximation, screened hybrid functional and spin-orbit coupling effects can not only predicts an accurate bandgap (2.63 eV), but also the correct position of valence and conduction band edges. The biaxial strain can tune its bandgap size in a wide range from 1 eV to 3 eV, which can be understood by the strain induced uniformly change of electric field between Pb and I atomic layer. The calculated imaginary part of the dielectric function of 2D graphene/PbI2 van der Waals type hetero-structure shows significant red shift of absorption edge compared to that of a pure monolayer PbI2. Our findings highlight a new interesting 2D material with potential applications in nanoelectronics and optoelectronics.