We use symmetry analysis to show that the G, C and A-type antiferromagnetic P nma perovskites can exhibit magnetoelectric (ME) responses when a ferroelectric instability is induced with epitaxial strain. Using first-principles calculations we compute the values of the allowed ME response in strained CaMnO3 as a model system. Our results show that large linear and non-linear ME responses are present and can diverge when close to the ferroelectric phase transition. By decomposing the electronic and ionic contributions, we explore the detailed mechanism of the ME response.Interest in magnetoelectric (ME) materials has increased over the last few years because of their cross coupling between the electric polarization and magnetization and their consequent potential for technological applications [1]. However, the search for good MEs is facing difficulties: compounds with the required symmetry (breaking of both the time and space inversion) are uncommon, and when these requirements are met, it is often at low temperatures. In addition, although the magnitude of the response is in principle bounded by the product of the dielectric and magnetic permeabilities (α ≤ √ µ), in practice it tends to be much smaller than this value. One promising direction in the search for improved magnetoelectrics is the exploration of multiferroic materials, since the presence of multiple ferroic orders often presents the desired coupling properties for large ME responses [1][2][3]. Another route is the engineering of artificial heterostructures with specific chemistries and symmetries [4][5][6]. Here we demonstrate from symmetry considerations that the P nma G, C or A-type antiferromagnetic perovskites, which are not multiferroic and do not allow a ME response in their bulk form, can become ME when a polar distortion is induced using thin film heteroepitaxial strain. Then, using first-principles calculations for a model example -P nma CaMnO 3 -we show that particularly large ME responses can be achieved in the vicinity of the ferroelectric phase transition.Technical details-We performed all calculations within density functional theory as implemented in the VASP code [7,8]. Since the purpose of the present study is to provide a model example, we restricted ourselves to the Local Density Approximation (LDA) functional and intentionally avoided studying the U and J dependence of the LDA+U functional [9]. We oriented the P nma unit cell with the longest axis along the b direction, and applied a cubic epitaxial strain on the a and c directions by imposing a = c and relaxing the b direction. To compare with previous studies on strained CaMnO 3 [10], we performed all calculations at the calculated generalised gradient approximation (GGA) PBEsol volumes for each strain (see Supplemental Information). The non-collinear properties and the ME response were well converged at a plane wave cutoff of 550 eV and a 4×2×4 k-point grid. To compute the ME responses, we applied a Zeeman magnetic field to the spins with the spin-orbit coupling included and ex...