We have developed a quantum molecular dynamics simulation method which can widely search hydrogen systems from an isolated hydrogen molecule to condensed phase hydrogens exhibiting strong nuclear quantum effects such as zero-point energy and nuclear delocalization. We will briefly summarize derivations and advantages of our method and report our recent results on liquid, solid and supercooled hydrogens, suggesting new molecular pictures and guiding insights useful for future experiments of these mysterious quantum molecules.