Polyunsaturated lipids are an essential component of biological membranes, influencing order and dynamics of lipids, protein-lipid interaction, and membrane transport properties. To gain an atomic level picture of the impact of polyunsaturation on membrane properties, quantum mechanical (QM) and empirical force field based calculations have been undertaken. The QM calculations of the torsional energy surface for rotation about vinyl-methylene bonds reveal low barriers to rotation, indicating an intrinsic propensity toward flexibility. Based on QM and experimental data, empirical force field parameters were developed for polyunsaturated lipids and applied in a 16 ns molecular dynamics (MD) simulation of a 1-stearoyl-2-docosahexaenoyl-sn-glyerco-3-phosphocholine (SDPC) lipid bilayer. The simulation results are in good agreement with experimental data, suggesting an unusually high degree of conformational flexibility of polyunsaturated hydrocarbon chains in membranes. The detailed analysis of chain conformation and dynamics by simulations is aiding the interpretation of experimental data and is useful for understanding the unique role of polyunsaturated lipids in biological membranes. The complete force field is included as Supporting Information and is available from http://www.pharmacy.umaryland.edu/faculty/amackere/research.html.