A boson two-leg ladder in the presence of a synthetic magnetic flux is investigated by means of bosonization techniques and density matrix renormalization group (DMRG). We follow the quantum phase transition from the commensurate Meissner to the incommensurate vortex phase with increasing flux at different fillings. When the applied flux is ρπ and close to it, where ρ is the filling per rung, we find a second incommensuration in the vortex state that affects physical observables such as the momentum distribution, the rung-rung correlation function and the spin-spin and chargecharge static structure factors.A remarkable characteristic of charged systems with broken U(1) global gauge symmetry such as superconductors is the Meissner-Ochsenfeld effect [1]. In the Meissner phase, below the critical field H c1 , a superconductor behaves as a perfect diamagnet, i.e. it develops surface currents that fully screen the external magnetic field. In a type-II superconductor, for fields above H>H c1 , an Abrikosov vortex lattice phase is formed in the system, where the magnetic field penetrates into vortex cores. In quasi one-dimensional systems, analogues of the Meissner and Abrikosov vortex lattice have been predicted for the bosonic two-leg ladder [2-5], the simplest system where orbital magnetic field effects are allowed. It was shown that in this model, the quantum phase transition between the Meissner and the Vortex phase is a commensurate-incommensurate (C-IC) transition [6][7][8]. For ladder systems at commensurate filling, a chiral Mott insulator phase with currents circulating in loops commensurate with the ladder was obtained [9][10][11][12]. Initially, Josephson junction arrays [13][14][15][16] were proposed as experimental realizations of bosonic one-dimensional systems [17,18]. However, Josephson junctions are dissipative and open systems [19-21] that cannot be described using a Hermitian manybody Hamiltonian in a canonical formalism. Moreover, the quantum effects in the vortex phase of the Josephson ladder are weak [22]. Fortunately, with the recent advent of ultracold atomic gases, another route to realize low dimensional strongly interacting bosonic systems has opened [23][24][25]. Atoms being neutral, it is necessary to find a way to realize an artificial magnetic flux acting on the ladder. Alternatively, one can consider the mapping of the two-leg ladder bosonic model to a two-component spinor boson model in which the bosons in the upper leg become spin-up bosons and the bosons in the lower leg spin-down bosons. Under such mapping, the magnetic flux of the ladder becomes a spin-orbit coupling for the spinor bosons. Theoretical proposals to realize either artificial gauge fields or artificial spin-orbit coupling have been put forward [26,27], and an artificial spinorbit coupling has been achieved in a cold atoms experiment [28,29]. Recently, the Meissner effect and the formation of a vortex state have been observed for non-interacting ultracold bosonic atoms bosons on a two leg ladder in artificial gauge fi...