The Astropy project supports and fosters the development of open-source and openly-developed Python packages that provide commonly-needed functionality to the astronomical community. A key element of the Astropy project is the core package Astropy, which serves as the foundation for more specialized projects and packages. In this article, we provide an overview of the organization of the Astropy project and summarize key features in the core package as of the recent major release, version 2.0. We then describe the project infrastructure designed to facilitate and support development for a broader ecosystem of inter-operable packages. We conclude with a future outlook of planned new features and directions for the broader Astropy project.

Through synthetic observations of a hydrodynamical simulation of an evolving star-forming region, we assess how the choice of observational techniques affects the measurements of properties which trace star formation. Testing and calibrating observational measurements requires synthetic observations which are as realistic as possible. In this part of the paper series (Paper I), we explore different techniques for how to map the distributions of densities and temperatures from the particle-based simulations onto a Voronoi mesh suitable for radiative transfer and consequently explore their accuracy. We further test different ways to set up the radiative transfer in order to produce realistic synthetic observations. We give a detailed description of all methods and ultimately recommend techniques. We have found that the flux around 20 microns is strongly overestimated when blindly coupling the dust radiative transfer temperature with the hydrodynamical gas temperature. We find that when instead assuming a constant background dust temperature in addition to the radiative transfer heating, the recovered flux is consistent with actual observations. We present around 5800 realistic synthetic observations for Spitzer and Herschel bands, at different evolutionary time-steps, distances and orientations. In the upcoming papers of this series (Paper II, Paper III and Paper IV), we will test and calibrate measurements of the star-formation rate (SFR), gas mass and the star-formation efficiency (SFE) using our realistic synthetic observations.

Weierstrass elliptic and related functions have been recently shown to enable analytical explicit solutions to classical problems in astrodynamics. These include the constant radial acceleration problem, the Stark problem and the two-fixed center (or Euler's) problem. In this paper we review the basic technique that allows for these results and we discuss the limits and merits of the approach. Applications to interplanetary trajectory design are then discussed including low-thrust planetary fly-bys and the motion of an artificial satellite under the influence of an oblate primary including $J_2$ and $J_3$ harmonics.

We present for the first time an explicit, complete and closed-form solution to the three-dimensional problem of two fixed centres, based on Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our solution is exact, valid for all initial conditions and physical parameters of the system (including unbounded orbits and repulsive forces), and expressed via a unique set of formulae. Various properties of the three-dimensional problem of two fixed centres are investigated and analysed, with a particular emphasis on quasi-periodic and periodic orbits, regions of motion and equilibrium points.

While the constant radial acceleration problem is known to be integrable and has received some recent attention in an orbital mechanics context, a closed form explicit solution, relating the state variables to a time parameter, has eluded all researchers so far. It is here shown how such a solution exists and is elegantly expressed in terms of the Weierstrass elliptic and related functions. Previously known facts can be derived from the new explicit solution and new insights are revealed.

We present a new general, complete closed-form solution of the Stark problem in terms of Weierstrass elliptic and related functions. With respect to previous treatments of the problem, our analysis is exact and valid for all values of the external force field, and it is expressed via unique formul\ae{} valid for all initial conditions and parameters of the system. The simple form of the solution allows us to perform a thorough investigation of the properties of the dynamical system, including the identification of quasi-periodic and periodic orbits, the formulation of a simple analytical criterion to determine the boundness of the trajectory, and the characterisation of the equilibrium points.

We revisit the relativistic restricted two-body problem with spin employing a perturbation scheme based on Lie series. Starting from a post-Newtonian expansion of the field equations, we develop a first-order secular theory that reproduces well-known relativistic effects such as the precession of the pericentre and the Lense-Thirring and geodetic effects. Additionally, our theory takes into full account the complex interplay between the various relativistic effects, and provides a new explicit solution of the averaged equations of motion in terms of elliptic functions. Our analysis reveals the presence of particular configurations for which non-periodical behaviour can arise. The application of our results to real astrodynamical systems (such as Mercury-like and pulsar planets) highlights the contribution of relativistic effects to the long-term evolution of the spin and orbit of the secondary body.