In this work, noncovalent interactions including hydrogen bonds, C···C, N···O, and van der Waals forces between paracetamol and formaldehyde were investigated using the second-order perturbation theory MP2 in conjunction with the correlation consistent basis sets (aug-cc-pVDZ and aug-cc-pVTZ). Two molecular conformations of paracetamol were considered. Seven equilibrium geometries of dimers were found from the result of the interactions with formaldehyde for each conformation of paracetamol. Interaction energies of complexes with both ZPE and BSSE corrections range from -7.0 to -21.7 kJ mol-1. Topological parameters (such as electron density, its Laplacian, and local electron energy density at the bond critical points) of the bonds from atoms in molecules theory were analyzed in detail. The natural bond orbital analysis showed that the stability of complexes was controlled by noncovalent interactions including O-H···O, N-H···O, C-H···O, C-H···N, C-H···H-C, C···C, and N···O. The red- and blue-shifted hydrogen bonds could both be observed in these complexes. The properties of these interactions were also further examined in water using a polarized continuum model. In water, the stability of the complex was slightly reduced as compared to that in the gas phase.