A detailed study of vortices is presented in Ginzburg-Landau (or Abelian Higgs) models with two complex scalars (order parameters) assuming a general U(1)×U(1) symmetric potential. Particular emphasis is given to the case, when only one of the scalars obtains a vacuum expectation value (VEV). It is found that for a significantly large domain in parameter space vortices with a scalar field condensate in their core (condensate core, CC) coexist with Abrikosov-Nielsen-Olesen (ANO) vortices. Importantly CC vortices are stable and have lower energy than the ANO ones. Magnetic bags or giant vortices of the order of 1000 flux quanta are favoured to form for the range of parameters ("strong couplings") appearing for the superconducting state of liquid metallic hydrogen (LMH). Furthermore, it is argued that finite energy/unit length 1VEV vortices are smoothly connected to fractional flux 2VEV ones. Stable, finite energy CC-type vortices are also exhibited in the case when one of the scalar fields is neutral.In a considerable number of physical theories describing rather different situations, vortices play often an essential rôle to understand key phenomena. In gauge field theories spontaneously broken by scalar fields, the vortex of reference is undoubtedly the celebrated Abrikosov-NielsenOlesen (ANO) one [1] associated to the breaking of a U(1) gauge symmetry by a complex scalar doublet. ANO vortices correspond to the planar cross-sections of static, straight, magnetic flux-tubes, with an SO(2) cylindrical symmetry. Their magnetic flux is quantized as Φ = nΦ 0 , where Φ 0 is an elementary flux unit and n is an integer. The integer n can be identified with a topological invariant, the winding number of the complex scalar, which is also responsible for their remarkable stability. Rotationally symmetric ANO vortex solutions form families labelled by n and by the mass ratio β = m s /m v , where m s , resp. m v denote the mass of the scalar resp. of the vector field.The ubiquity of vortices in different branches of physics, ranging from condensed matter systems, such as superfluids, superconductors [2,3,4] to cosmic strings in high energy physics [5,6] greatly contributes to their importance. By now models of superconductors with several order parameters (scalar doublets) have become the subject of intense theoretical and experimental studies [7,8]. Under extremely high pressure liquid metallic hydrogen (LMH) is expected to undergo a phase transition to a superconducting state, where two types of Cooper pairs are formed, one from electrons and another one from protons [9,10,11,12,13,14]. For experimental data on the existence of liquid metallic hydrogen see Refs. [15,16], for numerical simulations, 1