Convolutional Neutral Networks have been successfully applied in searching for strong lensing systems, leading to discoveries of new candidates from large surveys. On the other hand, systematic investigations about their robustness are still lacking. In this paper, we first construct a neutral network, and apply it to $r$-band images of Luminous Red Galaxies (LRGs) of the Kilo Degree Survey (KiDS) Data Release 3 to search for strong lensing systems. We build two sets of training samples, one fully from simulations, and the other one using the LRG stamps from KiDS observations as the foreground lens images. With the former training sample, we find 48 high probability candidates after human-inspection, and among them, 27 are newly identified. Using the latter training set, about 67\% of the aforementioned 48 candidates are also found, and there are 11 more new strong lensing candidates identified. We then carry out tests on the robustness of the network performance with respect to the variation of PSF. With the testing samples constructed using PSF in the range of 0.4 to 2 times of the median PSF of the training sample, we find that our network performs rather stable, and the degradation is small. We also investigate how the volume of the training set can affect our network performance by varying it from 0.1 millions to 0.8 millions. The output results are rather stable showing that within the considered range, our network performance is not very sensitive to the volume size.

Gravitational and plasma lensing share the same mathematical formalism in the limit of geometrical optics. Both phenomena can be effectively described by a projected, two-dimensional deflection potential whose gradient causes an instantaneous light deflection in a single, thin lens plane. We highlight the differences in the time-delay and lensing equations that occur because plasma lensing is caused by a potential directly proportional to the deflecting electron number density and gravitational lensing is caused by a potential related to the deflecting mass density by a Poisson equation. Since we treat plasma and gravitational lensing as thin-screen effective theories, their degeneracies are both caused by the unknown distribution of deflecting objects. Deriving the formalism-intrinsic degeneracies for plasma lensing, we find that they are analogous to those occurring in gravitational lensing. To break the degeneracies, galaxies and galaxy-cluster scale strong gravitational lenses must rely on additional assumptions or complementary observations. Physically realistic assumptions to arrive at self-consistent lens and source reconstructions can be provided by simulations and analytical effective theories. In plasma lensing, a deeper understanding of the deflecting electron density distributions is still under development, so that a model-based comprehensive lens reconstruction is not yet possible. However, we show that transient lenses and multi-wavelength observations help to break the arising degeneracies. We conclude that the development of an observation-based inference of local lens properties seems currently the best way to further probe the morphologies of plasma electron densities. Due to the simpler evidence-based breaking of the lensing degeneracies, we expect to obtain tighter constraints on the local plasma electron densities than on the gravitationally deflecting masses.

In this paper we investigate the relation between the potential, and the geometrical time delays in the gravitational lensing. In the original paper of Shapiro (1964), it is stated that there is a time delay in the radar signals between Earth and Venus, passing near a massive object (the Sun), relative to a case of the absence of such mass. The reason is connected with the influence of the gravity on the speed of a light rays in a gravitational potential. The contribution from the change of the path length, which happens to be of a second order, is considered as negligible. Nevertheless, in the gravitational lens theory the geometrical delay, related to the change of path length, is routinely taken into account along with the potential term. Inspired by this apparent discrepancy, in this work we address the contribution of the geometric part of the time delay in different situations. We introduce a unified treatment with two limiting regimes of lensing. One of them corresponds to the time delay experiments near the Sun where geometrical delay is shown to be negligible. The second corresponds to the typical gravitational lens scenario with multiple imaging where the geometrical delay is shown to be significant. In particular, it is shown that the geometrical delay is non-negligible in the case of good alignment between source, lens and observer, because in such a case it becomes a first order quantity.

Spherical plasma lens models are known to suffer from a severe over-pressure problem, with some observations requiring lenses with central pressures up to millions of times in excess of the ambient ISM. There are two ways that lens models can solve the over-pressure problem: a confinement mechanism exists to counter the internal pressure of the lens, or the lens has a unique geometry, such that the projected column-density appears large to an observer. This occurs with highly asymmetric models, such as edge-on sheets or filaments, with potentially low volume-density. In the first part of this work we investigate the ability of non-magnetized plasma filaments to mimic the magnification of sources seen behind spherical lenses and we extend a theorem from gravitational lens studies regarding this model degeneracy. We find that for plasma lenses, the theorem produces unphysical charge density distributions. In the second part of the work, we consider the plasma lens over-pressure problem. Using magnetohydrodynamics, we develop a non self-gravitating model filament confined by a helical magnetic field. We use toy models in the force-free limit to illustrate novel lensing properties. Generally, magnetized filaments may act as lenses in any orientation with respect to the observer, with the most high density events produced from filaments with axes near the line of sight. We focus on filaments that are perpendicular to the line of sight that show the toroidal magnetic field component may be observed via the lens rotation measure.

Radio signals are delayed when propagating through plasma. This type of delay is frequency-dependent and is usually used for estimating the projected number density of electrons along the line of sight, called the dispersion measure. The dense and clumpy distribution of plasma can cause refractive deflections of radio signals, analogous to lensing effects. Similar to gravitational lensing, there are two contributions to the time delay effect in plasma lensing: a geometric delay, due to increased path length of the signal, and a dispersive delay due to the change of speed of light in a plasma medium. We show the delay time for two models of the plasma distribution, and point out that the estimated dispersion measure can be biased. Since the contribution of the geometric effect can be comparable to that of the dispersive delay, the bias in the measured dispersion measure can be dramatically large if plasma lensing effects are not taken into account when signals propagate through a high-density gradient clump of plasma.

Plasma lensing is the refraction of low-frequency electromagnetic rays due to free electrons in the interstellar medium. Although the phenomenon has a distinct similarity to gravitational lensing, particularly in its mathematical description, plasma lensing introduces other additional features, such as wavelength dependence, radial rather than tangential image distortions, and strong demagnification of background sources. Axisymmetrical models of plasma lenses have been well-studied in the literature, but density distributions with more complicated shapes can provide new and exotic image configurations and increase the richness of the magnification properties. As a first step towards non-axisymmetrical distributions, we study two families of elliptical plasma lens, softened power-law and exponential plasma distributions. We perform numerical studies on each lens model, and present them over a parameter space. In addition to deriving elliptical plasma lens formulae, we also investigate the number of critical curves that the lens can produce by studying the lens parameter space, in particular the dependence on the lensing ellipticity. We find that the introduction of ellipticity into the plasma distribution can enhance the lensing effects as well as the complexity of the magnification map.

Minkowski functionals (MFs) quantify the topological properties of a given field probing its departure from Gaussianity. We investigate their use on lensing convergence maps in order to see whether they can provide further insights on the underlying cosmology with respect to the standard second-order statistics, i.e., cosmic shear tomography. To this end, we first present a method to match theoretical predictions with measured MFs taking care of the shape noise, imperfections in the map reconstruction, and inaccurate description of the nonlinearities in the matter power spectrum and bispectrum. We validate this method against simulated maps reconstructed from shear fields generated by the MICE simulation. We then perform a Fisher matrix analysis to forecast the accuracy on cosmological parameters from a joint MFs and shear tomography analysis. It turns out that MFs are indeed helpful to break the $\Omega_{\rm m}$--$\sigma_8$ degeneracy thus generating a sort of chain reaction leading to an overall increase of the Figure of Merit.

In contrast to the converging, achromatic behaviour of axisymmetric gravitational lenses, diverging frequency-dependent lensing occurs from refraction due to a distribution of over-dense axisymmetric plasma along an observer's line of sight. Such plasma lenses are particularly interesting from the point of view of astronomical observations because they can both magnify and dim the appearance of background sources as a function of frequency. Plasma lensing is believed to be involved in a number of separate phenomena involving the scintillation of radio pulsars, extreme scattering events of background radio sources and may also play a role in the generation of fast radio bursts. These lensing phenomena are believed to occur in dense environments, in which there may be many density perturbations between an observer and background source. In this work we generalize individual plasma lens models to produce dual component lenses using families of plasma lens models previously studied in the literature, namely the exponential and softened power-law lenses. Similar to binary gravitational lens models, these dual component plasma lenses feature a rich and complex critical and caustic morphology, as well as generate more complicated light curves. We map the number of criticals formed for a given component separation and angular size, and highlight a relevant degeneracy between two particular models. This work provides an argument in favor of close monitoring of extreme scattering events in progress in order to break such model degeneracies.

Accurate shape measurements are essential to infer cosmological parameters from large area weak gravitational lensing studies. The compact diffraction-limited point-spread function (PSF) in space-based observations is greatly beneficial, but its chromaticity for a broad band observation can lead to new subtle effects that could hitherto be ignored: the PSF of a galaxy is no longer uniquely defined and spatial variations in the colours of galaxies result in biases in the inferred lensing signal. Taking Euclid as a reference, we show that this colourgradient bias (CG bias) can be quantified with high accuracy using available multi-colour Hubble Space Telescope (HST) data. In particular we study how noise in the HST observations might impact such measurements and find this to be negligible. We determine the CG bias using HST observations in the F606W and F814W filters and observe a correlation with the colour, in line with expectations, whereas the dependence with redshift is weak. The biases for individual galaxies are generally well below 1%, which may be reduced further using morphological information from the Euclid data. Our results demonstrate that CG bias should not be ignored, but it is possible to determine its amplitude with sufficient precision, so that it will not significantly bias the weak lensing measurements using Euclid data.

The unprecedented quality, the increased dataset, and the wide area of ongoing and near future weak lensing surveys allows to move beyond the standard two points statistics thus making worthwhile to investigate higher order probes. As an interesting step towards this direction, we expolore the use of higher order moments (HOM) of the convergence field as a way to increase the lensing Figure of Merit (FoM). To this end, we rely on simulated convergence to first show that HOM can be measured and calibrated so that it is indeed possible to predict them for a given cosmological model provided suitable nuisance parameters are introduced and then marginalized over. We then forecast the accuracy on cosmological parameters from the use of HOM alone and in combination with standard shear power spectra tomography. It turns out that HOM allow to break some common degeneracies thus significantly boosting the overall FoM. We also qualitatively discuss possible systematics and how they can be dealt with.

The gravitational coupling of a long-wavelength tidal field with small-scale density fluctuations leads to anisotropic distortions of the locally measured small-scale matter correlation function. Since the local correlation function is known to be statistically isotropic in the absence of such tidal interactions, the tidal distortions can be used to reconstruct the long-wavelength tidal field and large-scale density field in analogy with the cosmic microwave background lensing reconstruction. In this paper we present the theoretical framework of cosmic tidal reconstruction and test the reconstruction in numerical simulations. We find that the density field on large scales can be reconstructed with good accuracy and the cross-correlation coefficient between the reconstructed density field and the original density field is greater than 0.9 on large scales $(k\lesssim0.1\ h/\mathrm{Mpc})$, with the filter scale $\sim1.25\ \mathrm{Mpc}/h$. This is useful in the 21cm intensity mapping survey, where the long-wavelength radial modes are lost due to a foreground subtraction process.

In the standard gravitational lensing scenario, rays from a background source are bent in the direction of a foreground lensing mass distribution. Diverging lens behaviour produces deflections in the opposite sense to gravitational lensing, and is also of astrophysical interest. In fact, diverging lensing due to compact distributions of plasma has been proposed as an explanation for the extreme scattering events (ESEs) that produce frequency-dependent dimming of extra-galactic radio sources, and may also be related to the refractive radio-wave phenomena observed to affect the flux density of pulsars. In this work we study the behaviour of two families of astrophysical diverging lenses in the geometric optics limit, the power-law and the exponential plasma lenses. Generally, the members of these model families show distinct behaviour in terms of image formation and magnification, however the inclusion of a finite core for certain power-law lenses can produce a caustic and critical curve morphology that is similar to the well-studied Gaussian plasma lens. Both model families can produce dual radial critical curves, a novel distinction from the tangential distortion usually produced by gravitational (converging) lenses. The deflection angle and magnification of a plasma lens varies with the observational frequency, producing wavelength-dependent magnifications that alter the amplitudes and the shape of the light curves. Thus, multi-wavelength observations can be used to physically constrain the distribution of the electron density in such lenses.

Gravitational lensing has long been considered as a valuable tool to determine the total mass of galaxy clusters. The shear profile as inferred from the statistics of ellipticity of background galaxies allows to probe the cluster intermediate and outer regions thus determining the virial mass estimate. However, the mass sheet degeneracy and the need for a large number of background galaxies motivate the search for alternative tracers which can break the degeneracy among model parameters and hence improve the accuracy of the mass estimate. Lensing flexion, i.e. the third derivative of the lensing potential, has been suggested as a good answer to the above quest since it probes the details of the mass profile. We investigate here whether this is indeed the case considering jointly using weak lensing, magnification and flexion. We use a Fisher matrix analysis to forecast the relative improvement in the mass accuracy for different assumptions on the shear and flexion signal - to - noise (S/N) ratio also varying the cluster mass, redshift, and ellipticity. It turns out that the error on the cluster mass may be reduced up to a factor 2 for reasonable values of the flexion S/N ratio. As a general result, we get that the improvement in mass accuracy is larger for more flattened haloes, but extracting general trends is a difficult because of the many parameters at play. We nevertheless find that flexion is as efficient as magnification to increase the accuracy in both mass and concentration determination.

Large area lensing surveys are expected to make it possible to use cosmic shear tomography as a tool to severely constrain cosmological parameters. To this end, one typically relies on second order statistics such as the two - point correlation fucntion and its Fourier counterpart, the power spectrum. Moving a step forward, we wonder whether and to which extent higher order stastistics can improve the lensing Figure of Merit (FoM). In this first paper of a series, we investigate how second, third and fourth order lensing convergence moments can be measured and use as probe of the underlying cosmological model. We use simulated data and investigate the impact on moments estimate of the map reconstruction procedure, the cosmic variance, and the intrinsic ellipticity noise. We demonstrate that, under realistic assumptions, it is indeed possible to use higher order moments as a further lensing probe.

In gravitational lensing, the magnification effect changes the luminosity and size of a background galaxy. If the image sizes are not small compared to the scale over which the magnification and shear vary, higher-order distortions occur which are termed differential magnification. We give an approximation of the magnification gradient for several halo models. Assuming a symmetric distribution of source brightness, estimates for the differential magnification are obtained and then tested with simulations. One of the main uncertainties of our estimators comes from the finite resolution of the image. We study the strength of our method with the resolution of current and future telescopes. We point out that out method is a potential approach to estimate the first flexion, and can be used to study galaxy and cluster mass profiles.

We study gravitational lensing when plasma surrounds the lens. An extra deflection angle is induced by the plasma in addition to the deflection generated by gravity. An inhomogeneous plasma distribution generates a greater effect than a homogeneous one, and may cause significant effects to be detected in low frequency radio observations (a few hundred MHz). In particular, the lensed image positions will be different between optical and radio observations. The change of position due to a plasma can reach a few tens of milli-arcsec, which is readily detectable. One can use the position difference in different frequencies to estimate the density of plasma in the lens. The magnification ratios between multiple images are mainly determined by other properties of the lens, and are only weakly affected by the plasma. More importantly, we find that the strong lensing time delay will be affected by the plasma. Estimation of the Hubble constant from the time delay in low radio frequency observation may be slightly biased due to plasma in the lens. Unfortunately, the ionosphere of the Earth strongly affects low frequency radio observations. Thus our ability to detect the effect depends on how well we are able to calibrate out the ionosphere.

We study the lensing magnification effect on background galaxies. Differential magnification due to different magnifications of different source regions of a galaxy will change the lensed composite spectra. The derived properties of the background galaxies are therefore biased. For simplicity, we model galaxies as a superposition of an axis-symmetric bulge and a face-on disk in order to study the differential magnification effect on the composite spectra. We find that some properties derived from the spectra (e.g., velocity dispersion, star formation rate and metallicity) are modified. Depending on the relative positions of the source and the lens, the inferred results can be either over- or under-estimates of the true values. In general, for an extended source at strong lensing regions with high magnifications, the inferred physical parameters (e.g. metallicity) can be strongly biased. Therefore detailed lens modelling is necessary to obtain the true properties of the lensed galaxies.

We study the effects on the number counts of sub-millimetre galaxies due to gravitational lensing. We explore the effects on the magnification cross section due to halo density profiles, ellipticity and cosmological parameter (the power-spectrum normalisation $\sigma_8$). We show that the ellipticity does not strongly affect the magnification cross section in gravitational lensing while the halo radial profiles do. Since the baryonic cooling effect is stronger in galaxies than clusters, galactic haloes are more concentrated. In light of this, a new scenario of two halo population model is explored where galaxies are modeled as a singular isothermal sphere profile and clusters as a Navarro, Frenk and White (NFW) profile. We find the transition mass between the two has modest effects on the lensing probability. The cosmological parameter $\sigma_8$ alters the abundance of haloes and therefore affects our results. Compared with other methods, our model is simpler and more realistic. The conclusions of previous works is confirm that gravitational lensing is a natural explanation for the number count excess at the bright end.

In a systematic search over 11 cluster fields from Cluster Lensing And Supernova survey with Hubble (CLASH) we identify ten passively evolving massive galaxies at redshift z~2.We derive the stellar properties of these galaxies using HST WFC3/ACS multiband data, together with Spitzer IRAC observations. We also deduce the optical rest-frame effective radius of these high redshift objects. The derived stellar masses and measured effective radii have been corrected by the lensing magnification factors, which are estimated by simply adopting the spherical NFW model for the foreground cluster lens. The observed near-IR images, obtained by HST WFC3 camera with high spatial resolution and lensed by the foreground clusters, enable us to study the structures of such systems. Nine out of ten galaxies have on average three times smaller effective radius than local ETGs of similar stellar masses, in agreement with previous works at redshift 1.5 < z < 2.5. Combined with literature data for z~2, we find that the mass-normalized effective radius scales with redshift as re/M^0.56 \propto (1 + z)^{-1.13}. We confirm that their size distribution shows a large scatter: from normal size to ~5 times smaller compared to local ETGs with similar stellar masses. The 1-{\sigma} scatter {\sigma}_{logre} of the size distribution is 0.22 and 0.34 at z~1.6 and z~2.1,respectively.The observed large size scatter has to be carefully taken into account in galaxy evolution model predictions.

N-body simulations predict that dark matter halos with different mass scales are described by a universal model, the Navarro-Frenk-White (NFW) density profiles. As a consequence of baryonic cooling effects, the halos will become more concentrated, and similar to an isothermal sphere over large range in radii ($\sim 300$ $h^{-1}$kpc). The singular isothermal sphere model however has to be truncated artificially at large radii since it extends to infinity. We model a massive galaxy halo as a combination of an isothermal sphere and an NFW density profile. We give an approximation for the mass concentration at different baryon fractions and present exact expressions for the weak lensing shear and flexion for such a halo. We compare the lensing properties with a Singular Isothermal Sphere and NFW profiles. We find that the combined profile can generate higher order lensing signals at small radii and is more efficient in generating strong lensing events. In order to distinguish such a halo profile from the SIS or NFW profiles, one needs to combine strong and weak lensing constraints on small and large radii.

Embedded clusters are ideal laboratories to understand the early phase of the dynamical evolution of clusters as well as the massive star formation. An interesting observational phenomenon is that some of the embedded clusters show mass segregation, i.e., the most massive stars are preferentially found near the cluster center. In this paper, we develop a new approach to describe mass segregation. Using this approach and the Two Micron All Sky Survey Point Source Catalog (2MASS PSC), we analyze eighteen embedded clusters in the Galaxy. We find that eleven of them are mass-segregated and that the others are non-mass-segregated. No inversely mass-segregated cluster is found.

Estimating the ellipticity of dark matter haloes at the galaxy or galaxy-cluster scale can provide important constraints on the formation of galaxies or clusters, as well as on the nature of dark matter. We show in this paper that the spin-3 gravitational flexion can add useful information on the ellipticity of lensing haloes. We introduce a general formalism to decompose fields with arbitrary spin into radial and tangential components. The ratio of the tangential and radial flexion components directly estimates the lens ellipticity. We point out that any centroid offset will significantly bias our estimate, which on the other hand can be used to determine the centre of the lens halo.

We constrain properties of cluster haloes by performing likelihood analysis using lensing shear and flexion data. We test our analysis using two mock cluster haloes: an isothermal ellipsoid (SIE) model and a more realistic elliptical Navarro-Frenk-White (eNFW) model. For both haloes, we find that flexion is more sensitive to the halo ellipticity than shear. The introduction of flexion information significantly improves the constraints on halo ellipticity, orientation and mass. We also point out that there is a degeneracy between the mass and the ellipticity of SIE models in the lensing signal.

The cold dark matter theory predicts triaxial dark matter haloes. The radial distribution of halo ellipticity depends on baryonic processes and the nature of dark matter particles (collisionless or collisional). Here we show that we can use lensing flexion ratios to measure the halo ellipticity as a function of radius. We introduce a weight function and study the relationship between the first and second order statistics of flexion ratios, both of which can be used to reduce the bias in the estimate of ellipticity. we perform numerical tests for our method, and demonstrate that it can reduce the bias and determine the halo ellipticity as a function of radius. We also point out that the minimum mean flexion ratio can be used to trace the centres of galaxy clusters.

Aims. The predictions of the ellipticity of the dark matter halos from models of structure formation are notoriously difficult to test with observations. A direct measurement would give important constraints on the formation of galaxies, and its effect on the dark matter distribution in their halos. Here we show that galaxy-galaxy flexion provides a direct and potentially powerful method for determining the ellipticity of (an ensemble of) elliptical lenses. Methods. We decompose the spin-1 flexion into a radial and a tangential component. Using the ratio of tangential-to- radial flexion, which is independent of the radial mass profile, the mass ellipticity can be estimated. Results. An estimator for the ellipticity of the mass distribution is derived and tested with simulations. We show that the estimator is slightly biased. We quantify this bias, and provide a method to reduce it. Furthermore, a parametric fitting of the flexion ratio and orientation provides another estimate for the dark halo ellipticity, which is more accurate for individual lenses Overall, galaxy-galaxy flexion appears as a powerful tool for constraining the ellipticity of mass distributions.

Galaxy clusters are considered as excellent probes for cosmology. For that purpose, their mass needs to be measured and their structural properties needs to be understood. We propose a method for galaxy cluster mass reconstruction which combines information from strong lensing, weak lensing shear and flexion. We extend the weak lensing analysis to the inner parts of the cluster and, in particular, improve the resolution of substructure. We use simulations to show that the method recovers the mass density profiles of the cluster. We find that the weak lensing flexion is sensitive to substructure. After combining the flexion data into the joint weak and strong lensing analysis, we can resolve the cluster properties with substructures.

An exact analytic calculation is presented for the spectrum of relic gravitational waves in the scenario of accelerating Universe $\Omega_{\Lambda}+\Omega_m = 1$. The spectrum formula contains explicitly the parameters of acceleration, inflation, reheating, and the (tensor/scalar) ratio, so that it can be employed for a variety of cosmological models. We find that the spectrum depends on the behavior of the present accelerating expansion. The amplitude of gravitational waves for the model $\Omega_{\Lambda}=0.65$ is about $\sim 50 %$ greater than that of the model $\Omega_{\Lambda}=0.7$, an effect accessible to the designed sensitivities of LIGO and LISA. The spectrum sensitively depends on the inflationary models with $a(\tau) \propto |\tau|^{1+\beta}$, and a larger $\beta$ yields a flatter spectrum, producing more power. The current LIGO results rule out the inflationary models of $\beta \geq -1.8$. The LIGO with its design sensitivity and the LISA will also be able to test the model of $\beta=-1.9$. We also examine the constraints on the spectral energy density of relic gravitational waves. Both the LIGO bound and the nucleosynthesis bound point out that the model $\beta=-1.8$ is ruled out, but the model $\beta=-2.0$ is still alive. The exact analytic results also confirm the approximate spectrum and the numerical one in our previous work.

In this paper we briefly present our works on the relic gravitational waves (RGW) and the CMB polarization in the accelerating universe. The spectrum of RGW has been obtained, showing the influence of the dark energy. Compared with those from non-accelerating models, the shape of the spectrum is approximately similar, nevertheless, the amplitude of RGW now acquires a suppressing factor of the ratio of matter over dark energy $\propto \Omega_m/\Omega_{\Lambda}\sim 0.4$ over almost the whole range of frequencies. The RGW spectrum is then used as the source to calculate the spectra of CMB polarization. By a two half Gaussian function as an approximation to the visibility function during the photon decoupling, both the "electric" and "magnetic" spectra have been analytically derived, which are quite close to the numerical ones. Several physical elements that affect the spectra have been examined, such as the decoupling process, the inflation, the dark energy, the baryons, etc.

In weak gravitational lensing, the image distortion caused by shear measures the projected tidal gravitational field of the deflecting mass distribution. To lowest order, the shear is proportional to the mean image ellipticity. If the image sizes are not small compared to the scale over which the shear varies, higher-order distortions occur, called flexion. For ordinary weak lensing, the observable quantity is not the shear, but the reduced shear, owing to the mass-sheet degeneracy. Likewise, the flexion itself is unobservable. Rather, higher-order image distortions measure the reduced flexion, i.e., derivatives of the reduced shear. We derive the corresponding lens equation in terms of the reduced flexion and calculate the resulting relation between brightness moments of source and image. Assuming an isotropic distribution of source orientations, estimates for the reduced shear and flexion are obtained; these are then tested with simulations. In particular, the presence of flexion affects the determination of the reduced shear. The results of these simulations yield the amount of bias of the estimators, as a function of the shear and flexion. We point out and quantify a fundamental limitation of the flexion formalism, in terms of the product of reduced flexion and source size. If this product increases above the derived threshold, multiple images of the source are formed locally, and the formalism breaks down. Finally, we show how a general (reduced) flexion field can be decomposed into its four components: two of them are due to a shear field, carrying an E- and B-mode in general. The other two components do not correspond to a shear field; they can also be split up into corresponding E- and B-modes.